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SUMMARY 


This report presents results of an investigation concerning the 
analysis, design, and testing of short curved wall axially symmetrical 
subsonic diffusers utilizing suction through slots in the diffuser walls 
to prevent flow separation. 

A potential flow analysis was made and a digital computer program 
was written for determining the diffuser wall geometry for prescribed 
boundary conditions. The prescribed boundary conditions included the 
following for the fluid velocity at the diffuser walls (both inner and 
outer walls for the annular diffusers): a region of nearly constant high 

velocity from the diffuser entrance to the slot; a region of rapid deceleration 
at the slot, and a region of nearly constant low velocity downstream from 
the slot to the diffuser exit. Branch flow was included in the flow field 
and this allowed the suction slot geometry to be a part of the design. 

Four diffuser geometries were designed, fabricated and tested: 

(1) a 7“inch (18 cm) exit diameter bell shaped diffuser with an area 

ratio of 3:1 and a length to inlet radius ratio of 3*5:1* 

(2) a 9“inch (23 cm) exit diameter annular diffuser with an area ratio 

of 3:1 and a length to inlet width ratio of 5*0:1. 

(3) a 9“inch (23 cm) exit diameter annular diffuser with an area ratio 

of 2.5:1 and a length to inlet width ratio of 3*9:1* 

(4) an 18-inch (46 cm) exit diameter annular diffuser with an area ratio 
of 3:1 and a length to inlet width ratio of 3*6:1* 

For inlet air velocities in the 200-250 ft/sec (61-76 m/sec) range, 
metastable operation required 6.5% suction for the bell shaped diffuser while 
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stable operation required 12% suction. Corresponding suction rates for the 
9“inch (23 cm) annular diffuser (3:1 area ratio) were 8.5% and 16.0% for 
metastable and stable operation respectively. For the 9-inch (23 cm) annular 
diffuser with 2.5:1 area ratio, 6 . 3 % and 12% suction rates were required. For 
the 18-inch (46 cm) annular diffuser, only 7% suction was required for meta- 
stable operation and 10% for stable operation. 

Suction rate requirements were found to increase as inlet velocity was 
decreased. Thus for stable operation of the ]8-inch (46 cm) diffuser the 
required suction rate was about 20% at 100 ft/sec (30 m/sec) inlet velocity. 

In all cases the diffuser effectiveness was above 95% based on the 
conventional definition, and from 81% to 94% when the suction loss was 
accounted for. The exit velocity profiles were virtually flat with no 
more than - 9% variation over 95% of the exit area when operated with 
sufficient suction to prevent separation. 

Non-uniform inlet tests and tests with periodic shedding of wakes upstream 
of the diffuser inlet were conducted on the 9~inch (23 cm) annular diffuser 
with AR=3:1* Non-uniform inlet conditions caused no significant increase 
in the required suction rates. However, periodic shedding of wakes upstream 
of the diffuser inlet caused the suction requirement for stable operation 
to be increased from 16.0% to 18 . 0 %. 

A downstream obstruction in the form of a perforated metal sheet at 
the diffuser exit slightly increased the minimum suction rate requirement at 
low inlet velocities and had no significant effect upon the suction rate 
requirement at inlet velocities above 150 ft/sec (46 m/sec). 
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SECTION I 
INTRODUCTION 


1.1 Background 

In an internal flow system, a diffuser is a transitional section 
which connects a flow passage having a smaller cross-sectional area to 
a flow passage of larger area. Two desirable characteristics of a 
diffuser are that it will provide a nearly uniform velocity distribution 
at its exit plane and that it will transform most of the kinetic energy 
of the fluid at entrance into potential energy in the form of a higher 
static pressure at the exit plane. 

A previous report by the authors [1]* concerned the feasibility of 
designing short two-dimensional diffusers utilizing suction through the 
diffuser wall slots. The design philosophy was to prescribe the diffuser 
wall velocity so as to have essentially constant high velocity upstream 
of a slot, constant low velocity downstream of the slot, and thus have 
all the deceleration occur across the slot. By removing through the slot 
those fluid particles having insufficient kinetic energy to get across 
the adverse pressure gradient in the region of the slot, flow separation 
could be avoided. The concept of having deceleration occur abruptly and 
applying slot suction to this narrow region was originally suggested by 
A. A. Griffith as reported by Lachmann [2] and is henceforth referred to 
the Griffith diffuser. 

The two-dimensional Griffith diffusers designed and tested were only 
partially successful. When sufficient suction was applied to prevent 

' Numbers in brackets refer to the references at the end of this report. 



separation, the diffuser exit plane velocity distribution was nearly uniform 
and very low total pressure losses were achieved. However, the suction rates 
required for stable operation varied from 16% for 2:1 area ratio to 

for A:1 area ratio. Problems were encountered in designing the optimum 
slot geometry and also in controlling the sidewall boundary layer. 

In 1952, L. R. Manoni [ 3 ] designed several bell channel diffusers 
based upon the Griffith concept. He used an electrical analogy tank to 
determine the wall geometry required for the specified wall velocity dis- 
tribution. While his design method was only approximate and somewhat 
cumbersome, he was able to demonstrate that such a concept would enable the 
design and operation of a relatively short diffuser without separation 
provided that sufficient slot suction could be applied. 

In order to design planar flow channel geometries, J. D. Stanitz [A] 
in 1953 performed a double transformation of the describing equations into 
a coordinate system where the geometry did not need to be known for the 
equations to be solved. After the equations were solved on the transformed 
plane, the inverse transformation yielded the desired geometry. 

In 1971, C. D. Nelson [ 5 ] extended the work of Stanitz by developing a 
method for transforming the flow equations of an axially symmetrical 
coordinate system. Nelson also devised an approximation to account for the 
branch flow due to the slot by solving the channel design problem with a 
more complicated set of boundary conditions. 
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w 1.2 Objective 

The objective of the research effort described in this report was to 
examine the concept and the practicality of a short axially symmetrical 
subsonic curved wall diffuser using suction through slots in the diffuser 
walls to prevent flow separation. 

1 . 3 Scope 

In order to achieve the objective a research program was carried out to 

(1) Develop a digital computer program which will determine the wall 
shape of an axially symmetrical or a two-dimensional diffuser for prescribed 
boundary conditions. 

(2) Utilize the program to design one bell-shaped diffuser and three 
annular diffusers with area ratios of about 3:1. 

(3) Fabricate and test each of the above diffusers to determine 

(a) The minimum suction requirement for unseparated flow as a 
function of inlet air velocity; 

(b) Performance characteristics such as diffuser effectiveness, 
total pressure loss, and exit plane velocity distribution; 

(c) Effect of non-uniform inlet velocity distributions; 

(d) Effect of periodic shedding of wakes upstream of the diffuser 
inlet ; 

(e) Effect of imbalance between inner and outer wall suction rates 
per unit slot length; 

(f) Effect of downstream blockage in the form of perforated steel 
plates. 
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SECTION 2 

ANALYSIS AND DESIGN 


2.1 Analysis for Determination of Diffuser Geometry 

The describing equations of incompressible irrotational motion and 

continuity in the axially symmetrical coordinate system are 

2 2 

tQ _ Q = 0 

9R 9X 3X 8R 


a <p , a 2 ^ j_ a 

ax 2 aR 2 R 3R 


( 2 ) 


The geometry of the channel being designed might appear as shown 
Figure 1 . 



n 


Figure 1. Possible flow channel geometry on the 
physical plane. 

Transforming equations (l) and (2) into natural coordinates by using 

ds = cos a dX + sin a dR (3) 

dn = -sin a dX + cos a dR (b) 

and introducing the stream function and velocity potential in natural coordinates 

df = R Q dn (5) 

d<J> = Q ds (6) 
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results in the equations 




9 In Q 
9n 

9a f\ 

- 97 = 0 

(7) 

and 






3 In Q 
9s 

1 9R 9o n 

R 9s 9n ° 

(8) 

The 

channel geometry on the 

natural coordinate plane might then appear 

as shown 

in Figure 2. 

i 

n 





<p=0 


<j>=4> 




ip=0 



Figure 2. 

Possible flow channel geometry 
coordinate system. 

on natural 

Now 

using equations 

(5) and 

(6) to transform from natural coordinates 

to a <t>-¥ 

coordinate system results in the following equations of irrotational 

motion and continuity. 





R 

9 In Q 

3 jp 

- 2* = o 

9<j» 

(9) 

and 

1 9 
R 

In Q 

9<j> 

1 9R 3o ' q 

R 2 9* dip 

(10) 


The flow field on the (f)-'!' plane might appear as shown in Figure 3- 





■ rv s i n a , , . cos a , lM 

dR = — d<|> + ( 1 *0 

Since in eqn. (11) the radial coordinate R remains as a dependent 
variable, one method of solution consists of assuming values for the radial 
coordinates and then alternately solving for In Q and R in an iterative 
manner. 


As a first trial, assume R—l which reduces equation (11) to 

3 2 In Q 3 2 In Q _ 

2 2 - U 

3<T 3<(> 


( 15 ) 
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which is the same as the corresponding plane flow equation of Stanitz [4]. The 
term In Q may be found for a given set of boundary conditions by utilizing some 
numerical scheme such as the Gauss-Siedel method. Boundary conditions are 
prescribed along the wall streamlines and along the inlet and outlet 
equ i potential lines. 

The first approximation to the velocity vector angles a are then com- 
puted from 


a 


2 


- j 


3 In 


d<j> 




f Lin q 

J 3<J> 

<t> 


d’/' 


( 16 ) 


where the path of integration is chosen to be along equi potent i al lines and 
streaml i nes . 

The R coordinates are then computed using the line integrals 


and 


2 

-2 



sin a 

Q. 


dip 


cos a 

Q 


dV 


( 17 ) 

( 18 ) 


These R coordinates are treated as a second approximation to solve 
equation (ll) again for new velocities. New flow angles are computed using 
the line integrals 


3 In Q 

3iJ> 


dip 


, r (] 3 In Q . 1 3R ) 

and V a, = -J — -j 


d^ 


(19) 

( 20 ) 


Corrected R coordinates are again computed from equations (17) and ( 1 8) and the 
entire process is reiterated until the R coordinates have converged sufficiently. 
The axial coordinates are computed next by using the line integrals 


cos a 


( 21 ) 


X 

and X 


2 

2 



Q 
s i n 


R Q 


d<f> 


d\fi 


( 22 ) 


In order to account for the branch flow due to slot suction, a more 
complicated set of boundary conditions are necessary. The physical flow 
system may now appear as in Figure 4. The formerly simple channel wall now 
involves curves B, C, D, and E with some of the flow leaving across C. This 
introduces a stagnation point at A. 



Figure 4. Possible geometry with branch flow. 

On streamlines B, D, and E, the value of the desired velocity is prescribed 
as usual, but on line C. (which is chosen to be a constant potential line), 


3 In Q _ 

3<J> 


= 0 


. aR 

and 3* 


0 


are specified which in essence prescribes a parallel flow across C. 

The stagnation point A introduces a singularity into equation (11) since 
In Q is not defined when Q = 0. In order to circumvent the difficulty introduced 
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by this singularity, a portion of the potential flow solution for two- 
dimensional flow around a wedge was patched into a small region around the 
stagnation point. The introduction of this two-dimensional wedge flow patch 
yields only an approximate solution which becomes more accurate as the radial 
dimension of the patch becomes small relative to the distance of the stagna- 
tion point from the axis of symmetry. 

2.2 Digital Computer Design Program 

The digital computer program which was used to design the wall contours 
of the Griffith diffusers described in this report has the following features. 

(1) It may be used to design either planar (two-dimensional) or axially 
symmetrical (bell shaped or annular) flow channel^. 

(2) Branch flow due to slot suction is incorporated into the design 
program. 

(3) The boundary conditions are flexible; i.e., at either entrance 

to or exit from the flow channel, the velocities may be prescribed 
(Dirichlet boundary conditions) or the velocity gradients may be 
set to zero (Neumann boundary conditions) so as to result in 
unidirectional flow at either entrance or exit. The wall velocities 
are prescribed at up to 1A8 points along each channel wall (both 
inner and outer walls if annular). 

(A). Up to. 1A8 mesh points may be used, in the phi direction and up to 
A8 points may be used in the psi direction. 

(5) Printed output, includes the coordinates (R and X), velocit ies (Q), 
and flow angles (a) for each mesh point and a summary of the wall 
coordinates of the channel design. 

(6) Graphical output includes all the stream! i nes and potential flow 
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lines on the X-R (or X-Y) plane. Also, the inlet and outlet velocity 
profiles are plotted as a function of the R (or Y) coordinate. 

Using a 48 x 135 grid, a typical axially symmetrical design requires about 
45 minutes on an IBM System 3&0 Model 50 computer. A typical two-di mensi onal 
design of comparable grid points requires about 20 minutes. 

The program listing is given in Appendix A. Appendix B describes briefly 
the procedure used in obtaining the 18- inch (46 cm) ..annular diffuser design. 

2.3 Griffith Diffuser Geometry 

Figures 5, 6, and 7 show the geometries of the bell shaped diffuser, 
the 9“inch (23 cm) annular diffuser, and the 18-inch (46 cm) annular dif- 
fuser respectively. The diffusers were designed to have an area ratio of 
about 3:1 since this represents a typical diffuser application. 

The diffuser geometries were analyzed using the Douglas Neumann technique 
[6] to predict the velocity distributions along the diffuser walls and at the 
diffuser inlet and exit planes. A comparison was made between the velocities 

A 

predicted by the analysis program and those prescribed to the design program. 

In each case the agreement was good although slight disagreements occurred in 
the vicinity of the slot. 

The bell shaped diffuser had the simplest geometry and was designed and 
tested first. Only one curved wall had to be considered and its slot was 
relatively far away from the axis of symmetry. 

The second geometry designed and tested was the 9~inch (23 cm) annular 
diffuser. Since two slotted walls were involved with the inner wall slot being 
closer to the axis of symmetry, this case provided a more severe test of the 
slot design procedure using the two-dimensional wedge flow patch. 

The third diffuser used the same outer body as the second diffuser, but a 
smaller inner body was designed so as to change the area ratio from 3:1 to 
1*4 





2.5:1. This diffuser was to demonstrate that a common outer body could be 
used with more than one inner body thus providing more than one area ratio 
without a completely new diffuser. 

The last diffuser designed and tested was the 1 8- i nch (46 cm) annular 
diffuser. This diffuser was to verify that a reduction of slot suction can be 
achieved by reducing the ratio of 6/W. 6 is the boundary layer displacement 

thickness just upstream of the slot, and W is the annular width of the diffuser 
at inlet. 

SECTION 3 

INSTRUMENTATION AND APPARATUS 

3. 1 Test Faci 1 i ty 

Two primary air supplies were used in testing the diffusers. The 18" 

3 

annular diffuser was tested using a 30,000 CFM (14 nrr/sec) industrial fan 
delivering air through a 42-inch (106 cm) by ,48- i nch (122 cm) duct 30 ft. 

(9.1 m) long. 

All other diffusers were tested using a 10,000 CFM {k.72 m /sec) fan 

delivering air through a 20- i nch (51 cm) circular duct 20 ft. (6.1 m) long. 

Flow straightener tubes and fine mesh screens were installed inside the duct 

to provide a uniform steady flow to the test section. A shutter type damper 

at the fan inlet was used to regulate the flow rate. 

Two diesel superchargers were installed in parallel providing for up to 
3 

900 CFM (0.^3 m /sec) suction air through the diffuser slots. Seven gate 
valves located in six suction lines and an air bleed line provided a flexible 
setup for varying both the total amount of slot suction and the distribution of 
suction between the inner wall and the outer wall of the annular diffusers. 
Plastic pipe lengths and fittings were attached between the diffuser and the 
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suction control valves and all joints were sealed with tape to eliminate air 
^leakage. 

The diffusers were fabricated from aluminum stock. The wall contours were 
cut on a pattern lathe. Figure 8 shows the 18-inch (46 cm) annular diffuser 
attached to the end of the large a i r supp 1 y due t . Figure 9 shows the bell 
shaped diffuser attached to the small' air supply duct. 

3.2 Instrumentation and Measurements • ; 

The stagnation pressure at the end of each duct and the static pressure 
at the diffuser inlet were measured with a micromanometer having resolution of 
0.001 inch of water. WaTl static pressure taps were installed in the diffusers 
by glueing 1/1 6- i n ch o. d. brass tubes /into holes drilled in the diffuser walls 
and smoothing the ends to conform to the d i f fuser contour . The pressure signals 
were monitored on an inductive type 'pressure transducer through a rapid scan- 
ning mechanism. Up to 48 ^pressure taps could be scanned in 30' seconds. The 
output of the pressure transducer was , recorded/ on paper tape with a digital 
recorder. •*. /* r , 

Velocity profiles were measured' wi th a : constant temperature hot wire probe. 
The voltage output of the probe, was’ displayed by a multirange digital voltmeter 
with a selective damping device. The probe was positioned in the diffuser 
with a traversing mechan ism ; adapted from a milling machine-bed. 

Depending upon the size of diffuser being tested either five or six flow 
meters were used to measure the slot suction rates. A manometer bank with 20 
legs each 60-inches high employ ing ei/ther water or carbon tetrachloride was used 
to obtain the flow meter pressure differentials. A copper-cons tantan thermo- 
couple was used to measure the. air temperature, and a well-type barometer was 
used to obtain atmospheric pressure. 
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Figure 8. 1 8" (46 cm) annular diffuser attached to the large air supply duct 
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Figure 9. Bell shaped diffuser attached to the small air supply duct 





SECTION 4 


TEST CONDITIONS AND PROCEDURE 

4.1 Test Cond i t i ons 

As reported' in reference [l], the Griffith diffuser operates with either 
very high effectiveness (when unseparated flow' exists) or very low effective- 
ness (when separated flow exists). The minimum suction rate required for 
unseparated flow is thus . a very important diffuser characteristic. 

Experimental' work showed that two modes of operation with unseparated 
flow can exist. In. the stable mode a flow disturbance may cause temporary 
flow separation, but r unseparated flow will return when the disturbance is 
removed. In the metastab le mode, unseparated flow persists even if a dis- 
turbance suf fi c ientf to cause separation is removed.. Such factors as inlet 
velocity, area ratio, suet i on s lot width, and d i f fuser downstream cond i t i ons 
have an effect upon' the minimum suction required for both metastable and for 
stable unseparated flow; - - 

An extens i ve se t of -tests were made, to determi ne vthe minimum suction 
requirements for unseparated flow for each diffuser. Inlet air velocities 
were varied from 50/. ft /sec (15.2 m/sec):-: to 280 ft/sec (85 m/sec) . For the 
9-inch (23 cm) annular diffuser,, slot widths of 1/16-inch (O'. ,1 59 cm) to 1/4-inch 
(0.635 cm) were used. For the 18-inch (46. cm) annular diffuser, slot widths of 
1/8-inch (0.317 cm) to . 1/2-inch (1.27 cm) were used. Area ratios of 2.5:1 and 
3:1 were used . 

The effect of downstream conditions was investigated in two ways. 

(l) A perforated sheet metal plate with 51% free area was attached to 
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the diffuser exit plane of the 9“ inch (23 cm) and of the 18-inch 
(46 cm) annular diffuser. The 1 8- i nch (46 cm) diffuser plate had 
5/8-inch (1.59 cm) holes while the 9" i nch (23 cm) diffuser plate had 
5/l6-inch (.79 cm) holes. 

(2) The usual annular diffuser requires struts which may cause flow 
separation. Three struts located 120° apart were mounted in the 
downstream constant velocity zone of the 9~inch (23 cm) annular 
diffuser (AR = 2.5:1) to s imulate this structural problem. 

The effect of. an unsteady inlet velocity was examined by installing 
a motor driven propeller just upstream of the inlet to the 9“inch (23 cm) 
diffuser. The effect upon suction requirements when driven at 500 RPM 
was exami ned. 

In addition to tests made to determine minimum suction requirements, 
performance tests were made to determi ne di f fuser effectiveness, total 
pressure loss, wall pressure distribution, exit plane velocity d i str i but i on ^ 
and center line velocity distribution for the four d if fuse rs when operat i ng 
with about 200 ft/sec (60 m/sec) inlet velocity and with a suction rate 
sufficient to maintain unseparated flow. 

4.2 Test Procedure 

In determining the minimum suction requirement for stable operation, 
the pressure differential between the end of the supply duct and the end 
of the diffuser inlet section was observed in order to adjust the damper so ^ 
as to give the desired inlet velocity. The suction blowers were started 
with the air bleed valve wide open and the slot suction valves nearly closed, 
thus having a deficiency of suction. The suction valves were slowly 


21 



opened while observing the venturi meter manometers so as to keep the proper 
ratio of outer wall to inner wall suction rate. This ratio was approximately 

the same as the slot length ratio. The air bleed valve was slowly closed so 

\ 

as to increase the suction capacity either simultaneously with or subsequent 
to the opening of the suction valves. Increasing the suction rate sufficiently 
would cause the flow to become unseparated, and this condition was considered 
to be the minimum for stable operation. 

In determining the minimum suction requirement for metastable operation, 
the suction rate was decreased from the amount required for stable operation 
(or from any excess suction rate) until separation occurred. It was again 
necessary to observe the flow manometers as the suction rate was reduced so 
as to keep the suction rate per unit length of slot about the same on 
the inner and outer walls. 

For tests involving performance characteristics such as effectiveness, 
velocity distribution, and wall pressure distribution, three static pressure 
measurements were made with the micromanometer: (l) at the end of the air 

supply duct (2) at the diffuser inlet (3) differential between (l) and 
(2). Flow meter manometers were read and recorded. Wall pressure readings 
were recorded on paper tape. 

For velocity traverses, the reference position of the probe traversing 
mechanism was observed. Measurements were then made at the various desired 
stations by manipulating the traversing mechanism to the desired position and 
recording the hot-wire voltage. 
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SECTION 5 


DISCUSSION OF RESULTS 

The results will be presented in the following arrangement. The 
minimum suction requirements for unseparated flow will be shown first. Then 
measured and predicted velocity distributions at the exit planes, the walls, 
and the passage centerlines will be Shown. Further the effects of downstream 
blockage, of non-uniform inlet velocities, and of unsteady inlet velocity 
will be discussed. Finally a table of di ffuser performance characteristics 
will be presented. 

5.1 Minimum Suction Requirement 

Figure 10 shows the minimum suction requirements for unseparated flow 
through the four diffusers when approximately uniform inlet velocities exist. 

Both stable and metastable requirements are shown when operated with the design 
slot width (l/8-inch (.32 cm) for 18-inch (46 cm) diffuser, and 1/16-inch (.16 ci») 
for the other diffusers) and, for two of the diffusers, when operated with a 
slot width several times wider than the design width. A wider slot enables 
a- given suction rate to be removed with a smaller pressure drop through the 
slot. In each case the design suction requirement is shown for an inlet 
velocity of 200 ft/sec (61 m/sec). The design suction was determined by 
assuming the one-seventh power law for the boundary layer velocity profile 
and applying Taylor's criterion [7]. This criterion gives an estimate of the 
amount of fluid wich must be removed to prevent flow separation due to an 
unfavorable pressure gradient across the suction slot. The boundary layer 
fluid from the solid surface up . to where the velocity ratio 



23 



SUCTION RATE , % SUCTION RATE 




0.18 - INCH (46 cm) ANNULAR DIFFUSER AR- 3: 1 


b. 9-INCH (23 cm) ANNULAR DIFFUSER AR« 3/1 



DIFFUSER INLET VELOCITY 


■ 

1 


A 1/16" S 
A " s 

LOT, ST A 
■ , ME1 

BLE 

'ASTABLE 

■ 

■ 




■ 

■ 

■ 

■ 

■ 


■ 




19 

m 

■ 

1 





Si 

■ 


rr- , / 

m 

■ 

■ 



a L 



1 




DESIGN SI 

ICTION RA 

E 


0 50 100 150 200 250 FT/SEC 

' 1 i 1 . i . 1 

0 20 4 0 6 0 8 0 M/SEC 


DIFFUSER INLET VELOCITY 


C.9-INCH (23 cm) ANNULAR DIFFUSER 

Figure 10. Minimum suction 
inlet velocity. 


AR-2.5U d.BELL SHAPED DIFFUSER AR« 3 •: I 

requirement for unseparated flow when using uniform 


2k 







exists must be removed according to Taylor's criterion. Uj and are the 
velocities just upstream and downstream of the slot at the outer edge of the 
boundary layer, and u is the velocity within the boundary layer upstream of the 
slot along that streamline which ends in a stagnation point just downstream of 
the slot. 

The results indicate that the design suction rate was very close to that 
required for metastable operation for the 9~inch (23 cm) annular diffuser 
(AR = 2.5:1)* For the 9~inch (23 cm) annular (AR = 3:1) and for the bell 
shaped diffusers the design suction rate was slightly below that required for 
metastable operation. For the l8-inch (46 cm) diffuser, the design rate of 
5% was k0% lower than the 1% actually required for metastable operation. 

The suction rates required tended to drop asymptotically as the inlet 
velocities were increased. Since the boundary layer thickness tends to decrease 
with increased inlet velocities, this trend was expected. 

The bell shaped diffuser required slightly smaller suction rates than the 
annular diffusers having the same area ratio because only one wall requires 
suction. Dropping the area ratio from 3:1 to 2.5:1 also reduced the suction 
requirement as expected. 

The suction requirement for stable operation was of course above the 
metastable operation requirement but the differences tended to become less at 
higher inlet velocities and with wider slots. The wider slots also tended to 
make the suction requirement curves flatter or less dependent upon inlet 
velocities. 

For inlet air velocities in. the 200-250 ft/sec (61-76 m/sec) range, metastable 
operation required 6.5% suction for the bell shaped diffuser while stable 
operation required 12% suction. Corresponding suction rates for the 9- inch (23 cm) 
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annular diffuser (3:1 AR) were 8 . 5 % for metastable and 16.0? for stable 
operation. For the 9-inch (23 cm) annular diffuser with 2.5*' 1 area ratio, 

6.3$ and 12% suction rates were required respectively for metastable and stable 
operation. For the most recently designed diffuser, the 18-inch (46 cm) 
annular (AR = 3:1), 7% suction was required for metastable operation and 10% 
for stable operation at an inlet velocity of 230 ft/sec (70 m/sec). For lower 
inlet velocities the suction rate increased as shown in figure 10a. 

5.2 Velocity Distributions 

Figure 11 shows the exit plane velocity profiles for the four diffusers 
when an approximately uniform inlet velocity of 200 ft/sec (61 m/sec) is used 
and unseparated flow exists. Three sets of velocity measurements resulting from 
three circumferential locations 120° apart are shown for each diffuser. The 
velocity distribution predicted by the digital computer analysis program is also 
shown. 

The dominant feature shown by this set of results is the virtually flat exit 
velocity profile of each of the four diffusers. For the bell shaped diffuser, 
the velocity ratio (local exit velocity/maximum exit velocity) was O .98 - 0.02 
over 95% of the exit area. For the 9“ inch (23 cm) annular (AR = 3:1) the 
exit velocity ratio was 0.96 - 0.04 over 95% of the exit area. For the 
9-inch (23 cm) annular (AR = 2.5:1) the exit velocity ratio was 0.93 - 0.07 
over 95% of the exit area. For the 18-inch (46 cm) annular diffuser the 
exit velocity ratio was 0.91 i 0.09 over 95% of the exit area. The computer 
analysis program predicted a virtually flat exi t. velocity profile for 
each diffuser as shown. 

Figure 12 shows the "wall velocity" distribution for unseparated flow 
through the four diffusers when an approximately uniform inlet velocity of 
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200 ft/sec (61 m/sec) exists. The "wall velocities" were not measured directly 
but were computed from pressure coefficients obtained from wall pressure taps. 
Thus they are velocities near the wall just outside the boundary layer. For 
the three annular diffusers velocity distributions are shown for both the inner 
and the outer wall. The wall velocity distribution predicted by the computer 
analysis program is also shown for each diffuser. There is generally good agree- 
ment between the wall velocity distribution predicted and that determined from 
wall pressure measurements. 

The "wall velocity" d i s t r i but ions show the design philosophy of the Griffith 
type diffuser. There is virtually no deceleration along the diffuser walls and 
thus the "wall velocity" must change drast ical ly at the slot. 

One may note the difference between the "wall velocity" distribution 
and the "centerline velocity" distribution shown by Figure 13- The centerline 

profiles were determined from velocity measurements along a constant radius near 
the middle of the flow passage. The deceleration at the center of the flow 
passage is seen to occur over a significant part of the diffuser length in con- 
trast with the wall deceleration. For the centerline profile the axial position 
was taken from the diffuser inlet (smallest cross-sectional area), while for the 
wall velocity profile, the axial length included the nozzle or accelerating 
portion of the inlet section. 

If the velocity distribution were precisely uniform at both the inlet and 
the exit planes, the velocity ratio along the passage centerline would drop 
from 1.0 at inlet to [(1.0 - FS/AR] at exit. For example, for the 9“inch 
(23 cm) annular diffuser (AR = 3-0 the velocity ratio would, with uniform 
inlet and exit velocity distribution, drop from 1.0 to (1.0 - 0 . 2 1 ) / 3 - 0 or 
0.263 since 21% suction was applied. The measured velocity ratio dropped 
from 1.0 to 0.2A. A suction rate higher than the minimum required for stable 


30 



VELOCITY RATIO, LOCAL/MAXIMUM VELOCITY RATIO . LOCAL/ MAXIMUM 



1 . 1 8 - INCH (46 cm) ANNULAR DIFFUSER AR=3*I 



AXIAL POSITION /AXIAL LENGTH 


c. 9-INCH (23cm) ANNULAR DIFFUSER AR = 25 : I 


Figure 13. Axial velocity profile along a 
the flow passage. 


O 0.2 0.4 0.6 0.8 1.0 


b. 9 -INCH (23 cm) ANNULAR DIFFUSER AR*=3:I 



31 







operation was used in order to minimize the fluctuation of the hot wire 
measurement caused by wakes shed from the stem of the probe. 

\ 

5.3 Effect of Partial Blockage Downstream 

Figure 14 shows the minimum suction requirement for unseparated flow 
through the 1 8- i nch (46 cm) annular diffuser when a perforated sheet metal 
plate was attached to the diffuser exit plane. The plate had 5/l6~inch (0.825 
cm) diameter holes resulting in approximately 51% free area and simulated a 
partial blockage or obstruction downstream of the diffuser. There was a slight 

increase in the required suction rate at low inlet velocities when compared 

S 

with no downstream obstruction. At inlet velocities above 1 50 ft/sec (46 m/sec) 
no significant effect was observed. The velocity distribution just upstream of 
the perforated plate was determined to be very simitar to the exit plane velocity 
distribution when no downstream obstruction existed. 

Three struts located 120° apart were mounted in the downstream section of 
the 9-inch (23 cm). (AR = 2.5:1) annular diffuser to simulate flow interference 
due to structural members. Figure 15 shows the exit plane circumferential 
velocity profile when this diffuser was operated with an approximately uniform 
inlet velocity of 200 ft/sec (61 m/sec). Approximately 10% suction was used to 
prevent flow separation. The velocity ratio drops to about 0.92 as the strut 
location is approached. This amount of variation was no more than the radial 
velocity variation when no struts were present. By taking advantage of the 
nearly constant velocity zone for the location of the struts, the problem of 
flow separation was avoided. 

5.4 Effect of Non-Uniform Inlet Velocities 

A series of tests were made to determine the effect of non-uniform diffuser 
inlet velocities on the 9"inch (23 cm) annular diffuser (AR = 2.5:1). Figure 16 
shows the non-uniform duct discharge velocity profiles when a plate in the 
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shape of a half circle was mounted in the supply duct. These duct discharge 
velocity profiles were determined before the diffuser was attached to the 
end of the duct. Subsequent diffuser exit velocity measurements were made 
after installing the diffuser at the end of the partially blocked air duct. 
Figure 17 shows the resulting exit radial velocity distribution at two 
circumferential locations 90° apart. The velocity ratio variation was from 
a' maximum near the inner wall to about 82 % of maximum near the outer wall. 

This variation was only slightly more than when a uniform inlet velocity 
was used. 

Figure 18 shows the non-un i form duct discharge velocity profile when a 
14-inch (35 cm) circular plate was installed in the supply duct so as to 
yield a radially non-uniform duct discharge. Figure 19 shows the diffuser 
exit velocity profiles which were subsequently determined. As with the 
ci rcumferent i a 1 distortion, the diffuser exit radial velocity distribution 
was determined at two locations 90° apart. The velocity ratio variation was 
from a maximum near the inner wall to about 83% of maximum near the outer wall. 
The suction rates required for unseparated flow were approximately the same 
as with previous tests using uniform inlet velocities. 

5-5 Effect of Unsteady Inlet Velocity 

As previously mentioned, the effect of an unsteady inlet velocity was 
examined by installing a motor driven propeller which caused wakes to be shed 

upstream of the 9~inch (23 cm) annular diffuser inlet. 

Figure 20 shows the minimum suction requirement for stable operation without 
and with the propeller. The shedding of wakes from the propeller caused an 
increase in the suction requirement, but the difference tends to become less as 
the inlet air velocity is increased. At 200 ft/sec ( 6 1 m/sec) the suction 

requirement for stable operation was increased from 16% to 18% by the propeller. 
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5.6 Diffuser Performance 


Tests were run to determine the effectiveness of the various Griffith 
diffusers. The effectiveness of the Gr i f f i th-type diffuser should be high 
when operated without flow separation because it is a very short diffuser and 
thus friction and turbulence losses should be very low. The price which must 
be paid is the suction loss, and for the effectiveness to have real meaning, 
this loss should be accounted for. 

For this reason, this report will show both the conventional diffuser 
effectiveness, n (= actual pressure rise/ideal pressure rise) and a more 
meaningful effectiveness called e. e = (flow rate out) (actual pressure rise)/ 
(flow rate in) (ideal pressure rise). Appendix C derives the following 
expressions used in computing n and e. 


n (convent i ona 1 ) = 


(P - P .) 
s ,e s , i 


(24) 


(^rj 


e (modified) = 


(1 .0 - FS) (P - P 
s ,e 


s , \ 


(25) 


f 

(1.0 - FS ) *1 

l“i - V 

( ar ) J 


where FS is the fraction removed by suction and k is the kinetic energy coefficient 
to account for the non-uniformity of the flow at the diffuser inlet and exit 
planes. Typical values of k. and were 1.025 and 1.011 respectively. 

The logic of modifying the conventional effectiveness by the ratio of 
the exit flow rate to the inlet flow rate becomes evident when one considers 
that the conventional definition represents a ratio of energy terms per unit 
mass flowing. When no suction or injection occurs, flow rates in and out 
cancel. 


The total pressure loss (see Appendix C) was computed from 

1.0 - FS 2 - 


P t , - P t - 0-n) p . . fie 

t,i t,e d,i [ 


. “ K 

i e 


f I .0 - FS \ 1 

( ar ) j 


( 26 ) 


and non-dimensional i zed to the inlet dynamic pressure P 


d,i 


Table I shows the results of the performance tests which were made on the 
four diffusers. In the 150-250 ft/sec (46-76 m/sec) range of inlet velocities, 
the diffuser effectiveness values based on the conventional definition were from 
97% to 33 % . When slot suction was properly charged to the diffusers, the effective- 
ness values were about 86% for the 9“ inch (23 cm) annular (AR = 3*1)> 88% for 
the 9-inch (23 cm) annular (AR = 2.5:1), 89% for the 18-inch (46 cm) annular, 
and 91% for the bell shaped diffuser. A typical total pressure loss was 2% 
of the inlet dynamic pressure. Fluctuations from this value were probably 
due more to slight pressure measurement errors than to actual changes in the 
total pressure loss. 
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TABLE I. GRIFFITH DIFFUSER PERFORMANCE 

\ 


Test 

No. 

Inlet 

Ft ./Sec 

Je loci ty . 
M/Sec 

Inlet 
Reynolds 
Numbe r 

Slot 
Suet i on 
% 

P ressu re 
Recove ry 
Coefficient 
C 

P 

Effectiveness Total Pressure 

J /ojX Loss , % of 

(g. 25) (Eq 24) ' n,et ° ynamiC 

(Eq. 26) 


A. 

Bell Shaped Diffuser; AR=3 • 1 ; Slot 

Width = 

1/16 in. (.159 cm) 

1 

82 

25 

m 

7.3 

0.92 

91.5 

98.8 l.l 

2 




6.2 

0.92 

93.6 

99-0 0.9 

3 

235 

72 

12.0 x 10 4 

6.3 

0.91 

92.0 

98.2 1.7 

B. 

9 in. 

(23 cm) Annular 

Diffuser; AR=3:1; 

Slot Width = 1/16 in. (.159 cm) 

4 

79 

24 

2.0 x lO* 

16.6 

0.90 

81.3 

95-2 4.5 

5 

142 

43 

3.6 x 10** 

12.3 

0.91 

85.5 

97.5 2.3 

6 

236 

72 

6.0 x 10^ 

1 1 .2 

0.91 

86.5 

97.5 2.3 

c. 

9 i n 

(23 cm) Annular Diffuser; AR=2.5: 

1; Slot 

Width = 1/16 in. (.159 cm) 

7 

82 

25 

2.8 x 10^ 

12.8 

0.89 

86. 1 

98.8 1.1 

8 

147 

45 

5-0 X 10^ 

1 1 .0 

0.87 

86.6 

97.4 2.2 

9 

237 

72 

8.0 x to* 1 

9-5 

0.88 

89.5 

99.0 0.9 

D. 

18 in. (46 cm) Annular Diffuser; AR=3: 

; Slot 

Width = 1/2 in. (1.27 cm) 

10 

84 

26 

4.3 x 1 0 1 * 

10.3 

0.90 

86.7 

96.7 3.1 

1 1 

144 

44 

4 

7.3 x 10 4 

8.3 

0.92 

90.9 

99-0 0.9 

12 

229 

70 

11.7 x lO* 1 

10.8 

0.93 

88.5 

99-0 0.9 
















SECTION 6 
CONCLUSIONS 


The following conclusions were reached based upon the results of the 
.research program. 

(1) The concept of the Griffith diffuser is a reasonable approach to designing 
short axially symmetrical diffusers for special applications. 

(2) The potential flow design program utilizing the two dimensional wedge 
flow patch to approximate the branch flow from the suction slot yields 

a workable curved wall diffuser. This is not to rule out the possibility 
of improvement by using a more sophisticated design program such as an 
improved branch flow patch and the accounting for boundary layer 
d i spl a cement th i ckness . 

(3) Reasonably good agreement was achieved among the velocity distributions 
prescribed to the design program, predicted by the analysis program, and 
determined experimentally. 

(4) The slot suction rate required to prevent sepa rat i on w i th stable flow 
was found to decrease asymptotically with diffuser inlet velocity and 
radial dimensions. Thus the minimum suction rate was about 10% at 230 fps 
(69 m/sec) for the largest diffuser tested. 

(5) High performance can be expected from the Griffith diffuser when 

sufficient suction is used to prevent flow separation. An effectiveness 
88% is a typ ical value if slot suction is accounted for. 

(6) A reasonably uniform exit velocity profile results from the Griffith 
type diffuser when sufficient suction to prevent flow separation is 
used . 
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APPENDIX A - PROGRAM DESCRIPTION AND LISTING 


The main program is the controller for the program. It calls each of 
the subroutines as needed, and it contains the main loop for the iteration of 
the radial coordinates. 

SETUP is a subroutine which computes coefficients which are needed to per- 
form Simpson's Rule integration. 

S I MP is a subroutine for integration using Simpson's Rule. 

PART I is a subroutine which reads in, checks, and outputs the input data. 

Obvious errors in input data terminates the program after error explanation 
is printed. Certain constants to control subsequent calculations are com- 
puted and initialization of arrays Q and R are also performed in PART I. 

CQEF is a subroutine which computes constants used in the solution for the 
velocity distribution. 

RELAX is a subroutine which uses the Gauss-Siedel Technique with overrelaxation 
to solve for In Q over the region. Where derivatives of the velocity are 
prescribed as zero, this subroutine corrects those boundary conditions accordingly. 
PART II is a subroutine which performs the inverse transformation to get the 
flow angles, a, the R coordinates, and finally the X coordinates if the iteration 
on R is on its last pass. This subroutine also corrects the boundary conditions 
on R where required. The integration is carried out first along the middle 
streamline and then along each potential line individually. 

OUTPUT is a subroutine which prints out the R array upon request during the 
iteration on R. At the completion of all computations, it outputs the Q, R, X, 
and a arrays and a summary of the wall coordinates of the finished channel 



GRAPH I C is a subroutine which plots out the streamlines and equ i potent i a 1 


lines for the final channel design. Inlet and outlet channel velocity pro- 
files are optionally output. 


INPUT DATA PREPARATION INSTRUCTIONS 


The following is a description of the input data, their meanings, and how 
they must appear on the data cards. 


Control Data 
CARD 1 


Column Variable Description 


1 FLAG 1 


2 FLAG2 


3 FLAG3 


A FLAGA 


FLAG1 is the input flag for the plane flow case. 

FLAG1 = 0 for axially symmetric flow. 

FLAG1 = 1 for plane flow. 

FLAG2 is the boundary condition flag. 

FLAG2 = 0: Dirichlet B.C. are used on the inlet and outlet 
potential lines and the velocities are set equal to the 
first and last wall velocities input respectively. 

FLAG2 = 1: Neumann B.C. are used on the inlet and outlet 
potential 1 i nes . 

FLAG2 = 2: Dirichlet B.C. are used on the inlet potential 
line with the velocities being input there and Neumann 
B.C. are used on the outlet potential line. 

FLAG2 = 3 - Neumann B.C. are used on the inlet potential 
line and Dirichlet B.C. are used on the outlet potential 
line with the velocities being input there. 

FLAG2 = A: Dirichlet B.C. are used on the inlet and outlet 
potential lines with the velocities being input on both. 

FLAG3 is the bell diffuser flag. 

FLAG3 = 0 if the diffuser is plane or annular. 

FLAG3 = 1 if the diffuser is a bell diffuser. 

FLAGA is the graphical output flag. 

FLAGA = 0 i f no graphical output is desired. 

FLAGA = 1 if graphical output is desired. 


A2 



Col umn 


Variable 


Description 


7-66 

HEDR 

69-70 

DATE1 

71-72 

DATE2 

73-74 

DATE 3 

77-80 

CASE 

CARD 2 

Col umn 

Var iabl 

1-10 

RMULT 

11-20 

W 

21-30 

TOLSYS 

31-40 

TOLDE 

41-50 

MSTAGU 

51-60 

MSTAGL 

CARD 3 

Col umn 

Var iabl 


1-10 M 


HEDR is a 60 character comment that will be printed 
out at the top of each page of output. 

DATE1 is the two digit number giving the month of the 
date . 

DATE2 is the two digit number giving the day of the 
date. 

DATE3 is the two digit number giving the year of the 
date. 

CASE is a four digit case number that will be printed 
out at the top of each page of output. 


Descript? on 

RMULT is a constant, by which all computed X and Y 
coord i nates wi 1 1 be multiplied at the completion of the 
program. If RMULT = 0.0 no multiplication will take 
place. 

W is the relaxation factor. Although the optimum value 
varies with the grid size, a value of 1.9 has been found 
to work well in planar cases, 1.9 in bell cases, and 1,7 
in annular cases. 

TOLSYS is the accuracy to which the radial coordinates 
are to be computed. This accuracy is determined as the 
largest change of any radial coordinate from the pre- 
vious i terat ion. 

TOLDE is the accuracy to which the velocity is to be 
computed. This accuracy is determined as the largest 
change of any velocity from the previous iteration. 

This is the number of the nodal point in the (j> direction 
at which the upper stagnation point (slot downstream 
lip) i s to be located . \ 

This is the number of the nodal point in the <f> direction 
at which the lower stagnation point is to be located. 

For plane and bell designs this item is not used. 


Description 

The number of mesh points to be used in the <f> direction 
(The maximum number 148.) 


Col umn 


Vari able 


Description 


11-20 

21-30 

31-40 

41-50 

51-60 

CARD 4 
Col umn 
1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


N 

MX I SYS 


MX IDE 


NSTAGU 

NSTAGL 


The number of mesh points to be used in the ip direction. 
(The maximum number is 48.) 

The maximum number of iterations allowed to find the 
radial coordinates. If FLAG1 = 1 then MXISYS = t. 

MXISYS = 10 gives good results for annular and bell 
d i f fusers . 

The maximum number of iterations allowed to find the 
velocities; 300-400 usually gives good results for the 
2-D case and 50 gives good results for the axially 
symmetric case. 

This is the number of nodal points in the \p direction 
at which the upper stagnation point (slot downstream 
lip) is located. N-NSTAGU cannot be less than 3- 

This is the number of nodal points in the \p direction 
at which the lower stagnation point is located. For 
plane and bell designs this item is not used. NSTAGL 
cannot be less than 4. 


Vari able Description 


RAD IN 


MIDJ 


IPRSYS 


RADIN is the radial coordinate at the nodal point 
numbered (1,1) for annular design. If left blank its 
value is assumed to be 1.0. This item is not used in 
plane or bell designs. 

The number of the streamline for an annular diffuser. 

Usually chosen as N/2. Care should be exercised in the 
choise of MIDJ since an improper choice can cause an 
annular design to cross the X axis. This results in a 
divide check, overflow or underflow and a consequent 
termination of the program. This item is only used for 
annular designs. 

The frequency at which intermediate values of the coordinates 
of the diffuser are printed out during iteration on R. 

If no intermediate values are desired, make IPRSYS greater 
than MX I DE. 


IPRDE 


The frequency at which intermediate values of the 
velocity are printed out. 


MSL0TU 

MSL0TL 


The number of the potential line at the bottom of the 
upper slot. This number is usually only slightly larger 
than MSTAGU. 


The number of the potential line at the bottom of the 
lower slot. This item is only used for annular designs. 
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Spacing and Boundary Conditions 

The following variables are input six to a card in columns 1-10, 
11-20, etc. 


Variable 


Descri pt ion 


PHI 
Q (♦, 

upstream 

a (♦, v 

downstream 
QSLOTU (<f>) 

Q (<J>, 0) 

upstream 

Q (♦, v,) 

downstream 
QSLOTL (<p) 

PS I 

Q. (0, ip) 

Q (*, 4 >) 


The values of the velocity potential along a streamline 
at each nodal point. There should be M values. 

The values of the velocity along the upper wall upstream 
of and along the upstream wall inside of the slot. There 
should be MSLOTU of these values. 


The values of the velocity along the upper wall downstream 
of the slot. There should be M-MSTAGU of these values. 

The values of the velocity along the downstream wall 
inside of the slot. There should be MSLOTU-MSTAGU of 
these values. 

The values of the velocity along the lower wall upstream 
of and along the upstream wall inside of the slot. There 
should be MSLOTL of these values. Omit for plane and 
bel 1 des i gns . 

The values of the velocity along the lower wall downstream 
of the slot. There should be M-MSTAGL of these values. 

Omit for plant and bell designs. 

The values of the velocity along the downstream wall inside 
of the slot. There should be MSLOTL-MSTAGL of these values. 
Omit for plane and bell designs. 

The values of the stream function along a potential line 
at each nodal point. There should be N values. 

The value of the velocity at nodal points along the inlet 
potential line. There should be N values. Values are 
input only if FLAG2 = 2 or 4. 

The values of the velocity at nodal points along the outlet 
potential line. There should be NSTAGU values for plane 
and bell designs and NSTAGU-NSTAGL + 1 values for annular 
designs. Values are input only if FLAG2 = 3 or 4. 




OUTPUT EXPLANATION 


On the first page of output the control data is printed out and the 
options called for by the control data are listed. Following the printout of 
the control data comes a printout of the wall boundary conditions and spacing 
in the <p direction. The wall boundary conditions are given in three parts: 
the inlet wall data, the downstream wall data and data inside the slot. The 
first column of each set of data gives the arc length, S, computed along the 
wall streamline. The second column labeled Q gives the velocity that is input 
along the wall streamline. The third column labeled PHI gives the values of 
<p as input. The fourth column labeled DPHI gives A<j>, the spacing as computed 
by subtracting adjacent values of <j>. For annular designs three more sets of 
wall boundary conditions are output for the diffuser lower wall. 

Following the wall streamline boundary conditions the stream function 
distribution is printed out along with the spacing computed from it. The 
numbers to the left of these outputs give the number of the node to which the 
information corresponds. 

Following the input data comes a list of the convergence variables ITER, 
JTER, EPS and DEL for each iteration on R (these are printed out for the axially 
symmetric case only). ITER is the number of the iteration on R that was then 
in progress. ■ JTER is the number of iterations on the velocity that have been 
completed. EPS is the accuracy to which R has been computed. DEL is the 
accuracy of solution of the InQ. equation. 

The arrays of Q, R, X and ALPHA are output after the last iteration. 

The arrays QSLOTU, QSLOTL, RSLOTU, RSLOTL, XSLOTU, XSLOTL, ASLOTU, and ASLOTL, 
the values of the corresponding variables inside the slot on the downstream 
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wall, are also output. The values of ITER, JTER, EPS and DEL are output 
in the heading of each array. 

A summary of the wall coordinates is output last along with the 
velocity that was prescribed along the wall, the arc length computed along 
the wall and the tangent angle along the wall. This output is presented in 
three parts as was the input data along the walls. For annular channels 
another set of output is presented for a lower wall. 

The following 45 pages is a program listing for the design program. 
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********** main program ********** 

CLEMSON UNIVERSITY MECHANICAL ENGINEERING DEPARTMENT 
PROGRAM 70-02 

AXIALLY SYMMETRIC AND 2-D BRANCHED FLOW CHANNEL DESIGN 


DIMENSION DX(150> , ESIMP! 150),FS IMP! 150) ,GS IMP! 150) , 
GRAND! 150) , G I FT( 1 50) , PHI ( 150 > , PS l ( 150 ) , 

DP SI! 150) ,DPH 1(150 ),S( 150) ,SBAR( 150) , SI! 150) , 
SBAR1 (150) , S2( 150) , SBAR2 ( 150) , HEDR ( 15 ) , 

EPHI! 150) ,FPHI (150) ,GPHI( 150) ,EPSI( 150) , 

FPSI( 150) ,GPSI(150),RSL0TU(i50) ,RSLOTL< 150). 
QSLOTU! 150) .QSLOTL! 150) .ASLOTU! 150) , 
ASLOTU150) ,XSL0TU(150),XSL0TL(150) 

PRECISION P 

FLAG1, FLAG2, FLAG3.FLAG4, CASE, DATE 1.0ATE2.0ATE3 
/Cl / IMIN.SU, IMAXSU.IMINSI .IMAXSI .GRAND, GIFT, 
DX,ESIMP,FSIMP,GSIMP ^ 

COMMON / C2 / FLAG1 ,FLAG2»FLAG3»FLAG4»HEDR»0ATE1» DATE2, 
DATE 3, CASE ,MXIDE,ICKDE,IPRDE,IMIN, I MAX, 
JMIN, JMAX.NUT »M,N,TOLDE»DEL»EPS,JTER, 
ITER, DP HI, DP SI, S, SBAR ,RMULT» S1.S2.SBAR1, 
SBAR2,MSTAGU,NST AGU, MSLOTU »MST AGL ,NST AGL , 
MSLOTL,RADIN,W,PHI,PSI,RSLOTU,RSLOTL, 
QSLOTU, QSLOTL, ASLOTU, ASLOTL , XSLOTU , 
XSLOTL,G,MIDJ,TOLSYS 

EPHI ,FPHI , GPHI , EPS I ,FPSI .GPSI.MXISYS, 
IPRSYS 


DOUBLE 

INTEGER 

COMMON 


COMMON / C3 / 


CALL PARTI 
II = 1 

DO 100 I TER=1 , MXI SYS 
CALL COEF 
CALL RELAX 
CALL PART2 

IF (EPS .LE. TOLSYS) GO TO 200 
IF (ITER ,LT. I I* I PRSYS ) GO TO 100 
NUT = l 
II = II * 1 
CALL OUTPUT 
100 CONTINUE 
200 NUT = 0 

CALL OUTPUT 

IF (FLAG4 .LE. 0) STOP 

CALL GRAPIC 

end 


********** SUBROUTINE SETUP ********** 


COMPUTES COEFFICIENTS FOR USE IN SIMPSON RULE 
INTEGRATION 


SUBROUTINE SETUP 

DIMENSION DX(150),ESIMP(150),FSIMP(150),GSIMP( 150), 

1 GRAND! 150) ,GIFT(150) 

COMMON / Cl / IMINSU, IMAXSU.IMINSI, IMAXSI, GRAND, GIFT, 
1 OX,ESIMP,FSIMP,GSIMP 

DO 1000 l=IMINSU» IMAXSU 
0X1 = DX(I) 

DX2 = DX( 1+1) 
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dxiq = dxi * oxi * oxi 

0X1 SQ = 0X1 * 0X1 
0X2Q = DX2 * 0X2 * DX2 
0X2 SO = DX2 * 0X2 

XNUM = 2. * 0X10 + 3. * DX ISO * 0X2 - DX2Q 
DEN = 6 . * 0X1 * ( 0X2 + 0X1 ) 

E S I M P ( I ) = XNUM / 0EN 

XNUM = 0X20 ♦ 0X10 *■ 3 . * 0X2SQ * DXI * 3. * 0X1SQ * 0X2 
DEN = 6 . * 0X1 * 0X2 
FSIMP( I ) = XNUM / DEN 

XNUM = 2. * DX2Q ♦ 3. * DX2SQ * DXI - DX10 
OEN = 6. * 0X2 * ( 0X2 ♦ DXI ) 

1000 GSIMP(I) = XNUM / DEN 
RETURN 
END 

*#** 

SIMP 

********** SUBROUTINE SIMP ********** 

THIS SUBROUTINE OOES SIMPSON'S RULE INTEGRATION 
OF A FUNCTION. AND RETURNS THE INTEGRATED FUNCTION 

SUBROUTINE SIMP 

DIMENSION DXI 150) , ESIMP 1 1 50 ) , FS IMP 1 150 ) , GS I MP ( 150) , 

1 GRAND (150) »GIFT(150) 

COMMON / Cl / IMINSU, IMAXSU, IMINSI , IMAXSI , GRAND, GIFT, 

1 OX, ESIMP, FSIMP.GSIMP 

GIFT ( IMINSI ) = 0.0 

II = IMINSI + 2 

COMPUTE MAIN INTEGRAL 

00 1000 I«I I, IMAXSI, 2 

III = I 

SUM = GIFT! 1-2 ) + ESIMPII-2) * GRAND! I -2 ) 

SUM = SUM + FSIMPII-2) * GRAND ( I- 1 ) 

SUM = SUM ♦ GS I MP ( 1-2) * GRAND! I) 

1000 G I FT 1 1 ) = SUM 
* 

COMPUTE INTERMEDIATE VALUES 
IUP = IMAXSI 

IF (IMAXSI - III .EQ. 1) IUP * IMAXSI - 2 
II = IMINSI + 1 

00 2000 1*1 I,IUP,2 
D1 = OX! 1-1 ) 

D2 = DX! I » 

GIFT(I) = 01 * ! 2. * D1 + 3. * 02 I * GRAND! I — 1 • / 

1! 6. * 02 + 6. * D1 ) 

GIFT ( I ) = GIFT ( I ) + D1 * ( D1 ♦ 3. * 02 ) * GRANO(I) / 

1! 6. * D2 ) 

2000 GIFT(I) = GIFT! I) - D1 * D1 * D1 * GRAND! I +1 ) / 

1! 6. * D2 * 02 ♦ 6. * D2 * D1 ) «■ GIFT C 1-1 » 

IF (IMAXSI - III . EQ. 0) RETURN 

1 = IMAXSI 

GIFT! I > = ( GRAND! I -1 ) * GRAND! I ) ) * DX(I^l) / 2. ♦ GIFTCI-1) 

RETURN 

END 

********** SUBROUTINE PARTI ********** 

THIS SUBROUTINE READS INPUT DATA AND PREPARES 
ARRAYS FOR RELAXATION 
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SUBROUTINE PARTI 

DIMENSION DXll 501 »ESIMP( 150) » FS IMP! 150) »GSIMP( 150) » 
l GRAND! 150) , GIFT! 150) ,PHI ( 150) , PS 1(1 50) , 

l DPS I (150) f DPHI(150),S(150),SBAR(150),Sl(150), 

3 S8AR1 (150) , S 2 ( 150 ) , SBAR2 ( 1 50 1 ,HEDR(15) , 

i EPHI( 150 ) ,FPHI (150»,GPHIC 150) ,EPSI < 150) , 

3 F P SI ( 1501 ,GPSI(15Q),RSLQTU(15Q1, RSLOTH 1501, 

S OSLOTUI 150) tQSLOTL ( 150) ,ASLOTU( 150) i 

7 ASL0TL(150) , XSLOTU ( 150) »XSL0TL(150) 

DIMENSION Q( 150,50) ,R( 150,50) ,X< 150,50) 

DOUBLE PRECISION P 

INTEGER FLAGl, FLAG2 , FLAG3, FLAG4, CASE , DATE 1 , 0ATE2 , 0ATE3 
COMMON / Cl / IMINSU»IMAXSU,IMINSItIMAXSIt GRAND, GIFT, 
l DX,ESIMP,FSIMP,GSIMP 

COMMON / C2 / FLAGl, FL AG2 , FL AG3 , FL AG4, HEDR , DATE1 , DATE2 , 
L DATE3,CASE,MXIDE,ICK0E, IPRDE.IMIN, I MAX, 

? JMIN,JMAX,NUT,M,N, TOLDE »DEL » EPS » JTER » 

3 ITER, DPH I ,DPSI ,S, SBAR, RMUL T , S 1 , S2 , S8AR 1 , 

i SBAR2,MSTAGU,NSTAGU, MSLOTU , MSTAGL , NST AGL, 

5 MSLOTL » RAD IN ,W » PH I,PSI ,RSLOTU,R SLOTL , 

5 QSLOTU.QSLOTL.ASLOTU, ASLOTL, XSLOTU, 

7 XSLOTL,G,MIDJ,TOLSYS 

COMMON / C3 / EPHI ,FPHI,GPHI »EPSI»FPSI ,GPSI , MX I SYS, 

L IPRSYS ' 

REWIND A 
REWIND 5 
REWIND 12 


READ CONTROL DATA 


! READ( 1,10) FLAGl, FLAG2.FLAG3, FLAG4, HEDR ,DATE1 , OAT E2, 

1 DAT E3, CASE, RMULT, W.TOLSYS, TOLDE, MSTAGU, MSTAGL, 

2 M,N,MXISYS,MXIDE,NST AG U, NST AGL ,RADIN,M IDJ , 

3 IPRSYS, IPROE .MSLOTU, MSLOTL 

10 FORMAT ( 4 1 1 , 2X , 15A4 , 2X , 3 12 , 2X , I4/4FI 0 .0, 2 1 10/61 10/ 

1 F10. 0,5110) 

OUTPUT CONTROL DATA AND OPTIONS 


WRITE(3,20) 

20 FORM AT ( 1H1//T57, ’CLEMSON UNI VERS 1 T Y «// T49 , 

1 * MEC HAN I CAL ENGINEERING DEPARTMENT*) 

WRITE (3, 30)HEDR,DATE1,DATE2,DATE3,CASE ,FLAGl , RMULT, 

1 M, RAD IN, FLAG2 » W »N »MID J »FL AG3 » TOL SYS »MXI SYS, 

2 IPRSYS , FL AG4, TOLDE , MX IDE, IPRDE,MSTAGU,NSTAGU, 

3 MSLOTU, MSTAGL .NSTAGL, MSLOTL 
30 FORMAT! //T59,’ PROGRAM 70-02*//T38, 

1 • AX I ALL Y SYMMETRIC AND 2-0 BRANCHEO CHANNEL’, 

2T85 , • DESIGN’ ///Til ,’*****• ,T21, 15A4.T84, I2.T86 ,’ /• , 
3TS7, I2.T89, ’/• , T90 , 1 2 , T96 , • CASE NO. *,T107,I4, 

'4Tllfe, •**♦**’ ,/////T51, • CONTROL DATA • 

5/////T16, ’FLAGl = • ,T2 3 , 1 1 , T 36, • RMULT = *,T45,F11.6, 
6T66 , * M = * ,T75 , 1 10.T96, ’RADIN = • , T 105 , F 10. 5/ 

7T 16 , • FLAG2 = • , T2 3 , 1 1 , T3 6, • W = • , T45 , F 1 1 . 6 , T66 , 

8 ’ N = • ,T75, I 10.T96, ’MIDJ = • , T 105 , I 10/T16, 

9 • FL AG3 = ’,T23,I1,T36,’T0LSYS = * , T45 , F 1 1 . 6 , T66 , 
A’MXISYS = ’ ,T75,I 10.T96, • IPRSYS = • , T 1 05 , 1 1 0/T16 , 
B’FLAG4 = • ,T23 , II ,T36, ’ TOLDE = • , T45, FI 1 . 6 , T66 , 
C’MXIDE = • ,T75,I 10.T96, • IPRDE = • , T1 05 , I 10/T36 , 
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D * MST AGU = * ,T46, I 10,T66, *NSTAGU = *,175,110, 

ET96,*MSL0TU = 1 ,T105,110/T36, 'MSTAGL = *,T46 f llO, 

FT66 , ' NS T AGL = 1 , T 75 , I 10 , T96 , * MSLOTL = * , T 1 05 , I 1 0// / / 

GT55 , • OPTIONS •) 

EPS = 0.0 

IF (FLAGi .GT. 0) WRITE(3,40> 

40 FORMAT ( / //T 57 * •2-DIMENSIONAL* ) 

IF (FLAGI .LE. 0) WRITE(3,50) 

50 FORMAT! ///T57, 'AXIALLY SYMMETRIC*) 

IF ( FLAG2 .EQ. 0 .OR. FLAG2 .EQ. 4 ) WRITE(3,55) 

55 FORMAT(/T57, •DIRICHLET 80UNDARY CONDITIONS 1 ) 

IF ( FLAG 2 • F.0 . 1 ) WRITE(3,60) 

60 FORMAT (/T57 , ‘NEUMANN BOUNDARY CONDITIONS*) 

IF ( FLAG2 .EQ. 2 ) WRITE(3,65) 

65 FORMAT ( /T57 , •DIRICHLET B. C. AT INLET AND NEUMANN *, 

1 T94 , * B . C. AT OUTLET* ,/T60, 'WITH VELOCITY INPUT AT INLET*) 

; IF ( FLAG2 .EQ. 3 ) WRITE(3,70) 

70 FORMATI/T57, ‘NEUMANN B. C. AT INLET AND DIRICHLET *, 

1 T94 , • B • C. AT OUTLET* ,/T60,* WITH VELOCITY INPUT AT OUTLET*) 

IF ( FLAG3 .GT. 0) WRITE(3,80) 

8 C FORMAT! /T57, ‘BELL DIFFUSER - NOZZLE*) 

IF (FLAG3 .LE. 0 .AND. FLAGI .LE. 0) WRITE(3,90) 

90 FORMAT ( /T57 t * ANNUL AR DIFFUSER - NOZZLE*) 

IF (FLAG4 .GT. 0) WRITE(3,100) 

100 F0RMAT(/T57, 'GRAPHICAL OUTPUT*) 

IF ( FLAG4 .LE. 0) WRITE(3,110) 

110 FORMAT (/T57 , *N0 GRAPHICAL OUTPUT SPECIFIED*) 

OUTPUT ERRORS DETECTED 

IF (FLAGI .LE. 0) GO TO 130 
IF ( FLAG3 .LE. 0) GO TO 130 
/ WR I TE ( 3 , 120 ) 

120 FGRMATdHlt Tit *OIF FLAGI GREATER THAN ZERO FLAG3 MUST BE ZERO*) 
STOP 

130 IF (FLAGI .LE. 0) GO TO 150 
IF (MXISYS .EQ. 1) GO TO 150 
WRI T E ( 3 , 140 ) 

140 FORM AT ( 1 H 1 , T1 , • 0 I F FLAGI GREATER THAN ZERO THEN* , T32 , 

1* MXISYS MUST BE ONE*) 

STOP 

150 NM1 = N - 1 

NP1 = N + 1 

MM1 = M - 1 

MP 1 - M + 1 

NM2 - N - 2 

NP2 = N ♦ 2 

MM2 — M — 2 

MP2 = M + 2 
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NSTUP1 = NSTAGU ♦ 1 
NSTUP2 = NSTAGU ♦ 2 
DO 155 J=1»NP2 
DO 155 I=1»MP2 

0(1 t J) = 1.0 

155 R( I , J) = 1.0 
DO 157 1=1, MP2 
RSLOTU(I) = 1.0 
RSLOTL(I) = 1.0 
QSLOTU(I) = 1.0 
157 QSL0TL ( I ) = 1. 0 

READ WALL VELOCITY DISTRIBUTION 

READ (1,160) ( PH I < I ) , I =2 , MP1 ) 

READ( 1, 160) (Q( I ,NP1) ,I=2,MSLUP1 ) 

READ! 1, 160) (0(1 , NSTUP 1 ) , I=MSTUP2 , MP 1 ) 

READ( 1,160) ( QS LOTU ( I ) , I=MSTUP2 , MSLUP1 ) 

160 FORMAT(6F10.0) 

IF (FLAG1 .GT. 0 .OR. FLAG3 . GT . 0) GO TO 170 
READ( 1, 160) (Q( 1,2), 1 = 2 , MSLL PI ) 

REA 0(1,160) ( Q( I , NSTL P 1 ) , I=MSTLP2 i MP 1 ) 

R EAD ( 1,160) ( QS LOT L ( I ) , I = MSTL P2 , MSLL PI ) 

170 DO 100 1=2, M 

160 DPHMI) = PHHl + 1) - PH I ( I ) 

DPHI ( 1) = DP H I (2) 

OPHI(MPl) = DPHI (M ) 

READ STREAM FUNCTION SPACING 

RE AD (1,1 60) (PS Id ),I = 2,NP1) 

DO 190 1=2, N 

190 DPSI (I) = PSH I ♦ 1 ) - PSIU ) 

DPSI (1) = OPSI <21 
OPSI(NPl) = DPSI(N) 

COMPUTE ARC LENGTHS ALONG THE WALL 

- DO 200 1=2, M 
200 nx( I ) = DPH I ( I ) 

IMINSU = 2 
IMAXSU = MM 1 
CALL SETUP 
DO 210 I =2 , MSL UP1 
210 GRANnCU = 1.0 / Q ( I , NP 1 ) 

. IMINSI = 2 
I MAX S I = MSLUP1 
CALL SIMP 
DO 220 1=2, MSL UP1 
220 S ( I ) = GIFT ( I ) 

SKMSTUPl) = 0.0 
S2 ( M STUP 1 ) = 0.0 
Q ( MS TUP l , NS TUP 1 ) = 0.0 
QSLOTU(MSTUPl) = 0.0 

SKMSTUP2) = 2. * DPH I ( MSTUP 1 ) / Q ( MS TUP2 , N STUP 1 ) 
$2 ( M STUP2 ) = SKMSTUP2) 

DO 230 I = MS TUP 2 ,MP 1 
230 GRAND(I) = 1.0 / Q(I,NSTUP1> 

IMINSI = MSTUP 2 
IMAXSI = MP 1 
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CALL S T MP 

DO 240 I = MS TUP 3 ,MP1 
24 C SKI) - GIFT(I) + SKMSTUP2) 

DO 250 I=MSTUP2,MSLUP1 
25 G OR AND ( I ) = 1.0 / QSLOTU(I) 

1 Ml NS I = MS TUP 2 
IMAXSI = MSLUP1 
CALL SIMP 

DG 260 I=MSTUP3,MSLUP1 
260 S2« I I = GIFT(I) + S2(MSTUP2) 

IF (FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) GO TO 330 
DO 270 I =2 , MSL L PI 
270 GRAND! I ) = 1.0 / 0(1,2) 

IMINSI = 2 
IMAXSI = MSLLP1 
CALL SIMP 
DO 280 1 = 2, MSL L PI 
280 SB AR ( I ) = GIFT! I ) 

SBARi(MSTLPl) = 0.0 
SBAR 2 ( MSTLP1 ) = 0.0 
Q(MSTLP1,NSTLP1 ) = 0.0 
QSLOTL ( MSTL PI ) = 0.0 

S8AR 1 ( MSTLP2 ) = 2. * DPHI(MSTLPl) / Q ( MSTLP2 , NSTLP 1 ) 
SBAR 2 ( MSTLP2 ) = SB AR1 ( MSTL P2 ) 

DO 290 I=MSTLP2,MP1 
290 GRANO(I) = 1.0 / Q(I,NSTLP1) 

IMINSI = MSTLP2 
IMAXSI = MP1 
CALL SIMP 

DO 300 I=MSTLP3,MP1 

300 SBAR 1 ( I ) = GIF T ( I ) + SB AR1 ( MSTL P2 ) 

DO 310 I=MSTLP2,MSLLP1 
310 GRAND(I) = 1.0 / QSLOTL(I) 

IMINSI ^ MSTLP2 
IMAXSI = MSLLP1 
CALL SIMP 

DO 320 I-=MSTLP3»MSLLP1 
320 SBAR 2 ( I ) = GIFT(I) * S B AR 2 ( MSTL P2 ) 

OUTPUT WALL VELOCITY DISTRIBUTION ANO VELOCITY 
POTENTIAL DISTRIBUTION ON UPPER SURFACE 

33 C WR1 TE(3,20) 

WRITE (3,340 ) HEDR , DATE'l , DATE2 , DATE3 , CASE 
340 F0RMAT(//T2 ,15A4,T64,I2 ,T66, f /• ,T67, I2,T69, ■/• ,T70, 12, 
1T76, * C A S E NO. 1 , T87 , I 4/// T2 , • WALL VELOCITY', 

2T 1 5 , ' DISTRIBUTION AND POTENTIAL SPACING', 

3T50, • UPSTREAM OF SLOT ( UPPER SURFACE )'///T 16, ‘S', 
4T36, # Q' ,T54, 'PHI » , T74, 'DPHI •// ) 

ICOUNT = 1 
00 360 1=2, MSL UP 1 
ICOUNT = ICOUNT + 1 

WRITE (3, 350) I, S( I) ,Q( I ,NP1 ) ,PHI ( I ) , DPHI ( I) 

35 0 FORM AT (T2,I 3,T10,F10.5,T30,F10.5 ,T50,F10.5/T70,F10.5) 
IF (ICOUNT .LT. 23) GO TO 360 
ICOUNT = 1 
W R I T E ( 3 , 2 0 ) 

WR I TE ( 3, 340) HE DR, DATE 1, DAT E2,DATE3, CASE 
360 CONTINUE 

WR I TE ( 3 , 20 ) 
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WRITE(3,370)HE0R,DATE1,0ATE2,DATE3,CASE 
370 FORMAT ( / / T2 , 15A4,T64, 12 *T66 , • / • , T67 , I 2 , T69 , • / • ,T70, 12, 
1T76, 'CASE NO. • , T87 , I 4/ / / T2 , • WALL VELOCITY', 

2T 15, ' DISTRIBUTION. AND POTENTIAL SPACING', 

3T50, ' DOWNSTREAM OF SLOT ( UPPER SURFACE > • /// T16 , ' S ' , 
4T 36 * ' Q' , T 54 » 'PHI • , T74 , ' DPH I ' // ) 

I COUNT = 1 

DO 3 SO I = MS TUP 1 , MP I 

ICOUNT = ICOUNT + 1 

WRITE (3, 3 50) I , S 1 ( I I ,Q( I ,NSTUP1 ) , PH l ( I ) , DPH I ( I ) 

IF (ICOUNT .LT. 23) GO TO 380 
ICOUNT = 1 
W R I T E ( 3 , 20 ) 

WRITE(3,37O)HE0R,DATEljDATE2,DATE3,CASE 
380 CONTINUE 

WRITE(3»20) 

WRITE(3,390)HEDR,0ATE1,DATE2,DATE3,CASE 
39 0 FORM AT ( //T2, 1 5 A4 , T6 A , I 2 , T66 , • / • ,T67, 12,769, •/• ,770,12, 
1T76 , 'CASE NO. • ,T37 , 1 4///T2, ' WALL VELOCITY', 

2T 15 , ' DISTRIBUTION AND POTENTIAL SPACING', 

3T50, ' INSIDE SLOT ( UPPER SURFACE ) V//T16, • S' , 

4T36, 'Q* ,T54, 'PHI ' , T74, 'DPH I • // ) 

ICOUNT = 1 

DO 400 I=MSTUPl,MSLUPi 
ICOUNT = ICOUNT + 1 

WRITE (3, 350)1 , S2(I) ,OSLOTU( I ) , PH l ( I ) ,DPHI ( l ) 

IF (ICOUNT .LT. 23) GO TO 400 

ICOUNT = 1 

WRITE(3,20) 

WRITE(-3,390)HEDR,DATE1,DAT£2,DATE3,CASE 
40 C CONTINUE 

Q( MSTUPi , NSTtJP l ) = 1.0 
OSLOTU(MSTUPl) = 1.0 

OUTPUT WALL VELOCITY DISTRIBUTION AND VELOCITY 
POTENTIAL DISTRIBUTION ON LOWER SURFACE 

IF (FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) GO TO 470 
WR I T E ( 3 , 20 ) 

WRI TE (3,410 )HEDR, DA TE1 , DAT E2 , DAT E3 , CASE 
410 F0RMAT(//T2,15A4,T64,I2,T66,'/' , T67 , I 2 , T69 , • / • ,T70,I2, 
1 T 76 , 'CASE NO. • , T 87 , I 4/7 / T2 , • WAL L VELOCITY', 

2T 1 5 , ' DISTRIBUTION AND POTENTIAL SPACING', 

3T 50 , ' UPSTREAM OF SLOT ( LOWER SURFACE I '///T16, • S' , 
4T36, 'O' , T54, 'PHI • ,T74, • DPH I' // ) 

ICOUNT = 1 
DC 4-20 1=2, MSL L PI 
ICOUNT = ICOUNT +1 

WRITE (3 ,350 ) I , SBAR ( I ) ,U ( I ,2 ) , PHI ( I ) ,DPHI ( I ) 

IF (ICOUNT .LT. 23) GO TO 420 
ICOUNT = 1 
WR I TE ( 3 , 2 0) 

WRITE(3,410)HEDR,DATE1,DATE2,DATE3,CASE 
42 C CONTINUE 

WRITE! 3, 20) 

WRITE(3,430)HEDR,DATE1,DATE2,DATE3,CASE 
430 FORMAT (/ /T2, 15A4.T64, I2,T66,'/' ,T67,I2 ,T69, '/• ,T70, 12, 
IT 76, 'CASE NO. ' ,T37 , I4///T2, 'WALL VELOCITY', 

2T 15 , • DISTRIBUTION AND POTENTIAL SPACING', 

3T 50 , ' DOWNSTREAM OF SLOT ( LOWER SURFACE ) '///T16, ' S' , 
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4T36 * * 0’ »T54» ‘PHI* , T74 , • DPHI • // ) 

ICOUNT = 1 

DO 4 40 I=MSTLP1,MP1 

ICOUNT = ICOUNT «• 1 

WRITE (3, 350) It S8AR1 ( I) ,Q( I ,NSTLP1), PHI ( I > ,DPHI C I ) 

IF (ICOUNT .LT. 23) GO TO 440 
ICOUNT = 1 
WR t TE ( 3 » 2 0 ) 

WRITE(3,430)HEDR,DATE1,DATE2,DATE3,CASE 
44 C CONTINUE 

WRITF(3,20) 

WRI TE (3,450 )HE0R,0ATE1, OAT E2tOATE3, CASE 
450 FORMAT ( // T2 , 1 5 A4 , T64 , I 2 , T66 , 1 / ' , T67 , 1 2 , T69 , • / • ,T70, 12, 
1T76, ‘CASE NO. • , T87, I4///T2, • WALL VELOCITY*, 

2T15 , ' DISTRIBUTION AND POTENTIAL SPACING*, 

3T 50 , * INSIDE OF SLOT ( LOWER SURFACE > • ///T 16, • S ' , 

4T36 , *0* ,T54, ‘PHI* , T74 , * DPH I * // ) 

ICOUNT = 1 

DO 460 I=MSTLPI,MSLLP1 
ICOUNT » ICOUNT ♦ I 

WRITE (3, 350 ) I ,SBAR2( I ) ,QSLOTL< I ) , P H I ( I >,OPHI ( I ) 

IF (ICOUNT .LT. 23) GO TO 460 

ICOUNT =1 

WRITE(3,20) 

WRITE(3,450)HEDR,DATE1,DATE2,DATE3,CASE 
460 CONT INUE 

Q( MSTLP1 , NSTLP 1 ) = 1.0 
OSLOTL(MSTLPl) * 1.0 

OUTPUT STREAM FUNCTION DISTRIBUTION AND 
' SPACING 

470 WRI T E(3, 20) 

WRI TE( 3,430 )HEDR, DATE 1 , OAT E2 , DA TE3 , CASE 
480 FORMAT (//T2» l 5A4,Tfe4, I2,T66,* /• ,T67, l2,T69,*/‘ ,T70,I2» 

• 1 T 76 , ' CA SE NO. *, T87 , I 4///T2, • STREAM FUNCTION*, 

• 2 T IS » ’ OISTRIBUTI ON *,///T15,*PSI* , T35 , * OP SI * // ) 

ICOUNT =1 

DO 500 J=2 , NP1 

ICOUNT = ICOUNT +1 

WRI TE (3,490) J , PS I ( J ) , DP S I ( J ) 

490 FCRMAT(T2,I3,T10,F10.5,/T30,F10.5) 

IF (ICOUNT .LT. 23) GO TO 500 

ICOUNT =1 ' 

WRITE (3,20) 

WP 1 T E ( 3 , 480 ) HE DR , DATE 1 , D AT E2 , DAT E 3 , CA S E 
500 CONTINUE 

COMPUTE INTEGRATION CONSTANTS INVOLVING OPHI 

DO 510 1=2, M 
5 1 C 0X( I I = OPHI ( I ) 

IMINSU = 2 
IMAXSU = MM1 
CALL SETUP 
00 520 1=2, MM1 
E D H I (I) = ESIMP(I) 

F PH I ( I ) = FSIMP< I ) 

520 GPHI ( 1 ) = G S I M P ( I ) 
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. COMPUTE INTEGRATION CONSTANTS INVOLVING 0PSI 
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DC 530 J = 2 ♦ N 

PARI 

367 

530 

DX(J) = DPS I < J ) 

PARI 

368 


IMINSU = 2 

PARI 

369 


IMAXSU = NM1 

PARI 

370 


CALL SETUP 

PARI 

371 


DO 540 J=2,NM1 
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EPSI (J) = ES IMP ( J I 

PARI 
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FPSI ( J) = FSIMP(J) 
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540 

GPS I ( J) = GSIMP( J) 
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COMPUTE CONSTANTS TO CONTROL INTEGRATION 

PARI 
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I MI N = 3 

PARI 

379 


IMAX = M 

PARI 

380 


JMIN = 2 

PARI 

381 


JMAX = NP1 

PARI 

382 


IF { FLAG3 . GT . 0) JMIN=3 

PARI 

333 


IF ( FLAG2 .LE. 0) GO TO 550 
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384 


[MIN = 2 

PARI 

385 


IMAX = HP 1 
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INITIALIZE ARRAYS 
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389 

550 

DQ 560 J=2,NP1 

PARI 

390 


Q ( 2 * J ) = Q( 2 » N PI ) 

PARI 

391 

560 

Q(MP1,J> = 0< MPl tNSTUPl ) 

PARI 

392 


DO 570 J = 1 » NP2 

PARI 

393 


DO 570 1= 1 » MP 2 

PARI 

394 


P * G(I,J) 

PARI 

395 

570 

Q( I t J) = DLOGl P) 

PARI 

396 


DO 580 I = 1 1 MP2 

PARI 

397 


P = QSLOTUI I ) 

PARI 

398 


QSLOTU ( I > = DLOG( P) 

PARI 

399 


P = QSLGTL ( l > 

PARI 

400 

580 

QSLOTL(I) = DLOG(P) 

PARI 

401 
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COMPUTE VELOCITIES AROUND STAGNATION POINTS 

PARI 

403 
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404 


Q(MSTAGU,NSTUP1I - QS LOTU ( MSTUP2 1 

PARI 

405 


QIMSTUPI ,NSTUP2) = QSLQTU( MSTUP2 ) 

PARI 

406 


Q ( M S TUP 1 »NSTAGU) = QS LOTU < MSTUP2 ) 

PARI 

407 


IF ( FLAG 1 • GT. 0 .OR. FLAG 3 .GT. 0) GO TO 610 

PARI 

408 


Q(MSTAGLtNSTLPl) - QS LOT L < MSTL P2 > 

PARI 

409 


Q( MSTLP1 ,NSTLP2) - QSLOTL ( MSTL P2 ) 

PARI 

410 


Q(MSTLPlrNSTAGL) = QSLOTL ( MSTLP2 ) 

PARI 

411 

610 

CONT INUE 

PARI 

412 


DO 620 J = I ? NP 2 

PARI 

413 


DC 620 I = 1 1 MP2 

PARI 

414 

620 

X ( I f J ) = 0. 0 
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READ INLET AND OUTLET VELOCITY PROFILES (IF USED) 
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PARI 

418 


IF (FLAG2 .NE. 2 .AND. FLAG2 .NE. 4) GO TO 640 

PARI 

419 


I M I N = 3 

PARI 

420 


READ ( I v 160 ) (Q( 2 f J) «J-2«NP1) 
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421 


DC 630 J=2,NP1 

PARI 

422 


P = Q ( 2 t J ) 

PARI 

423 

63 C 

Q ( 2 » J ) = DLOG(P) 
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640 IF ( FLAG2 .NE. 3 .AND. FLAG3 .NE. 4) GO TO 660 PARI 

IMAX = M PARI 

NSTUP1 = NSTAGU ♦ 1 PARI 

IF (FLAG1 .GT. 0 .OR. FLAG3 .GT. 0 ) READ ( 1 , 160M Q ( MP1 , J > , J=2 , PARI 

1NSTUP1) . PARI 

NSTLP1 ~ NSTAGL + 1 PARI 

IF (FLAG1 .LE. 0 .AND. FLAG3 .LE. 0 ) RE AD ( 1 , 160 ) ( Q ( MP 1 f J ) , PARI 

1J=N$TLP1,NSTUP 1 ) PARI 

DO 650 J = 2 , NP1 PARI 

P = Q( MP 1 f J ) PARI 

650 'O(MPltJ) = OLOG(P) PARI 

660 CONTINUE PARI 

PARI 

WRITE HEADING FOR CONVERGENCE VARIABLES PARI 

PARI 

IF (FLAG1 .LE. 0) WR IT E ( 3 , 670 ) HEDR , DATE 1 , DATE2 , DATE3 , CA SE PARI 

670 FORMAT ( 1HI//T57, 1 CLEMSON UNIVERSITY* , / / T49 , * MECH AN I C AL E NGI NE E RI P AR1 
1NG DEPARTMENT * ,///T2, I 5A4 , T64, 1 2 , T66 f • / • ,T67, I 2 ,T69, •/• ,T70, 12 t PARI 
2 T 76 f • CA SE NO. • ,T87,I4///) PARI 

IF (FLAGI .LE. 01 WR I TE ( 3,680 ) PARI 

680 F0RMAT(T2,« IMPORTANT VARIABLES FOR EACH ITERATION OF • , PARI 

1 T43 , * THE R COORDINATES*//) * PARI 

WRITEJ4) ( ( Q ( I , J ) ,I=1,MP2),J=1,NP2) PARI 

WRITE (5) < (R ( I f J) , I=1,MP2) , J=1,NP2) PARI 

WRITE (12) (<X(I,J),I=1,MP2),J=1,NP2> PARI 

RETURN PARI 

END PARI 

********** SUBROUTINE COEF ********** COEF 

THI-S SUBROUTINE COMPUTES COEFFICIENTS TO BE USED COEF 

IN RELAXATION OF THE DIFFERENTIAL EQUATION COEF 

‘ ' COEF 

SUBROUTINE COEF COEF 

DIMENSION DPHI (150) , DPS I ( 150 ) , S ( 150 ) ,S8AR<150), COEF 

1 SI ( 150) , S2< 150 ) ,SBAR1 < 150) ,SBAR2< 150) , COEF 

2 PHI ( 1 50) , PSI ( 150) ,RSLOTU( 150) r RSLOTL( 150) f COEF 

3 QSLOTU( 150) ,QSLOTL( 150) ,ASLOTU( 150) , COEF 

4 ASLOT L ( 1 50 ) , X S LOTU ( 150 ), XSLOTL < 1 50 ), HEDR < 15) COEF 

DIMENSION A(150, 50) ,B( 150, 50) , C( 150 , 50) t 0( 150, 50) , COEF 

1 E( 150,50) ,R(150t50) COEF 

INTEGER FLAG1, FLAG2 ,FLAG3, FLAG4, CASE, DATE 1 , DAT E2, OAT E3 COEF 


COMMON / C2 / FLAG 1 , FLA G2 , FLAG3 , FLAG4 , HEDR , DATE 1 , DATE2 , COEF 



DAT E3 , CASE, MX IDE, ICKDE , I PRDE , I MI N , I MAX , 

COEF 

) 

JMIN, JMAX , NUT, M f N,TOLDE ,DEL, EDS , JTER, 

COEF 


ITER , DPHI ,DPSI ,S, SBAR , RMUL T , S 1 , S2 , SBAR1 , 

COEF 

► 

S BAR2, MSTA GU, NSTAGU, MSLOTU,MSTAGL, NSTAGL, 

COEF 

i 

MSLOTL,RADIN ,W,PHI ,PSI , RSLOTU , RSLOTL , 

COEF 

> 

OSLOTU,OSLOTL,ASLGTU, ASLOTL ,XSLOTU, 

COEF 


XSLOTL, GtMIDJfTOLSYS 

COEF 

REWIND 

5 

COEF 

REWIND 

6 

COEF 

PEW INO 

7 

COEF 

REWIND 

8 

COEF 

REWIND 

9 

COEF 

REWIND 

10 

COEF 

MP2 = M 

+ 2 

COEF 

NP2 = N 

♦ 2 

COEF 

NSTUM1 

= NSTAGU - 1 

COEF 

MSTUP2 

= MSTAGU + 2 

COEF 

MSTLP2 

= MS T AG L +2 

COEF 

MSLUP1 

- MS LOT U ♦ 1 

COEF 
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VSLLP1 - MS LOT L ♦ 1 
NSTLP1 = NSTAGL ♦ 1 
NSTLP2 * NSTAGL + 2 
NSTUP2 ~ NS TAG U + 2 
READ (51 ( <R< I,J I , Kl,MP2) , J=1,NP2) 

COMPUTE COEFFICIENTS THROUGHOUT DIFFUSER 

DO 100 J = J MI N t JMA X 
CD = DPSKJ-ll * DPSI(J) 

EE = DPSKJ-l) * DPSI(J) 

FF OPSI ( J-l) - OPSK J) 

DO 100 KIMIN, IMAX 

AA = DPMI ( I - 1 1 * DPHII I I 

BB = DPHKI-l) + DPHI(I) 

CC = DPHKI-l) -DPHICI) 

OERX = OPHKI-ll * / ( DPHKII * 88 I 

DEP.X = DERX - DPHKII * R(I-1,JI / ( DPHKI-1) * BB ) 

DERX = DERX - CC * R(I,J| / AA 

OERY = DPSKJ-ll * R(I«J + 1I / ( DPSI(J> * EE J 

DERY = DERY - DPSKJI * R(I»J-1I / < DPSKJ-i) * EE > 

DERY = DERY - FF * R(l f J) / DD 

DOER X = 2 . * R ( I ♦ 1 * J ) / ( DPHim * BB I 

ODER X = DDE RX + 2* * RU-1»JJ / < DPHI(I-l) * BB ) 

DDERX * DDE RX - 2. * R(I,J) / AA 

F = 2. * R ( I f J I * R ( I f J I * R ( I » J ) * R(ItJ) / DD 

F = F + R ( I » J ) * R ( I v J ) * R ( 1 1 J I * FF * DER Y / DD 

F * F ♦ 2 • * R(I,J) * R(lvJ) / AA 

F = F - R(I , J) * CC * DERX / AA 

A € I f J ) = 2. * R ( I « J ) - 0PHIU-1I * DERX 

A ( I ? J ) = A ( I t J ) * W * R ( I t J ) / ( F * DPMI ( I ) * BB ) 

B ( I * J ) = 2. * R(I,J) > DPHI(l) * DERX 

6(1, J) = B( 1 9 J ) * W * R ( 1 9 J ) / ( F * DPHKl-1) * BB ) 

R3 - R(I f J) * R(IiJI * Rd9JI 

C(ItJ) - 2. * R ( I i J ) + DPSI(J-l) * DERY 

C(ItJ) = C < I * J ) * W * R3 / ( F * DPSI(J) * EE) 

DllfJ) = 2. * R(I,J) - DPSI(J) * DERY 

D ( I 9 J ) = D ( I v J ) * w * R3/ ( F * DPSKJ-ll ♦ EE ) 

100 E ( I 9 J ) = ( R ( I t J ) * DDERX - 2. * DERX * DERX ) / F 

COMPUTE COEFFICIENTS INSIDE UPPER SLOT 

J = NSTAGU + 2 

DO = DPSKJ-ll * OPSI(J) 

EE = DPSI { J-l) ♦ DPS I ( J ) 

FF ^ DPS I( J-l ) - OPSI ( J ) 

DO 200 !=MSTUP2tM$LUPl 
A A = OPH I ( [-11 * DPH I ( I ) 

66 = DPH I ( 1-1) ♦ DPMI ( I-| 

CC = DPH I ( 1-1) - DPHI ( I ) 

DERX = DPHI(I-l) 4c R ( I + 1 f J ) / ( DPHKII 4c 3B ) 

DERX = DERX - DPHKII * R(I-1,J) / ( OPHKI-l) * BB ) 
DERX = DERX - CC * R(l f J) / AA 

DERY = DPSI(J-l) 4c R(I t j + l) / ( DPSKJI * EE ) 

DERY - OERY - DPSKJI * RSLOTU(I) / ( DPSI(J-l) * EE ) 

OERY = DERY - FF * R(I,J) / DD 

DDERX = 2. * R(i + l,J) / ( DPHKII * BB ) 

DDERX = DDERX + 2. * R(I-l f J) / ( DPHKI-ll * BB I . 
DOER X = DDERX - 2 . * R(I,J) / AA ' 

F = 2. * R(I>J) * RCItJI 4c R ( I t J ) 4c R ( 1 1 J ) / DD 
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F = F + R(I F J> * R ( I 9 J K * RCItJ) * F F * OERY / DD COEF 

F = F + 2* * R(I,J) * R(UJ\ / AA COEF 

F - F - R ( I f J ) * CC * OERX / AA COEF 

Ail i J) = 2. * RilrJ) - DPHII1-1) * DERX COEF 

Ail » J ) = A(ItJ) * W * R(ItJ> / ( F * OPHI (i) * B 8 > COEF 

B ( I v J ) = 2* ♦ R ( I F J ) ♦ DPHKI) * DERX COEF 

- B ( I » J l * W * R ( I » J ) / ( F * OPHUt-ll * BB ) COEF 

R 3 = R ( I 9 J ) * R < I t J ) * R ( I t J ) COEF 

C(I,JI = 2 . * RCItJ) ♦ DPS K J- 1 ) * DERY COEF 

CUtJI * CUtJ) * W * R3 /.( F * DPSI(J) * EE > COEF 

D(I>J) = 2* * R ( I 9 J ) -DPSKJ) * OERY COEF 

O(IiJ) = D< I » J 1 *W*R3/(F* OPSKJ-l) * EE ) COEF 

200 EU,J) = ( RlhJ) * DDE RX - 2 . * DERX * DERX ) / F COEF 

COEF 

COMPUTE COEFFICIENTS INSIDE LOWER SLOT COEF 

COEF 

IE (FLAG1 .GT. 0 .OR. FLAG3 *GT * 0) GG TO 1000 COEF 

J = NSTAGL COEF 

00 = DPSI(J-l) * DPSHJ) COEF 

EE = DPSKJ-1) + DPSI(J) COEF 

FF = DPS I ( J -1 ) - DPSI(J) COEF 

DO 300 I=MSTLP2,MSLLP1 COEF 

AA = DPHI(I-l) * DPHKI) COEF 

BB = DPH 1(1 — 1) + DPHKI) COEF 

CC = DPHKI-1) - DP H I ( I I COEF 

DERX = DPHIII-1) * RU + ltJ) / ( DPHKI) * BB ) COEF 

' OEP X = DERX - DPHKI) * R(I-lfJ) / ( DPHKI-1) * BB ) COEF 

DERX - DERX - CC * R < I , J I / AA COEF 

Oc»Y = OPSI(J-l) * RSLOTLII) / ( OPS I < J > * EE ) , COEF 

OERY = DERY - OPSKJI * R(ItJ-l) / ( DPSKJ-1) * EE ) COEF 

OERY = DERY - FF * R ( I » J ) / DD COEF 

OOERX = 2. * RK + ltJI / ( DPHI(I) * BB ) COEF 

DDERX = DDE RX + 2. * R(i-l f J) / ( DPHKI-1) * BB ) COEF 

ODER X = DDERX - 2. * RU»J) / AA COEF 

F = 2 • * R ( I 9 J ) * K ( I 9 J ) * RdtJ) * R ( 1 9 J ) / DO COEF 

F = F ♦ R ( I t J ) * R ( I t J ) * R < I i J ) * FF * DERY / DD COEF 

F = F + 2. * R(I»J) * R(I»J) / AA COEF 

F = F - P(I,J) * CC * DERX / AA COEF 

AU»J) = 2. * R ( I 9 J ) - DPHKI-1) * DERX COEF 

A ( I 9 J ) = A ( 1 9 J ) * W * R(I f J) / ( F * DPHI(I) * BB ) COEF 

B ( I » J ) = 2* * R ( I ♦ J ) ♦ DPHKI) * DERX COEF 

B(ItJ) * B(I»J) * W * R ( I 9 J ) / ( F * OPHHi-l) * BB ) COEF 

R3 = R < I 9 J ) * R ( I 9 J ) * R(I.J) COEF 

C 1 1 9 J ) = 2. * R ( I 9 J ) ♦ DPS It J- 1 ) * DERY COEF 

C(IfJ) - C ( 1 9 J ) * W * R3 / ( F * DPSKJ) * EE ) COEF 

LXItJ) « 2. * R ( 1 9 J ) - DPSKJ) * DERY COEF 

0(1 9 J) - D< I 9 J ) * W * R3 / ( F * DPS It J-l ) * EE > COEF 

300 EII 9 J) = ( RdtJ) * DDERX - 2. * DERX * DERX ) / F COEF 

1000 CONTINUE - COEF 

G = 1 • - W COEF 

WRITE (6 K (A(I f J)tI=l f MP2> f J*ltNP2> COEF 

WRI TE(7K (BU 9 J) 9 1 = 1 »MP2) ,J*1 9 NP2) COEF 

W R I T E ( 8 ) ( (C( I 9 J ) 9 I = 1 9 HP 2 ) 9J = l9NP2) COEF 

W R I T E ( 9 ) ( (D( I,J) »I = ltMP2) tJ"l»NP2) COEF 

WRITE (10) ( ( E< I f J) 9 1 = 1, MP2I , J=1,NP2) COEF 

RETURN COEF 

END * COEF 
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SUBROUTINE RELAX 


** 4c ****** * 


THIS SUBROUTINE RELAXES THE DIFFERENTIAL EQUATION 
USING COEFFICIENTS COMPUTED IN COEF 

SUBROUTINE RELAX 

DIMENSION DPHI f 150) t DPSI ( 150) t$f 150) , SBAR ( 150) , 

1 SI < 150) f S2C 150), SB ARK 150) , SBAR2C 150) * 

2 PHI (150) 9 PS I (150) i RSLO T U ( 150),RSL0TLC 150) , 

3 QSLOTUC 150) ,QSLOTL( 150) , ASLOTU ( 150 I , 

A A S LOT L (130) ♦ X S LOTU C 1 50 ) , XSLOTL ( 1 50 ) , HE DR ( 15) 

DIMENSION A < 150,5 0) ,B< 150,50) ,C< 150 ,50) ♦ DC 150,50) , 

1 E ( 150,50) ,Q(150,50) 

DOUBLE PRECISION P 

INTEGER FLAG1, FLAG2 , FLAG 3, FL AG4 , C AS E , DA TE 1 , DA TE 2 , DA TE3 
COMMON / C 2 / FLAG1 ,FLAG2,FLAG3,FLAG4,HEDR ,DATE1 , 

1 DATE2, DAT E3, CASE, MX IDE, ICKOE , IPRDE, I MIN, 

IMAX tJMINtJMAX. NUT tMfNf TOLD E, DEL * EPS t 
JTER , ITER ,OPHI ,DPSI , S, SBAR, RMULT, SI ,S2 , 
S8AR 1 * S3AR2, MSTAGUfNSTAGU, MSLOTUtMSTAGL , 
NST AGL , MS LOT L , RAD l N , W , PH I , P S I , RS LOTU , 
RSLOTL ,QSLOTU,QSLOTL , ASLOTU, ASLOTL , 
XSLOTU , XSLOTL ,G,MIDJ,TOLSYS 


READ (A) ( (Q( I , J ) , 1= 1, MP2) , J=1,NP2) 
READ (6) ( ( AC I ,J ) ,I = 1,MP2 ) ,J = 1,NP2) 
READ (7) ( (B( I , J ) ,1 = 1, MP2 ) , J=1,NP2) 
READ (?) ( (C( l , J ) , 1 = 1 ,MP2) , J=1,NP2 > 
READ ( 9 ) ( (0( I « j ) ,f-l »MP2) f J- It NP2) 
READ (10) ( (El l,J) ,1=1, MP2) ,J=1,NP2) 
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RELX 

35 

NM2 = N - 2 
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11 = 1 
J TER = 0 

00 3000 I TT T = 1 #MX 1 0E 

DEL = 0.0 

JTER = JTER + 1 

RELAX LOWER PORTION FOR 2-0 AND BELL 

IF (FLAG1 . L E • 0 .AND. FLAG3 .LE. 0) GO TO 30 
DO 10 J = 3 * N S T U M 1 
DO 10 I = I M I N » I MAX . 

AA = A( I , J) * G ( I + 1 1 J ) 

B8 = (3 ( I , J) * Q ( I-1,J) 

CC = C( I , J) * Q ( I t J +1 ) 

DD = D< I , J) * G( I , J-l ) 

EE = E < I r J ) + G * Q ( I » J ) 

QLAST = Q ( I » J ) 

0(1 t J) = AA + BB + CC + DD f EE 
DELSTR = ABS (QLAST - Q(I«j!) 

10 IF (DELSTR .GT. DEL) DEL = DELSTR 
IF (FLAG3 .GT. 0) GO TO 110 
J = 2 

DO 20 I = I MI N » I MAX 

•A A = A( I , J) * Q( I + 1 f J ) 

BB = 6 ( I, J) * Q ( I - 1 t J ) 

CC = C( I * J ) * Q( I t J + l I 
DO - 0( I , J) * 0(1 ♦ J-l) 

EE = E(!»J) ♦ G * O(IiJ) 

QLAST = Q ( I f J ) 

Q(ItJ) = AA * B8 + CC + DC ♦ EE 
DELSTR = A8 S ( QL AS T - Q(I,J)> 

20 IF (DELSTR .GT. DEL) OEL = DELSTR 
GO TO 110 

RELAX LOWER PORTION OF ANNULAR DIFFUSER 

3 C DO AO J = 3»N$TLM1 
DO 40 I-IMI N f MSLLP1 
A A = A( I f J) * Q ( I + 1 , J ) 

BB = 8( I , J) * 0( I-l, J ) 

CC = C( I, J) * 0 ( I »J+1) 

DD = 0(1 . J) * Q ( I , J-l) 

EE = £(IfJ) + G * Q(IfJ) 

QLAST = Q( I , J) 

Q ( I , J) = AA + BB + CC + DD + EE 
DELSTR = AB S ( QL AST - Q(I f J)) 

40 IF (DELSTR .GT. DEL) DEL = DELSTR 
J = NSTAGL 
DO 50 I — I M I NtMSTAGL 
AA = A( I , J) * Q< I+1,J ) 

BB = B( I ♦ J) * Q ( I-l f J ) 

CC = C( I , J) * Q< 1 v J + l ) 

DO = D( I t J) * 0( I , J-l ) 

EE = E(ItJ) + G * Of l t J) 

QLAST = Q ( I f J ) 

Q ( I t J ) = AA f BB + CC + DD + EE 
DELSTR = ABS(OLAST - Q ( I , J ) I 
5 C IF (DELSTR .GT . DEL) DEL = DELSTR 
DO 60 I = MST LP2 ? MSLL Pi 
AA = A( l • J ) * Of l + l « J ) t 
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BB = B ( I , Jl * 9(1-1, J) 

CC = C( I , J) * OSLUTLl 1 ) 

DC = D( I , J» * Oil , J-l ) 

EE = E( I, J) ♦ G * 0< I, J) 

QLAST = Q( I , J 1 

Q ( I , J I = AA + BB + CC + OD + EE 
OELSTR = A B S ( '1 L A S T - 0(I,J>> 

6 C IF ( DEL STR .GT, DEL ) DEL = DELSTR 
J = NSTAGL + 1 
, DO 70 I=IMIN,MSTLM1 
AA = A( I , J) * Q( [+1, J ) 

BB - R( I , J) * 0 C I - 1 , J ) 

CC = C.< I, J> * 0(1 , J+l) 

OD = D< I , J) * Q< I , J-l ) 

EE = E ( I , J ) ♦ G * 0(1, J) 

QLAST = 0(1 * J > 

Q( I , Jl = AA ♦ BB + CC + CD + EE 
DELSTR = ABS(OLAST - y<I,J)> 

7 C IF (DELSTR .GT, DEL) DEL = DELSTR 
J = NSTAGL <r 2 
DO 8 0 1= I MI NfMSTAGL 
A A = A ( I , J ) * Q ( I + 1 , J » 

RB = B ( I , J » * Q( I- 1 , J 1 

CC = C( I , J) * Q( I , J + l ) 

DO = D( I , J » * 0(1, J-i) 

EE = E ( I , J ) + G * Q ( I , J ) 

QLAST = Q( I , J) 

0< I , J) = AA + 89 CC + DD t EE 
DELSTR = AB S ( QLAST - 0(1, J)) 

80 IF (DELSTR .GT. DEL) DEL = DELSTR 
DC 90 I=RSTLP2,IMAX 


AA = 

A( I , J) 

* 

Q (I + 1 f J I 

P8 = 

B< It J) 

* 


CC = 

C( I , J) 

* 

QC I , J + l ) 

00 = 

D( I , Jl 

* 


EE = 

E ( I f J ) 

+ 

G * Q ( I f J ) 

QLAST = 0(1 

t J) 


Q ( I , 

J) = AA 

+ 

3R + CC + DD 


DELSTR = A B S ( Q L A S T - Q(I,J>) 

90 IF (DELSTR .GT . DEL) DEL - DELSTR 
00 100 I = I M I N , IMAX 
DC 100 J=N$TLP3,NSTUM1 
AA s A(I,J) * 0(1+1, J) 

BB = B( I , J) * Q( 1-1, J ) 

CC = C(I*J) * 0(1, J+l) 

00 = D( I , J) * 0(1 , J-l) 

EE - E ( I , J ) «• G * 0(1, J) 

QLAST = 0(1 ,J) 

0(1, J) = AA + 3B + CC ♦ DO + EE 
DELSTR = AB S ( OL AST - 0(1, J)) 

100 IF (DELSTR .GT. OEL) DEL = DELSTR 

RELAX UPPER PORTION FOR ALL KINDS OF DIFFUSERS 

110 J = NSTAGU 

DO 120 1= IMIN.MSTAGU 
AA = A ( I , J > * 0(1 + 1, J) 

BB = B( I , J) * 0(1-1, J) 

CC = C( I, J) * 0(1 ,J + 1 ) 

DO •= D(I,J) * 0(1, J-l) . 
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EE = E( I, J) + G * Q( I, J) 

QLAST = Q< I , J ) 

Q(I»J) = A A + 86 + CC + DO + EE 
DELSTR = AB S (QLAST - Q(i,J)> 

120 IF (DELSTR .GT. DEL ) DEL = DELSTR 
DO 130 I =MSTUP2 » I MAX 
AA = A( ! f J) * Q(I+L,J) 

BB = B( I , J ) * Q( 1-1, J) 

CC = C< I,J) * Q( If J+l ) 

00 = 0(l,J) * Q ( I ♦ J-l) 

EE = E(I»J) + G * 0(1, J) 

QLAST = Q( I , J) 

Q< I , J ) = AA + BB + CC + DD + EE 
DELSTR = AB S ( QLAST - 0(1, J)) 

1 3 C i F ( DEL STR .GT. DEL) DEL = DELSTR 
J = NSTAGU + 1 
' DO 140 I = 1 MI N, MSTUM 1 
AA = A ( I f J I * Q ( I + 1 , J ) 

BB = • B( I , J) * 0 ( I , J ) 

CC - C(IfJ) * 0(1, J+l) 

DD = D( I , J) * 0 ( I , J-l ) 

EE - E( 1« J) + G * 0( I , J ) 

QLAST = 0 ( I ♦ J ) 

Q( I , J) = AA + 33 .+ CC + DD + EE 
DELSTR = ABS( QLAST - Q(I,J)) 

14C IF (DELSTR .GT. DEL) DEL = DELSTR 
J = NSTAGU + 2 
DO 150 I = I M I N, MST AGU 
A A = A( I , J) * Q( 1 + 1, J) 

BB = B(.I, J) * G(I-ltJ) 

CC = C( I , J) * 0 ( I , J+l ) 

DD = D( I« J) * Q( I, J-l) 

EE - E ( I , J ) + G ♦ Q(I,J) 

QLAST = Q( I , J) 

Q ( I , J ) = AA + BB + CC + DD + EE 
DELSTR = ABS( QLAST - Q(I,JI) 

150 IF (DELSTR .GT . DEL) DEL = DELSTR 


DO 

160 I=MSTUP2,MSLUP1 

AA 

= A ( I , J > 

* 

Q( 1+1, J ) 

BB 

= B( I , J) 


Q< 1-1 ,J) 

CC 

= C( I, J) 

* 

Q( I t J +1 ) 

DD 

= 0(1, J) 

* 

OSLOTUI I ) 

EE 

= EC I, J) 


G * Q ( I , J ) 

QLAST = Q ( I , 

J) 


CJ ( I 

t J ) = AA 

♦ 

B B ♦ CC + 


DELSTR = A B S ( 0 L A S T - Q(I,JM 
160 IF (DELSTR .GT . DEL) DEL = DELSTR 
'DO 170 J=NSTUR3,N 
DO 1 7 C I = IMIN f MSLUP1 
AA = A( I f J) * Q (I + 1 , J ) 

B8 = 8 < I » J ) * 0< 1-1, J) 

CC - C( I,J) * 0(1 , J*1 ) 
no = D ( I , J ) ♦ Q( I , J-l ) 

EE = E C I , J ) + G * Q ( I , J ) 

QLAST = 0 ( I , J ) 

0(1 ,JI = AA + 88 + CC + DO + EE 
DELSTR = A B S ( 0 L A S T - Q(I,J)) 

170 IF (DELSTR .GT. DEL) DEL = DELSTR 

CORRECT BOUNDARY CONDITIONS 
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234 
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236 

237 

238 

239 
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00 180 1=2, MSL UP 1 
Y3 = 0(1, NP1 ) 

Y2 = 0(1 » N ) 

Y 1 = Q ( I , NM 1 ) 

01 = 0PSKNM1I 

02 = OPSI (N ) 

DER = ( Y3 - Y2 ) * < D1 * D2 ) / ( 01 * 02 ) 

OER = OER - ( Y3 - Y1 ) * 02 / ( 01 * 01 ♦ 01 * 02 ) 
180 Q ( 1 » NP2 ) - Y2 ♦ 1 . * OER * 02 

IF CFLAG1 .IE. 0 .AND. FLAG3 .LE. 0) GO TO 200 
JOT = 2 

DO 190 1=2, M P 1 
190 9(1,1) = 0(1,3) 

200 I F ( F L AG1 .GT. 0 .OR. FLAGS .GT. 0) GO TO 220 
JOT = NS TL PI 
DO 210 1 = 2 , MSLL PI 

Y 3 = 0(1,4) 

Y 2 = 0(1,3) 

Y 1 = 0(1,2) 

01 = OPSI (2) 

D2 = OPS I ( 3 ) 

OER = ( Y2 - Y1 ) * ( D1 ♦ D2 ) / < D1 * C2 ) 

DER = DER - ( Y3 - Y1 I * D1 / ( D1 * DZ «• D2 * D2 ) 

210 9(1 , 1) = Y2 - 2. * DER * D1 
220 00 230 J = 2 • NP 1 
Y3 = 0(4, J) 

Y 2 = Q( 3 , J) 

Y 1 = 0(2, J) 

01 = DPH I (2 ) 

D2 = DPH I ( 3 ) ' 

OER - I Y2 - Y1 ) * ( 01 ♦ D2 ) / « D1 * 02 1 

DER = DER - ( Y3 - Y 1 ) * 01 / ( D1 * D2 + D2 * D2 ) 

23 C Q ( 1 , J ) = Y2 - 2. * DER * 01 
DO 240 J=JQT,NSTAGU 
Y 3 = Q ( MP 1 , J ) 

Y 2 = 0 ( M i J ) 

Y 1 = Q(MMltJ) 

01 = OPHI(MMl) 

D2 = OPHI(M) 

DER = ( Y3 - Y2 ) * ( D1 «■ D2 I / ( D1 * D2 > 

DER = DER - ( Y3 - Y1 I ♦ 02 / I 01 * 01 + 01 * 02 ) 

240 U ( MP 2 , J ) = Y2 + 2. * DER * 02 

IF ( FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) GO TO 260 
00 250 J = 2 , NST AGL 
25 C Q(MSLLP2,J) = Q(MSLOTL,J) 

QSLOTL( MSLLP2) = 0 S LOTL ( MSLOTL ) 

26 C DO 270 J = NSTUP2,NP1 

27 C Q ( MSLUP2 , J ) = Q(M$LOTU,J) 

QSL0TU(MSLUP2) = US LOTU ( MSLOTU ) 


RE LX 241 
RELX 242 
REL X 243 
RELX 244 
RELX 245 
RELX 246 
RELX 247 
RELX 248 
RELX 249 
RELX 250 
RELX 251 
RELX 252 
RELX 253 
RELX 254 
RELX 255 
RELX 256 
RELX 257 
RELX 258 
RELX 259 
RFLX 260 
RELX 261 
RELX 262 
RELX 263 
RELX 264 
RELX 265 
RELX 266 
RELX 267 
RELX 268 
RELX 269 
RELX 270 
RELX 271 
RELX 272 
RELX 273 
RELX 274 
RELX 275 
RELX 276 
RELX 277 
RELX 278 
RELX 279 
RELX 280 
RELX 231 
RELX 282 
RELX 283 
RELX 234 
RELX 285 
RELX 236 
RELX 287 
RELX 288 
RELX 299 
RELX 290 



IF 

( FLAG3 .LE. 0) GO TO 205 





RELX 

29 1 


DO 

280 1 = 2, Mt>l 





RELX 

292 

280 

Q ( I 

,2) = Q( I ,3 » 





RELX 

293 

285 

IF 

(FLAG2 .LE. 01 GO TO 1000 





RELX 

294 


IF 

( FLAG1 .LE. 0 .AND. FLAG3 

.LE 

. 0) 

GO 

TO 310 

RELX 

295 


IF 

( FLAG2 « EQ. 2 .OR. ELAG2 . 

EQ. 

4) 

GO 

TO 295 

RELX 

296 


DO 

290 J=2, NP1 





RELX 

297 

29C 

0(1 

,J) = Q ( 3 , J ) 





RELX 

299 

295 

IF 

(FLAG2 .EQ. 3 .OR. FLAG2 . 

EQ. 

4) 

GO 

TO 1000 

RELX 

299 


DO 

300 J = 2 » NST AGU 





RELX 

300 
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300 Q(MP2tJ> = Q ( M , J ) 
GO TO 1000 


310 

IF 

( FLAG 2 . EQ. 

2 

.OR. 

FLAG2 .EQ. 

4) 

GO 

TO 

325 


DO 

320 J = 2 » NP1 








320 

Q ( 1 

, J) = 0 ( 3 t J ) 








325 

IF 

(FLAG2 .EQ. 

3 

• OR. 

FLAG2 .EQ. 

4) 

GO 

TO 

1000 


DO 

330 J=NSTLP2 

, 

NS T AGU 







33 C G(MP2,J) = Q(M,J) 

1000 CONTINUE 

CHECK FOR CONVERGENCE 

IF (DEL . LE • T0LDE) GO TO 4000 
IF ( JTER .LT. I I * I PR OE ) GO TO 3000 
11=11+1 

OUTPUT Q ARRAY 

DO 1010 I = 2 , M P 1 
DO 1010 J=2,NP1 
P = 0( I t J ) 

101C Q< I , J) = DE XP ( P) 

WRITEC 3 , 102 0 )HE0R, DATE 1, DATE 2, DAT E3, CASE 
1020 FORMAT ( 1H1//T57, • CLEMSON UNIVERSITY* , //T49, 

1 1 MEC HAN I CAL ENGINEERING DEPARTMENT 1 / //T2 , 1 5A4, T64, 
212 ,T66 , 1 /• t T67 f 12, T69 , • /• ,T70, 12, T76, •CASE NO. • , 
3T67, 14) 

OUTPUT ARRAY 

WRITE (3 f 1030 H TER t JTER, EPS f DEL 
1030 FORMAT( ///Tl, • OVELOCITY THROUGHOUT DIFFUSER'- NOZZLE 1 , 
lT46t * ITER = ST53, I4,T6l f • JTER = • , T6 8 , I 4 , T 76 , • EP S = • 
2tT82tE13*6tT99» •DEL = 1 , T105, E 1 3, 6, ///T7 , • PHI • , T66 , 

3 * PS I * ) 

WRITE (3, 1040) ( PS I ( IP ) , IP=2tNPl ) 

1040 FCRMAT(//(13X,5(7X»E13.6)n 
DO 1050 1=2, MP1 

1050 W R I T E I 3 * 1 06 0 ) P H I ( I ) t ( Q ( I,J) , J = 2 , N P 1 ) 

1060 FORMAT ( / / IX , E 1 3 . 6, 5 ( 7X » E 13 . 6 ) , / ( 14X » 5 ( 7X , E 13 . 6 ) ) ) 

CHANGE 0 BACK TO LOGARITHMIC FORM 

DO 1070 I = 2 t M P 1 
DO 1070 J=2,NP1 
P = Q( I , J) 

1070 Q< I , J) = DL CG ( P ) 

DO 1075 I =2 , MP 1 
P = QSLOT L ( I ) 

QSLOTK I ) = DEXP(P) 

P = QSLOTUt I ) 

1075 OSLO TU ( I ) = DE XP { P ) 

C 

C OUTPUT 0 VALUES ALONG THE SLOT INNER WALLS 

C 

WRIT E{ 3, 102 0IHEDR, DATE 1, 0ATE2, DAT E3, CASE 
WRITE (3, 1080 M TER, JTER tEPSt DEL 
1 08 C F0RMAT(///T1, •OVELOCITY INSIDE OF SLOTS*, 

IT 46 ITER = 1 ,T53»I4,T61, f JTER = • , T66 , I 4 , T 76 , • EPS = • 
2,T82,E13.6,T99,*DEL = # , T 1 05 , E 1 3 . 6 , / / / T7 , • PH I * , T66 , 
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RELX 351 
RELX 352 
RELX 353 
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RELX 355 
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RELX 357 
RELX 358 
RELX 359 
RELX 360 



ooooon non 


3' PS I • ) 

WRI TE! 3, 1040) ( PHI ( I P) , IP=2 ,MP1 ) 

IF { FLAG I .GT. 0 .OR. FLAG 3 .GT. 0) GO TO 1100 
WRITE! 3,1 09 0 I 

1090 FQRMAT!///T2,' LOWER SLOT') 

WRITE(3»106 0)PSI ( NSTL P 1 ) , (QSLOTL! I ) , 1 = 2, MP1 ) 
HOC WRITc(3,1110) 

1110 FQRMAT(///T2, ' UPPER SLOT') 

WRITE! 3,1 06 0)PSI(NSTUP1 ) ,( QSLOTU ( I I » I=2»MP1 ) 


CHANGE SLOT VALUES BACK TO LOGARITHMIC FORM 


DO 1120 1=2, MP1 
P = QSLOTU I ) 

QSLOTU I) = l)LOG(P) 

P = QSLOTU! I ) 

1 120 QSLOTU! I ) = DLOG! P > 

3000 CONTINUE 
4000 REWIND 4 

WRITE (4) ( !Q< I, J) ,I=1,MP2) ,J=1,NP2) 

RETURN 
END 

**:).**$** ** SUBROUTINE PART2 . ***$xi***<i* 

THIS SUBROUTINE COMPUTES- ALPHA , R AND X ARRAYS 
AND COMPUTES THE CHANGE OF THE R ARRAY FROM THE 
PREVIOUS ITERATION FOR AXISYMMETRIC COMPUTATIONS 


SUBROUTINE PART 2 

DIMENSION DX<150) ,ES IMP! 150), F SIMP! 150), GS IMP (150) , 

1 GRAND! 150) .GIFT! 150) , PHI I 150) ,PSI( 150) , 

2 DPSI l 150 ) ,DPHI ( 150) ,S( 150) , SBAR ( 150) , 

3 SI 1 150) , SRARH 150 ) ,S2< 150 ) ,SBAR2 ( 150) , 

4 HEOR! 15 ) ,.EPHI ( 150 ) ,FPHI! 150) ,GPHI ! 150) ,. 

,5 , . EPS I « 150) , ASLOTL! 150) , XSLOTUl 150) , XSLOTL! 150) 

6 , FPS1 ( 150) ,GPSI I 150) .RSLUTU! 150) .RSLOTU 150) , 

7 QSLOTU! 150) .QSLOTL! 150) ,ASL0TU(150) 

DIMENSION Q< 150,50) ,R( 150, 50) , X! 150,50) .ALPHA! 150,50 ) , 

1 R 1 < 1 5 0 , 5 0 I 

DOUBLE PRECISION P 

INTEGER F L AG1 , FLAG2 , FL AG 3, FL AG4 , C ASE , D ATE l , DATE 2 * DATE3 
COMMON / Cl / IMINSU, IMAXSU.IMINSI , IMAXSI .GRAND, GIFT, 

1 DX,ESIMP,FSIMP,GSIMP 


COMMON / C2 

1 . 

2 

3 • 

4 

5 
.6 

7 

COMMON t C 3 

1 


/ FLAG! ,FLAG2,FLAG3,FLAG4,HEDR,DATE1,DATE2 
.DATE 3, CASE, MX IDE , ICKDE, IPRDE, I MIN, I MAX, 
JMIN, JM4X ,NUT,M,N, TOLDE.OEL ,EPS, JTE'R, 

I TER, DP HI » OPS I » S » SBAR , RMULT , S 1 » S2 » SBAR 1 , 
SBAR2 ,MSTAGU,NSTAGU,MSLOTU,MSTAGL ,NSTAGL 
, MSLOTL , RADI, N,W, PHI, PSI, RSLCTU, RSLOTL ,' 
QSLOTU, QSLOTL, ASLOTU, ASLOTL, XSLOTU, 
XSLOTL ,G,MIDJ,TGLSYS 
/ EPHI ,FPHI , GPHI.EPS I ,FPSI , GPS I , MX I SYS , 
IPRSYS 


REWIND 4 
REWIND 5 
REWIND 11 
REWIND 12 
NM1 = N - 1 
N P 1 = N + 1 

MM1 = M - 1 
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RELX 

372 

RELX 

373 

RELX 

374 

RELX 

375 

RELX 

376 

RELX 
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RELX 

370 

RELX 

380 

RELX 

381 

RELX 

382 

PAR2 

1 

PAR2 

2 

PAR2 

3 

PAR2 

4 

PAR2 

5 

PAR2 

6 

PAR2 

7 

PAR2 

8 

• P A R 2 

9 

PAR2 

10 

PAR 2 

11 

PAR 2 

12 

PAR 2 

13 

PAR2 

14 

PAR2 

15 

PAR2 

16 

PAR? 

17 

PAR2 

18 

PAR2 

19 

PAR2 

20 

PAR2 

21 

PAR2 

22 

PAR 2 

23 

PAR2 

24 

PAR 2 

25 

PAR2 

26 

P AR2 

27 

PAR2 

28 

PAR2 

29 

PAR2 

30 

PAR 2 

31 

PAR 2 

32 

PAR2 

33 

PAR2 

34 

PAR 2 

35 

PAR2 

36 

PAR2 

37 

PAR2 

38 
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*P1 = M + 1 


PAR2 

39 

NP2 = N + 2 


PAR2 

40 

MP2 = M + 2 


PAR2 

41 

NM2 = N - 2 


PAR2 

42 

MM2 = M - 2 


PAR2 

43 

MSTUM1 = MSTAGU - 1 


PAR2 

44 

MSTUP1 = MSTAGU + 1 


PAR2 

45 

MSTUP2 = MSTAGU * 2 


PAR2 

46 

MSTUP3 = MSTAGU + 3 


PAR2 

47 

MSTLM1 = MSTAGL - 1 


PAR2 

48 

MSTLP1 = MSTAGL ♦ 1 


PAR2 

49 

MSTLP2 = MSTAGL ♦ 2 . 


PAR2 

50 

MSTLP3 = MSTAGL + 3 


PAR2 

51 

MSLUP1 = MS LOT U + 1 


PAR2 

52 

MSLLP1 = MS LOT L + 1 


PAR2 

53 

NSTUP1 = NS T AG U + 1 


PAR2 

54 

NSTUP2 = NSTAGU + 2 


PA P.2 

55 

NSTLP1 = NSTAGL +1 • 


P AR2 

56 

NSTLP2 = NSTAGL * 2 


PAR2 

57 

«EA0(4> ( ( 0( I ,J ) ,1 = 1 ,MP2 ) t J = 1,NP2) 


PAR2 

58 

. REAO(5M(R(I.J) »I = 1»MP2) ,J = l,NP2) 


PAR2 

59 

READ (12 ) ( ( X < I« J) t I = 1*MP2)»J=1»NP2) • • 


PAR2 

60 



PAR2 

61 

COMPUTE STARTING ALPHAS ALONG MIDDLE STREAMLINE 

P A R 2 

62 



PAR2 

63 

DO 5 I=1.MP1 


PAR 2 

64 

ASLOTIM U = 1.0 


PAR2 

65 

5 • ASLOTLC II = 1.0 


PAR2 

66 

.00 10 I =2 t MMl 


PAR2 

67 

E S I M P ( I ) = EPHI U ) 


PAR2 

68 

PSIMP( I) = FPHKII 


PAR2 

69 

1C . G S I M P < I ) = GPHI (I) 


PAR2 

70 

K =2 


PAR2 

71 

, • , If (FLAGl .LE. 0 .AND. ELAG3 .LE. 0) K =.MIDJ «- 1 

PAR2 

72 

/ J - K 


PAR2 

73 

D1 = OPSl ( J-l) 


, PAR2 

74 

02 = DPSI(J) 


PAR 2 

75 

'DO 20 I =2 » M P 1 


PAR2 

76 

OER = 01 * Q<I,J+1) / ( 02 * 01 + 02 . * 02 ) 


PAR2 

77 

D6R = OER - 02 * Q(ItJ-l) / ( 01 * 01 ♦ 01 * 

D2 ) 

PAR 2 

78 

OEK = OER - ( 01 - 02 I * Q ( l , J ) / ( 01 * D2 

) 

PAR2 

79 

2 C GRAND ( I ) = P < I , J) * OER 


PAR2 

60 

I MINS I = 2 


PAR2 

81 

IMAXSI = MP1 


PAR2 

82 

00 3 0 I = 2 t M P 1 


PAR 2 

83 

30 0X< I) = DPMI ( I ) 


PAR2 

84 

CALL SIMP 


PAR2 

35 

DO AO I = 2 , M P 1 


PAR2 

86 

4 C ALPHA! I t J | = GIFTU J 


PAR2 

37 



PAR2 

68 

COMPUTE ALPHA INTEGRANDS AND STORE THEM 

IN R 1 

PAR2 

39 



PAR2 

90 

JOT = 2 


PAR? 

91 

IF ( FLAG 3 • GT. 0) JOT = 3 


PAR2 

92 

DO 70 1-2, MP1 


PAR2 

93 

• 01 = OP H I ( I -11 


. PAR? 

94 

02 = DPMI ( I ) 


PAR 2 

9 5 

DO 7 0 J- J CT , NP 1 


P A R 2 

96 

OER = 01 * Q ( I + 1 1 J ) / < 02 * 01 ♦ 02 * 02 ) 


PAR? 

97 

JE W = OER - 02 * 0(1-1 » J) / ( 01 * 01 + 01 * 

02 ) 

PAP2 

98 


67 



non 


DER - DER — ( 01 — 02 ) * Q(IiJ) / ( Cl * 02 1 
Rl< I ,J) = - DER / R ( I t J ) 

DER = D1 * R ( I +1 , J ) / C 02 * 01 > D2 * D2 I 

DER = DER - 02 * R(I-1,J) / < D 1 * 01 * D 1 * D2 ) 

OER = DER - ( 01 - D2 ) * R(l,J) / < 01 * 02 ) 

7 C R 1 ( I , J ) = - DER / ( R(I,J) * R(ItJ) ) ♦ R1(I,J) 

IF (FLAG3 .LE. 0) GO TO 60 
DO 50 1=2, MP1 
50 R1C I 1 2 ) = 0*0 
60 CONTINUE 

COMPUTE ALPHAS IN THE UPPER PORTION OF THE DIFFUSER 

I = mstagu 

J = NSTUP1 
01 = DPHKI-1) 

D2 = DPH I ( I ) 

DER = ( 0(1 ,J) - Q< 1-1, J) ) / 01 
R1C I ,J> = - DER / R(I,JI 

DER * 01 * R ( I +1 » J I / i D2 * 01 ♦ D2 * 02 ) 

DER = DER - D2 * R(I-1,J) / ( D1 * 01 ♦ Dl* 02 ) 

DER = DER - I 01 - D2 ) * / ( 01 ♦ 02 I 

Rl(ItJ) = - DER / < R ( I » J ) * R ( I t J I ) ♦ R 1 ( I , J ) 

I = MSTUP 2 
J = NSTUP1 
01 = DPHKI-1) 

D2 = OPHI ( I ) 

DER = ( Q< I +1 , J ) - Q(ItJ) ) / D2 
Rl(IfJ) - - OER / R(l»J) 

DER =01* R(I+1,J) /( '02 * 01 * 02 * 02 I 

DER = DER - D2 * P(l-l,J) / ( 01*01*01*021 

DER = DER - ( 01 - 02 I * R(I,J) / ( D1 * D2 ) 

Rl(IfJ) = - DER / ( R ( I » J ) * R ( I » J ) ) + Rl(l,J) 

IF ( FLAG1 .GT. 0 .OR. FLAG3 .GT. 0 I GO TO 80 
I = MSTAGL 
J = N STL P 1 
D 1 = DPH I ( 1-1) 

• 02 = OPHI (I ) 

DER = ( Q ( I ♦ J ) - Q( 1-1, J) ) / 01 
R I ( I » J ) = - DER / R ( I t J ) 

DER = 01 * R ( I + 1 , J ) / ( 02 * Dl ♦ 02 * D2 ) 

DER = DER - D2 * R(I-1,J) / < D1 * D1 ♦ D 1 * D2 ) 

DER = DER - C 01 -‘02 ) * R(I,J) / < 01 * 02 ! 

R 1 ( I » J ) = - DER / ( R ( I 9 J ) * R ( I » J ) ) + R I ( I « J ) 

1 = MSTLP2 
J = N STL PI 

D 1 = DPH I (1-11 

02 = OPHI ( I ) 

DER ■= ( Q ( I * 1 , J ) - Q ( I f J ) ) / 02 

R I ( I # J ) = - DER / R( I , J ) 

DER •= D 1 * R(I+1,J) / ( 02 * 01 + D2 * D2 ) 

DER •= DER - 02 * R(I-1,J) / ( 01 * 01 ♦ 01 * 02 ) 

DER = DER - ( 01 - D2 ) * R(I,J) / ( 01 * D2 ) 

R 1 ( I » J ) = - DER / ( R ( I » J ) * R ( I t J ) ) + R 1 ( I t J ) 

80 DO 90 J=2,NM1 

ESIMP(J) = EPSI(J) 

FSIMP(J) = FPSI(J) 

90 G S I M P ( J ) = GPS I ( J ) 

DG 100 J = 2 1 NP 1 
100 DX( J ) = DPS I ( J ) 
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PAR2 116 
PAR2 117 
PAR2 118 
PAR 2 119 
PAR 2 120 
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IMINSI = K 

1 MAX S I = NP1 
DC 120 1=2, MP1 
DO 110 J = K , NP 1 

110 GRAND ( J ) = Rl< I , J ) 

IF (I .GT. MSTUP1) GR AND ( NSTUP2 ) = GRAND (NSTAGU ) 

CALL SIMP 

TEMP = ALPHA! I , K) 

DO 120 J=K,NP1 

12 0 ALPHA ( I , J I = GIFT(J) + TEMP 

COMPUTE ALPHAS IN UPPER SLOT 

DO 150 1=2, MM1 
E S I M P ( I ) = EPHI (I) 

FSIMP(I) = FPHI(I) 

15 0 GSl MP( I ) = GPHI ( I I 
DO 160 1=2, MPl 
160 DX( I ) = DPMI { I ) 

IMINSI = MSTUM1 
IMAXSI = MSLUP1 
J = NP 1 
D 1 = DPSI (N ) 

02 = D1 

DO 170 I=MSTUM1,MSLUPI 

O'ER = 01 * Q ( I , J+ 1 ) / ( 02 * Dl + D2 * D2 ) 

OER = DER - 02 * 0(1, J-l) / ( Di * D1 + D1 * D2 ) 

DER = OER - C 01 - 02 I * 0(1, J) / ( 01 * D2 ) 

170 GRAND ( I ) = R( I ,J) * DER 
CALL SIMP 

00 180 I=MSTAGU,MSLUP1 

18 0 ALPHA ( I , J ) = ALPHA ( MS TUM 1 f J ) + GIFT(I) 

1 = MSTUP2 
' J = NSTUP1 

Dl = DPMI (1-1) 

02 = OPHI ( I ) 

DER = (QSLOTU(Ul) - QSLOTU( I) ) / D2 
RKItJI = - OER / RSLOTU ( I ) 

DER = 01 * RSLOTU(H-l) / ( 02 * Dl ♦ D2 * D2 ) 

DER = DER - 02 * RSLOTU(I-l) / ( Dl * D 1 ♦ Dl * D2 ) . 

DER = DER - ( Dl - 02 ) * RSLOTUC I) / ( Dl * D2 ) 

R 1 ( I , J ) = - DER / ( RSLOTU(I) * RSLOTU( I ) ) +‘Rl(ItJ) 
DO 190 I=MSTUP3,MSLUP1 

01 = DPH I (1-1) 

02 = DPH I ( I ) 

DER = 01 * QSLOTUU + l) / ( D2 * Dl ♦ 02 * 02 ) 

DER = DER - 02 * QSLOTU(I-l) / ( 01 * 01 + Dl * 02 ) 

DER = DER - ( 01 - D2 I * USLOTIH II / < Dl * D2 ) 

R 1 ( I , J ) = - DER / RSLQTUC 1 ) 

DER = Dl * RSLOTU(Ul) / ( D2 * D 1 ♦ D2 ♦ D2 ) 

OER = DER - 02 * RSLGTU(I-l) / < Dl * Cl ♦ Dl * 02.) 
DER = DER - < 01 - D2 I * RSLOT U( l ) / ( Dl * D2 ) 

190 R1(I,J) = - DER / ( RSLOTU(I) * PSLOTU(I) ) ♦ R1(I,J) 
00 200 J = 2 , NM 1 
ESIMP( J) = EPS I ( J ) 

FSIMP( J) = FPSI ( J) 

20 C G S I M P ( J ) = GPS I ( J ) 

00 205 J=2 , N 
205 OX ( J ) -= DPSI ( J ) 

IMINSI = NSTUP1 
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I MAX S I = NP1 
DO 2 20 I = MS TAG U ? MSLUP1 
DO 210 J=N$TUP1,NP1 
210 GRAND ( J ) = Rl< I, J > 

CALL SIMP 

TEMP = ALPHA(1,NP1) - GIFT(NPl) 

00 215 J = NS TU P2 ,N P 1 

215 ALPHA!!, J) = GIFT(J) ♦ TEMP 
220 ASLOTUm = GIFT(NSTUPl) ♦ TEMP 

1 = MSTUP1 
J = NSTUP1 

ALPHA ( I , J) = A SLGTU ( I ) 

COMPUTE ALPHAS IN LOWER PORTION OF DIFFUSER 

IF ( FLAG1 .GT. 0 .OR. FLAG3 .GT. 0 ) GO TO 400 
DO 240 I =2 • M ST AGL 
DO 230 J = 2 r K 
230 GRAND ( J ) •= R1 ( I»J) 

KM 1 = K - 1 
I MINS ! = 2 
IMAXSI = K 
CALL SIMP 

TEMP = ALPHA(I f K) - GIFT(K) 

DO 240 J = 2 » K M 1 

240 ALPHA! I t J ) = GIFT ( J ) + TEMP 

1 = MSTLP1 

DO 260 J=NSTLP2,K 
260 GRANO( J ) = R1 ( I ,J) 

IMINSI = NSTLP2 
IMAXSI = K 
CALL SIMP 

TEMP = ALPHA(I,K) - GIFT ( K ) 

DO 270 J=NSTLP2tKMl 
270 ALPHA ( I , J ) -= GIFT(J) ♦ TEMP 
IMINSI = NSTLP1 
IMAXSI = K 
DO 300 I=MSTLP2»MP1 
00 290 J=NSTLP1»K 
290 GRAND! J) = Rl( I ,J) 

CALL SIMP 

TEMP = ALPHA! I , K ) - GIFT(K) 

DO 300 J=NSTLP1,KM1 
30 C AL PH A ( I 9 J ) = GIFT(J) ♦ TEMP 

COMPUTE ALPHAS INSIDE LOWER SLOT 

DO 310 1=2, MM1 
E SIMP! I ) = EPHI (I ) 

FSIMP! I ) = FPH HI) 

310 GSIMPi l ) = GPH I ( I ) 

DU 320 1=2, MPl 
32 C DX ( I ) = DPH I { I ) 

IMINSI = MSTLM1 
IMAXSI = MSLLP1 
J = 2 

D 1 = DPSI (2 ) 

02 = 01 

DO 3 30 I = M S T L M 1 1 M S L L P 1 

DER = 01 * Q ( I , J + 1 ) / ! C2 * Di ♦ 02 * D2) 
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DER = DER - 02 * Q(I,J-1) / ( 01 * D1 + 01 * D2 ) 

OER = DER - ( 01 - D2 ) * 0(1, J) / ( D1 * 02 ) 

330 GRAMOUI = RtltJ) * DER 
CALL SIMP 

DO 340 I=MSTAGL,MSLLP1 
340 ALPHA ( I , J ) = ALPHA(MSTLMI f J) ♦ GIFT(I) 

1 ^ KSTLP2 
J = NSTUPi 

Dl * DPHK 1-1) 

02 - DPHUII 

DER * ( QSLGTL ( 1 + 1 ) - Q SLOTH I ) ) / D2 
P1(I ,J> = - DER / RSLOTIU) 

DER - DI * RSLOTL(HU / ( D2 * Dl . ♦ 02 * 02 ) 

OER •-= OER - D2 * RSLOTL(I-l) / { 01 * 01 ♦ 01 * 02 ) 

OER = DER - ( Dl - D2 ) * RSLOTL(I) / ( 01 * 02 ) 

Rl(ItJ) = - OER / ( RSLOTL(I) * RSLCTL( I I ) + R1(I,J) 
. DO 330 I=KSTLP3,MSLLP1 

01 * OPHI ( I-l) 

02 - OPHI ( I ) 

OER - 01 * QSLOTL(Hl) / ( 02 * Dl + D2 * 02 ) 

DER - OER - D2 * QSLOTL(I-l) / ( Dl * Dl * Dl * 02 > 

OER = OER - < Dl - 02 > * OSLOTL(I) / ( 01 * 02 ) 

R1C I , J) = - OER / R SLOTL ( I ) 

OER = 01 * RSLOTL(H-i) / ( 02 * Dl + 02 * 02 ) 

OER = DER - D2 * RSLOTHI-l) / ( 01 * 01 + 01 * 02 ) 

: OER = DER - ( 01 - 02 ) * RSLUTL(I) / ( 01 * D2 > 

350 Rl(ftJ) - - OER / ( R SLOTL ( I ) * RSLOTHII ) + R1(I,J) 
on 360 J = 2 , NM 1 
esiMP(j) = EPSim 
FSIMP(J) = FPSIU) 

360 GSIMP(J) = GPS I(J) 

.00 370 J~2,N 
370 DX( J ) = OPS I < J ) 

. I MI MSI = 2 , 

>; I.MAXSI = MSTLP1 
,00 390 I=MSTAGL,MSLLP1 
DO 330 J=2,NSTLPI 
380 GRAND ( J ) = Rl< I,J) 

CALL SIMP 

. DO 385 J = 3 , NST AGL 

385 ALPHA(I,J) = A L PH A (1,2) + GIFT(J) 

390 ASLOTL(I) = GIFT(NSTLPl) + ALPHA(I,2) 

[ ^ MSTAGL 
J = NSTLPI 

ALPHAU , J) = A SLOTL ( I) 

ALPH A( HI , J ) = AL PH A ( I , J) 

40 C CONTINUE 

COMPUTE INLET R COOROI MATES 

DO 410 J = 2 , NM1 
ESIMP(J) = EPS I C J I 
FSIMP(J) - FPSI(J) 

410 GSIMP(J) = GPS I ( J> 

DO 420 J=2 , N 
420 DX ( J ) = OPS I ( J ) 

IMINSi = 2 
I MAX SI = NP1 

IF (RAOIN .GT. 0.0) R(2,2) = RAOIN 

IF (FLAGI .GT. 0 .OR. FLAG 3 .GT. 0) R(2,2) = 0.0 
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DO 430 J=2,NP1 

PARZ 

339 


P= 0(2, J) 

PAR 2 

340 

43 C 

GRAND ( J ) = 2. * COS ( ALPHA ( 2 , J ) I / DEXP(P) 

PAR2 

341 


CALL SIMP 

PAR2 

342 


POW = .5 

PAR2 

343 


DIV = 1.0 . 

PAR2 

344 


IF ( FLAG1 ,GT. 0) POW = 1.0 

PAR2 

345 


IF ( FLAG1 . GT. 01 DIV = 2.0 

PAR2 

346 


DO 440 J=3, NP1 

PAR2 

347 


RAT = ABS ( GIFT ( J ) / DIV ♦ R (2,2) * R(2,2) I 

PAR2 

348 

44 C 

R ( 2 , J I = ( RAT )**P0W 

PAR2 

349 



PAR2 

350 


COMPUTE R ALONG MIDDLE STREAMLINE 

PAR2 

351 



PAR2 

352 


DO 450 1=2, MM1 

PAR2 

353 


ES I MP ( I ) = EPH K I ) 

PAR2 

354 


FSIMP(I) = FPHI(I) 

PAR2 

35 5 

450 

GSIMP( I ) = GPH [III 

PAR2 

356 


DO 460 1=2, M 

PAR2 

357 

46 0 

DX( I ) = D PH I ( I I 

PAR2 

358 


I MI NS I = 2 

PAR2 

359 


IMAXSI = MP 1 

PAR2 

360 


J = K 

PAR2 

361 


DO 470 1=2, MP1 

PAR2 

362 


P = Q( I , J » 

PAR2 

363 

470 

GRAND ( l ) = S I N ( AL PH A( I , J ) ) / DEXP(P! 

PAR2 

364 


CALL SIMP 

PAR2 

365 


DO 490 1=3, MP1 

PARZ 

366 

480 

R U ,J) = R( 2 » J ) + GIFT ( I > 

PARZ 

367 



PARZ 

368 


COMPUTE R THROUGHOUT UPPER PORTION OF DIFFUSER 

PAR 2 

369 
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370 


DO 490 J = 2 , NM 1 

PARZ 

371 


E S I MP ( J I = EPS I ( J I 

PAR2 

372 


FSIMP(J) = FPSHJI 

PAR2 

373 

49 C 

GSIMP(J) = GPSI(J) 

PAR2 

374 


DO 500 J=2 , N 

PAR2 

375 

50 C 

DX( J 1 = DPSI ( J I 

PARZ 

376 


IMINSI = K 

PAR2 

377 


IMAXSI = NPl 

PAR2 

378 


DO 520 1=3, MP1 

PARZ 

379 


00 510 J = K , N P 1 

PAR2 

380 


P = Q( I , J) 

PAR2 

381 

510 

GRAND ( J ) = 2. * COS ( ALPHA( I , J I ) / DEXP(P» 

PAR2 

382 


IF (I .GT. MSTUP1) GRAND (NSTUP2 ) = GRAND ( NSTAGU ) 

PAR2 

383 


CALL SIMP 

PAR2 

384 


KP1 = K *■ 1 

PAR2 

385 


DO 520 J= K PI , NPl 

PAR2 

386 


RAT = AB S ( GIFT(J) / DIV + R(I,K> * R ( I , K » ) 

PAR2 

387 

520 

R ( I , J » = ( RAT )**POW 

PAR2 

388 



PAR2 

389 


COMPUTE R COORDINATES INSIDE UPPER SLOT 

PAR2 

390 



PAR2 

391 

56C 

DO 570 1=2, MM1 

PAR2 

392 


ESIMP(I) = EPH I ( I I 

PAR2 

393 


FSIMP(I) = F PH I ( I ) 

PAR2 

394 

570 

GSIMP(I) = GPH 1(1) 

PAR2 

39 5 


DO 590 1=2, M 

PAR2 

396 

580 

DX ( I ) = DPH I ( I ) 

PAR2 

397 


IMINSI = MSTUM1 

PAR2 

398 
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I MAX S I = MSLUPl 
J = NP1 

DO 590 I=MSTUM1, MSLUPl 
P - Q( I • J) 

5 9 C GRAND ( I ) = S I N ( AL PH A < I , J ) ) / OEXP(P) 

CALL SIMP 

DO 600 I=MSTAGU»MSLUP1 
600 R ( I t J ) = R ( MST UM 1 » J ) ♦ GIFT(I) 

DO 605 J = 2 * NM1 
ESIMPCJ) = E PS I ( J I 
t FS I MP ( J ) = FPSI(J) 

605 GSIMPC J) = GPS I ( J) 

DO 610 J = 2 t N 
610 DX ( J ) = DPSI (J ) 

I MI NS I = NS TUP 1 
I MAX S I •= NP1 
DO 650 I=MSTAGUt MSLUP1 
DO 620 J = NS TUP 2 » NP 1 
P ^ Q(ItJ) 

620 GRAND ( J ) = 2, * CO S ( ALPH A ( I # J ) ) / DEXP(P) 

P = QSLOTU ( I ) 

GRAND ( NSTUP 1 ) = 2. * COSCASLGTUUH / DEXP(P) 

CALL SIMP 

TEMP = R(ItNPl) * RIItNPlJ - GIFTINPl) 

DO 630 J = NSTUP2t N 

630 R ( I t J ) = ( ABS ( GIFT(J) + TEMP ) )**.5 

RSLOTUU) = ( ABS ( GIFT(NSTUPl) ♦ TEMP ) )**.5 

I F ( FLAG1 .LE. 0) GO TO 650 
TEMP = R ( I t NP1 ) - GIFTCNP1) 7 2. 

DO 640 J = NS TUP 2 t N 
640 R(I,J) = GIFT(J) / 2. ♦ TEMP 

RSLOTUtl) ^ GIFT(NSTUPl) / 2. ♦ TEMP 
. 650 CONTINUE 

I = MS TUP 1 
J = NSTUP1 
P = Q( I-1,J> 

R ( I t J ) = R ( I- 1 t J ) + SIN(ALPHA( I-lf J)> * DPHIU-1J / DEXP(P) 

COMPUTE R OVER LOWER PORTION OF DIFFUSER 

IF ( FLAG1 • GT. 0 .OR. FLAG3 .GT. 0) GO TO 840 
DO 670 I = 2 t MST AGL 
DO 655 J = 2 f K 
P = Q ( I , J ) 

655 GRAND ( J ) = 2. * COS ( ALPH A < I , J )) / OEXP(P) 

IMINSI = 2 
I MAX SI = K 
CALL SIMP 

TEMP = R ( I » K ) * R ( I » K ) - GIFT ( K ) 

DO 660 J = 2 t KM1 

660 R ( I f J ) = ( ABS ( GIFT(J) + TEMP ) >**.5 
670 CONTINUE 
I = MSTLP1 
DO 675 J = NS TLP 2 » K 
P = Q(IiJ) 

675 GRANO(J) = 2. * COS ( ALPHA! I * J ) ) / DEXP(P) 

IMINSI = NSTLP2 
I MAX S I = K 
CALL SIMP 

TEMP = R ( 1 1 K I * R ( I » K ) - GIFT(K) 
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00 680 J = NS TLP 2 , KM 1 

63 C R(IJ) = ( ABS< G I F T C J I + TEMP ) )**.5 

690 CONTINUE 

00 7C2 I=MSTLP2»MP1 

00 695 J = NS TLP 1 , K 

P = Q( I ? J ) 

695 .GRANO(J) = 2. * COS ( ALPHA! I , J ) ) / DEXP(P) 

IMINSI = NS TLP 1 
IMAXSI = K 
CALL SIMP 

TEMP = R(I,K> * R ( I f K ) - GIFT(K) 

DO 697 J = NS TLP 1 f KM 1 

697 K ( I f J ) = ( A3S ( GIET(J) + TEMP ) )**.5 
702 CONTINUE 

COMPUTE R INSIDE LOWER SLOTS 

DO 740 1=2, MMI 
ESIMPU ) = EPHKI ) 

FSIMP( l ) = FPH 1(1) 

740 GSI M P ( I ) = GPH I ( I ) 

DO 750 1=2, M 

75 C DX( I) = OPHI ( I ) 

DO 760 I=MSTLM1,MSLLP1 
P = Q( I ,2 ) 

76 C GRAND ( I ) = S I N ( AL PH A ( l , 2 ) ) / OEXP(P) 

IMINSI = MSTLM1 

IMAXSI = MSLLP1 

CALL SIMP 

DO 770 I=MSTAGL , MSLLP1 

77 C R ( I , 2 ) « R < MST LM1 , 2 ) ♦ GIFT(I) 

DO 780 J=2 , NMI 
ESIMPU) = EPSI(J) 

FSIMP(J) = FPS I ( J ) 

78 C GSIMP(J) = GPSI(J) 

DO 790 J = 2 , N 
790 DX( J ) = DPS I ( J ) 

IMINSI = 2 

IMAXSI = NSTLP1 

DO 830 I=MSTAGL,MSLLPI 
DO 800 J = 2 , NST AGL 
P = Q( I , J ) 

800 GRANO(J) = 2. * COS ( ALPHA! I , J ) ) / DEXP(P) 

P * OSLOTL(I) 

GRAND(NSTLPl) = 2. * COS ( ASLGTL ( I) ) / DEXP(P) 

CALL SIMP 

DO 810 J = 3 , NST AGL 

RAT = ABS( GIFT(J) + R C 1,2) * R(I,2) ) 

810 R( I , J) = ( RAT |**.5 

RAT = A9S( GIFT ( NSTLP 1 ) + R(I,2) * R(I,2) ) 

RSLOTL ( I ) = ( RAT )**.5 

IF ( FLAG1 .GT. 0) GO TO 830 
DO 820 U=3,NSTAGL 

82 C P(I,J) = R ( 1,2) + GIFT(J) / 2. 

RSLOTL(I) = R ( 1,2) ♦ GI F T{ NSTLP 1 ) / 2. 

830 CONTINUE 

1 = MSTLP1 
J = NSTLP1 

P = Q t I - 1 , J ) 

R ( I , J ) = R ( 1-1 * J ) + SIN( ALPHAC 1-1, J) ) * OPHKI-l) / DEXP(P) 
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R( I-lt Jl = R S L C T L ( 1-1 ) 

840 CONTINUE 

COMPUTE CONVERGENCE CRITERION 
REWINO 5 

READ (5I( (R1 U , J) • 1 = 1 »MP2 ) , J=1 ,NP2 ) 

EPS = 0.0 
DO 845 J =2 , NP 1 
DO 845 1=2, MP1 
EPSTR = ABS(RI ( I, J) - 
IF (EPSTR .LE. EPS) GO TO 345 
EPS = EPSTR 
845 CONTINUE 

COMPUTE INLET X COORDINATES 

IF (ITER .LT. MXISYS .AND. EPS .GT. TOLSYS) GO TO 3000 
DO 848 1=1, MP1 
- X SLOTU ( I > = 1. 0 
848 XSt 0 TL ( I ) = I. 0 
00 850 J=2 , NM1 
ESIMP(J) = EPSI(J) 

850 GSIMP(J) = GPS I ( J ) 

FSIMP(J) = FPSI(J) 

00 860 J = 2 , N 
860 0X( J ) = DPSI(J) 

IMINSI = 2 

1 MAX S I = NP 1 
JOT = 2 

IF ( FLAG 3 .GT. 0) JOT = 3 
DO 870 J = JOT , N P 1 

P = 0(2tJ) > 

GRAN D ( J ) = - SIN( ALPHA! 2t J) I / DEXP(P> 

870 IF ( FLAG1 .LE. 0) GRAND ( J ) = GRAND ( J ) / R(2,J) 

IF (FLAG3 .GT. 0) GRAND ( 2 ) = 0.0 
CALL SIMP 
DO 880 J=2 , NPI 
880 X ( 2 , J ) = GI FT ( J I 

COMPUTE X ALONG THE MIDDLE STREAMLINE 

DO 890 1=2, MM1 
ESIMP(I) = EPHI(I) 

F S I M P ( I ) = FPHI (I ) 

890 GS I M P ( I ) = GPHI ( I ) 

DO 900 1=2, M 
900 DX( I ) = DPH I ( I ) 

IMINSI = 2 
IMAXSI = MP1 
J = K 

DO 910 1=2, MP1 
P = 0( I , J) 

910 GRANDU) = COS ( AL PH A ( I , J ) ) / DEXP(P) 

CALL SIMP 
DO 920 1=3, MP1 

920 X( I , J) = X( 2,J ) ♦ G I F T ( I ) 

COMPUTE X THROUGHOUT UPPER PORTION OF DIFFUSER 
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PAR2 576 
PAR2 577 
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00 930 J=2,NM1 
E S I M P < J ) = EPSKJ) 

FSIMP(J) = FPSKJ) 

930 GS I M P ( J ) = GPS I ( J I 
00 940 J= 2 , N 
94 C 0X( J ) - DPSKJ ) 

IMINSI = K 
IMAXSI = NPI 
JET = K 

IF ( FLAG3 .GT. 0) JET = 3 
00 960 1 = 3, MPI 
0Q 950 J= JET, NPI 
P = G( I , J) 

GRAND ( J ) = - S I N( ALPHA ( I , J H / DEXP(P) 

950 IF ( FLAG1 .LE. 0) GRAND ( J ) = GRAND(J) / RU,J) 

IF { FLAG3 .GT. 01 GRAND ( 2 ) = 0.0 

IF (I .GT. MSTUP1) GRAND ( NSTUP2 ) = GR ANO ( N ST AGU ) 
CALL SIMP 
KP1 = K + 1 
DO 960 J=KP1,NP1 
960 X( I , J ) = X(I ,K ) + GIFT! J) 

COMPUTE X COORDINATES INSIDE UPPER SLOT 

00 990 1=2, MM1 
ESIMPH ) = EPHI l I ) 

FSIMP(I) = FPHIU) 

99 C G S I M P( I ) = GPHI (I ) 

00 1000 1=2, M 
100C DXUI a DPH I ( l ) 

IMINSI = MSTUM1 
IMAXSI = MSLUPl 
J = NPI 

DO 1010 I=MSTUM1, MSLUPl 
P = Q( I , J) 

1010 GRAND ( I ) = COS (ALPHA! I ,J )) / OEXP(P) 

CALL SIMP 

• DO 1020 I=MSTAGU, MSLUPl 
1020 X(I,J) = X ( MST UM 1 1 J ) + G I FT ( I I 
DO 1030 J = 2 , NM 1 
ESIMP(J) = EPSKJ) 

FSIMP(J) = FPSI(J) 

1030 GS I MP ( J ) = GPS I( J ) 

DO 1035 J = 2 , N 
1035 OX ( J ) = DPSKJ ) 

IMINSI = NS TUP 1 

IMAXSI = NPI 

00 1060 I=MSTAGU, MSLUPl 

00 1040 J = N STU P2 , NP 1 
P = Q( I , J) 

GRANO(J) = - SIN(ALPHA( I,J) ) / DEXP(P) 

1040 IF (FLAG1 .LE. 0) GRAND ( J ) = GRAND ( J ) / R(I,J) 

P = OSL OTU ( I ) 

GRAND ( NS TUP 1 ) = - S I N ( A SLQTU( I ) ) / DEXP(P) 

IF ( FLAG1 .LE. 0 ) GRAND(NSTUPl) = GR AND ( NSTUP 1 ) / 
1 R SLO TU ( I ) 

CALL SIMP 

TEMP = X ( I , N PI ) - G I FT ( NPI ) 

DO 1050 J=NSTUP2,N 
1050 X( I f J) = GI FT ( J ) + TEMP 
1060 XSLOTU(I) = GIFT(NSTUPi) * TEMP 

1 = MSTUPl 
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1065 

1 07 C 
1075 

1080 

1085 

109C 

1120 

1130 

1140 

1 150 

1160 

1170 


J = NSTUP1 






PAR2 

63 9 

P = O(I-lfJ) 






PAR2 

640 

XUtJJ = X(I-1,J) + COS ( ALPHA ( I - 

1, J ) > * 

DPHUI-1) / DEXP ( P \ 

PAR2 

641 







PAR2 

642 

COMPUTE X OVER LOWER PORTION 

OF DIFFUSER 

PAR2 

643 







PAR2 

644 

IF ( FLAGI .GT. 0 .OR. FLAG3 

• GT 

. 0) GG 

TO 

3000 

PAR2 

64 5 

DU 1070 I =2 , MS T AGL 






PAR 2 

646 

DO 1065 J = 2 , K 






PAR2 

647 

P = G(l,J> 






PAR2 

648 

GRAND ( J) = - S INI ALPHA ( I * J ) ) 

/ i 

DEXPI P> 

* 

R(I,J)I 

PAR2 

649 

l MI NS I = 2 






PAR2 

650 

I MAX S I = K 






PAR2 

651 

CALL SIMP 






PAR2 

652 

TEMP * X C I , K 1 - GIFT ( K ) 






PAR2 

653 

DO 1070 J=2,KM1 






PAR2 

654 

X( I , J) = GIFT! J ) ♦ TEMP 






PAR2 

655 

I = MSTLP1 






PAR2 

656 

DO 1075 J =NSTLP2,NP1 






PAR2 

657 

P = G < I * J > 






PAR2 

65 8 

GRAND ( J ) = - SIN! ALPHA! 1 ,J>> 

/ i 

DEXPl P) 

* 

RII , J) ) 

PAR2 

659 

IMINSI = NS TIP 2 






PAR2 

660 

IMA X SI = K 






PAR2 

661 

CALL SIMP 






PAR 2 

662 

TEMP = X< I , K ) - G I FT ( K ) 






PAR2 

663 

DO 1080 J=N STL P2» KM 1 






PAR2 

664 

X( I , J) = G I F T( J ) + TEMP 






PAR2 

665 

DO 1090 I=MSTLP2,MP1 






PAR2 

666 

DO 1085 J=NSTLP1,K 






PAR2 

667 

P = Q( I , J) 






PAR2 

668 

GRAND ( J ) * - SIN! ALPHA ( I » J ) I 

/ ( 

DEXPCP) 

* 

R ( I , J > > 

PAR2 

669 

IMINSI = NSTLP1 






PAR2 

670 

I MAX S I = K 






PAR2 

671 

CALL SIMP 






PAR2 

672 

TEMP = X( I , K) - GIFT(K) 






PAR2 

673 

DO 1090 J=NSTLP1,KM1 






PAR2 

674 

X( I , J) = GIFT! J) + TEMP 






PAR 2 

675 







PAR2 

676 

COMPUTE X INSIDE LOWER 

SLOTS 




PAR2 

677 







PAR2 

678 

DO 1120 I = 2 , MM 1 






PAR2 

679 

E S IMP( I ) = EPHI ( I ) 






PAR2 

680 

FS I M P< 1 1 =' FPH 1(1) 






PAR2 

681 

GSI M P ( I ) = --GPHI ( I ) 






PAR2 

682 

DO 1130 1=2, M 






PAR2 

683 

DX( I ) = DPHII I ) 






PAR2 

684 

DO 1140 I=MSTLM1,MSLLP1 






PAR2 

685 

P = Q ( I , 2 ) 






PAR2 

686 

GRAND! I) = COS ( AL PH A I 1,2) ) 

/ 

OEXP(P) 



PAR2 

687 

IMINSI = MSTLM1 






PAR2 

688 

IMAXSI = MSLLP1 






PAR2 

689 

CALL SIMP 






PAR2 

690 

DO 1150 l=MSTAGL,MSLLPl 






PAR2 

691 

X( I f 2) = X( MSTLM1 . 2 ) + GIFT(I) 




PAR2 

692 

DO 1160 J-2,NH1 






PAR2 

693 

ESIMP! Jl * EPSI ( J) 






PAR2 

694 

FSIMP(J) = FPSI(J) 






PAR2 

695 

GSIMP(J) = GPSHJI 






PAR2 

696 

DO 1170 J = 2 » N • 






PAR2 

697 

DX ( J ) = DPS I ( J ) 






PAR2 

698 


o no 


IMINSI = 2 
IMAXSI = NSTLP1 
DO 1200 I=MSTAGL,MSLLP1 
DO 1 1 BO J = 2 , NS T AGL 
P = Q( I t J ) 

GRAN D ( J ) * - S IN! ALPHA { I , J ) ) / DEXP(P) 

1180 IF (FLAG1 .LE. 0) GRAND(J) = GRAND ( J ) / R(I,J) 

P = QSLOT L ( I ) 

GRAND ( NS TLP 1 ) - - S I N ( ASLOTL ( I) ) / DEXP(P) 

IF !FLAG1 .IE. 0) GRAND! NSTLP1 ) = GRAND! NSTLP1 ) / 

1RSLOTL ( I ) 

CALL SIMP 
DO 1190 J = 3 , NS TAGL 
1190 X(ItJ) = X ( 1,2) + GIFT ( J ) 

1200 XSLOTL(I) = X < I » 2 I ♦ GI FT! NSTLP1 ) 

1210 I = MSTLP1 
J = NSTLP1 
P = QU-ltJ) 

X ! I , J ) = X(I-1,J) + COS! ALPHA! I — 1 , J ) ) * OPHI(I-l) / DEXP(P) 
3000 CONTINUE 

ASLOTU! MSTUP1) = ALPHA! MSTUP1,NSTUP1 ) 

RSLOTU(MSTUPl) = R ! MSTUP1 , NSTUP1 ) 

XSLOTU(MSTUPl) = X ( MSTU PI , NSTUP1 ) 

ASLOTL! MSTL PI ) = ALPHA! MSTLP1 , NSTLP1) 

RSLOTL! MSTLP1) = R ( MSTL PI , NSTLP 1 ) 

XSL0TL!MSTLP1) = XIMSTLPltNSTLPl ) 

SET UP BOUNDARY CONDITIONS FOR NEXT ITERATION 

IOT = MP1 

IF < FLAG1 .LE. 0 .AND. FLAG3 .LE. 0) IDT = MSLLP1 
DO 2010 1-2, IOT 
Y3 = R( I ,4) 

Y2 = R ( I , 3 ) 

Y 1 = R( I ,2) 

D1 * DP S I ( 2 ) 

02 = DPS 1(3) 

DER = ( Y2 - Y1 M ( D1 + 02 ) / ( D1 * 02 ) 

DER = DER - ( Y3 - Y1 ) * D1 / ( D1 * D2 «- 02 * D2 ) 

2010 R ( I , 1 ) = Y2 - 2. * DER * D1 
IOT = MSLUP1 
DO 2020 1 = 2, IOT 
Y3 = R( I , NP1 ) 

Y2 = R ( I , N ) 

Y 1 = R ( I , NM 1 1 

01 = DPSI ( NM1 ) 

02 = OPS I ( N I 

OER = ( Y3 - Y2 ) * ( 01 ♦ D2 ) / ( 01 * D2 ) 

DER = OER - ( Y3 - Y1 > * 02 / ( 01 * 01 ♦ 01 * D2 • 

2 02 C R(I f NP2l = Y2 «■ 2. * DER * D2 
DO 2030 J = 2 , NP 1 

Y 3 = R(4,J) 

Y 2 = R( 3, J) 

Y1 = R ( 2 , J ) 

D 1 = OPHI (2) 

D2 = DPHI (3) 

DER = ( Y2 - Y1 ) * ( D1 + D2 ) / ( 01 * D2 ) 

DER = DER - ( Y3 — Y1 ) * D1 / ( D1 * 02 + 02 * D2 ) 

2 03 C R ( 1 , J ) = Y2 - 2. * DER * 01 
DO 2035 J=2 , NP 1 

Y 3 = R(MPltJ) 
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Y2 = RCM, J) PAR2 

Y 1 = R ( MM 1 , J ) PAR2 

Dl = DPHKMM1) PAR2 

D2 = DPHI(M) PAR2 

DER = ( Y 3 - Y2 ) * ( Dl ♦ D2 I / ( ' Dl * D2 ) PAR 2 

DER = DER - ( Y 3 — Y 1 ) * D2 / ( Dl * Dl + 01 * 02 ) PAR2 

2035 R ( MP2 1 J ) = Y2 + 2. * DER * 02 PAR2 

IF ( FLAG2 .LE. 0) GO TO 2050 PAR2 

DO 2040 J = 2 * NP 1 PAR2 

R ( 1 f J ) = R ( 3 t J I PAR2 

204 C R ( MP 2 , J ) = R(M,J) PAR2 

2050 CONTINUE PAR2 

PAR2 

WRITE CONVERGENCE VARIABLES PAR2 

PAR2 

IE ( FLAG 1 .LE. 0) WRI T E ( 3 , 2060 II TER ♦ JTER ,EPS, DEL PAR2 

2060 F0RMAT(/T2, * ITER = * , T9 , 1 5 ♦ T25 , • JTER = 1 * T32 , I 5 , T50 , PAR2 

1'EPS = * ,T56,E13.6,T75,*DEL = *tT80 v E13.6) P AR2 

PAR2 

CHECK BELL CHANNEL SOLUTION PAR2 

PAR2 

IE ( FLAG3 .LE. 0) GO TO 2130 P AR2 

IE (EPS .GT. TOLSYS .AND. ITER .LT. MXISYS) GO TO 2130 PAR2 

IM1NSU •= 2 PAR2 

I MAX SU * NS T AG U - 1 PAR2 

1 MINS I = 2 PA R 2 

IMAXSI = NS TAG U ♦ 1 PAR2 

ERROR » 0.0 PAR2 

‘ 00 2090 I =2 » MP 1 PAR2 

00 2070 J-2 tNS TAGU PAR2 

207 C OX ( J ) = SORT ( ( ABS( X(i,J+l) - X ( I , J 1 ) )**2. PAR 2 

1* ( AB S ( R ( I , J + 1 ) - R ( I v J ) ) >**2. ) PAR2 

CALL SETUP PAR2 

NSTUP1 = NS T AGU ♦ 1 PAR2 

00 2080 J=2,NSTUP1 PAR2 

2 08 C GRAND ( J ) = R ( I ♦ J I * EXP(Q(J,J)) PAR 2 

CALL SIMP PAR2 

ERR = AB S ( PSI(NSTUPl) - GIFT(NSTUPl) ) PAR2 

IF ( ERR .GT. ERROR ) FNC - GIFT(NSTUPl) PAR2 

2090 IF ( ERR .GT. ERROR ) ERROR = ERR PAR2 

PC E = 100. * ERROR / PS I(NSTUPl) PAR2 

WRITE (3 , 2100 )HEUR , DA TE1 ,DATE2, DAT E3, CASE PAR 2 

2 10 0 FORM AT ( 1H1//T57, * CLEMS ON UNIVERSITY* , //T49 , • MECHAN I C AL E NG I NEERI P AR2 
1NG DEPARTMENT # »///T2, 15 A4,T 64 tI2»T 66 f l / , tT67 f I2,T69, •/• ,770,12 t PAR 2 
2T 76 , 'CASE NO. *,187,14///) PAR2 

WRITE(3, 2110) PAR2 

2110 FORMAT ( T2 , * THE ACCURACY OF THE SOLUTION FOR THE*, PAR 2 

1 T 39 , • BELL CHANNEL DESIGNED*) PAR2 

WRI TE ( 3 ,2120) ERROR, FNC , PC E PAR 2 

2 1 2 C FORM AT ( ////T2, * ERROR = *,F10.6//// PAR2 

1T2, * INCORRECT VALUE = *,F10.6//// PAR2 

2T2 , * PERCENT ERROR = *,F10.6) PAR2 

2 1 3 C CONTINUE PAR2 

REWIND 5 PAR2 

REWIND 12 PAR2 

WRITE (5) (<R(I, J) , I=l,MP2l f J=1,NP2) PAR 2 

WRITE! 11 ) ( ( ALPHA ( I , J) ,I=l f MP2) , J=l, NP2) PAR2 

WRITE(12)((X(I,J), 1=1, MP2) , J= 1 » NP2 ) PAR2 

RETURN PAR2 

’ END PAR2 
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********** SUBROUTINE OUTPUT ********** 

THIS SUBROUTINE PROVIOES PRINTER OUTPUT OF THE 
RESULTS OF THE MAIN PROGRAM 


SUBROUTINE OUTPUT 

DIMENSION DX< 150) , ESIMPl 150) ,FSIMP( 150) ,GSIMP( 150) , 

1 GRAND (150) ,GIFT< 150) ,PHI< 150) ,PSI ( 150) , 

2 DP SI ( 150) , DPH I ( 150) ,S< 150) ,SBAR(150) , 

3 SI (150) ,SBAR1 ( 150) , S2( 150) ,S8AR2< 150) , 

HEDR( 15) , EPHI ( 150) ,FPHI( 150) ,GPHI ( 150) , 

EP S I ( 150) ,ASL0TL( 150) ,XSL0TU( 150) ,XSLGTL( 150) 
,F PS I (150) ,GPSI (150) ,RSL0TU(150> ,RSL0TL( 150) , 
QSLOTU! 150) ,USL0TL( 1 50 ) , AS LOTU ( 1 50 > 

DIMENSION Q ( 150,50 ) ,R(150,50),X( 1 50, 50), ALPHA! 150,50) 
DOUBLE PRECISION P 

INTEGER FLAGi, FLAG2, FL AG3, FLAG* , C AS E , DATE 1 , DA TE2 , DA TE3 
COMMON / C2 / FLAGI ,FLAG2,FLAG3 .FLAG*, HEOR.OATEl ,0ATE2 
, DATE 3, CASE, MX IDE , ICKDE, IPRDE , IMIN, I MAX, 
JMIN, JMAX,NUT , M ,N , TOLOE , DEL , EPS , JTER, 

I TER, DPH I ,DPSI ,S, SBAR , RMUL T , SI ,S2,S8AR1, 
SBAR2,MSTAGU,NSTAGU,MSL0TU,MSTAGL,NSTAGL 
, MSLOTL, RADI N,W, PHI ,PSI , RSLOTU , RSLO TL , 
QSLOTU,QSLOTL,ASLOTU, ASLOTL , XSLOTU , 
XSLOTL , G , MI DJ 

CHANGE Q BACK TO VELOCITY AND OUTPUT 


MP1 = M + 1 

NP1 = N ♦ 1 

MM 1 = M - 1 

NM1 = N - 1 

MP2 = M + 2 

NP2 = N 2 

NSTLP1 = NSTAGL + 1 

NSTUP1 = NSTAGU + 1 

MSTUP1 = MSTAGU + 1 

MSTLP1 = MSTAGL * 1 

MSLLP1 = MSLOTL + 1 

MSLUP1 = MSLOTU + 1 

REWIND * 

REWIND 5 
REWIND 11 
REWIND 12 

READ (*) ( ( 0 C I , J ) , I=1,MP2 ) ,J=1,NP2) 
READ(5)((R(I,J),I=1,MP2),J=1,NP2) 

READ! 11 ) ( (ALPHA! I , J ) , I = 1 , MP2 ) , J = 1 , NP2 ) 
READ (12 ) ( (X! I , J ) ,1 = 1, MP2) ,J = i»NP2) 

IF ( NUT .EQ. 1 ) GO TO 1000 
00 10 1=2, MP1 
DO 10 J=2,NP1 
P = 0! I , J ) 

10 0! I,J> = DEXP(P) 

0( MSTUP1 , NS TUP 1 ) = 0.0 
Q(MSTLP1,NSTLP1 ) = 0.0 


OUTPUT Q ARRAY 


WRITE (3,20) HEDR,DATE1 ,DATE2 , DAT E3 , CASE 
20 FORMAT ( 1H1//T5 7 , 'CLEM SON UN I VER S I T Y • , // T49 , 
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OUT 

21 

OUT 

22 

OUT 

23 

OUT 

24 

OUT 

25 

OUT 

26 

OUT 

27 

OUT 

28 
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1 ' MECHANI CAL ENGINEERING DEPARTMENT* ,///T2,15A4,T64, 
2I2.T66, •/• ,T67, 12,T69, • /• ,T70, I2.T76, 'CASE NO. *, 
3T87,I4) 

WRITEI3, 30)1 TER, JTER, EPS, DEL 
3 C FORMAT ( ///T 1 , ' OVELOCI TY THROUGHOUT DIFFUSER*, 

1 T46 * • I TE R = ' , T53 , I4,T61 ,' JTER = • , T68 , I4,T76, • EPS = • 
2,T82,E13.6, T99, ‘DEL = * , T IC 5, E 1 3 . 6 , / // T7, * PHI • , T66 , 

3 * PS I ' ) 

WRITE! 3,40) (PS I! IP ) , IP = 2,NP1 ) 

90 FORMAT! //I13X, 5(7X,E13. 6) ) ) 

DO 50 1=2, MP1 

5 C WRI TE!3,60)PHI ( I ) , (Q( I , J) , J=2,NP1 ) 

toO FORMAT! //LX,E13.6,5(7X,E13.6),/(14X,5(7X,E13.6))) 

CHANGE Q BACK TO LOGARITHMIC FORM 

QIMSTUP1 ,NSTIJP 1 ) = 1.0 
Q(MSTLPl,NSTLPir = 1.0 
DO 70 1=2, MP1 
DO 70 J=2 , NP1 
P = Q ( I , J ) 

70 Q ( I , J ) = DLOG(P) 

OUTPUT VELOCITY IN SLOTS 

DO 75 1=2, MPl 
P = OSLOTL! I ) 

QSLOTL! I ) = DEXP(P) 

P = OSLOTU ( I) 

75 OSLOTU! I ) = DEXP! P) 

QSLOTU(MSTUPl) = 0.0 
OSLCTLt MSTLP1) = 0.0 

WRITE! 3,20)HEDR,DATE1,DATE2,DATE3,CASE 
WRITE (3, 80) ITER,JTER,EPS,DEL 
80 FORMAT! ///T2, * VELOCITY INSIDE SLOTS', T46, 

1 ' ITER = • ,T53, I4.T61,* JTER = • , T 68 , 1 4 , T76 , • EPS * *, 
2T82,E13.6,T99, 'DEL = • , T 105 , El 3 . 6 , / / /T 7 , • PS I • , T66 , 

3 ' PH I ’ ) 

WRITE! 3,40) (PHI! IP), IP =2, MPl) 

IF (FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) GO TO 90 
WR I T E < 3 , 85 ) 

85 FORMAT! ///T2, ' LOWER SLOT') 

WRI TE(3,60> PS I (NSTLP1 ) , ! OSLOTL < IP ) , IP = 2, MPl ) 

9C WRI TEI3, 100) 

IOC FORMAT! ///T2, ' UPPER SLOT') 

WRITE! 3, 60) PS I ( NSTUP1 ) , (OSLOTU! I P ) , I P=2 , MPl » 

OSLOTU! MSTUP1) = 1.0 
OSLOTL! MSTLP1 ( = 1.0 
DO 105 1=2, NP1 
P = QSLOTL! I ) 

QSLOTL ( I ) = DLOG(P) 

P = OSLOTU! I ) 

105 OSLOTU! I ) =. DLOG!P) 

OUTPUT OF R-COORDINATES < Y-COORD INATES FOR 2-D ) 

1000 WRI TE ( 3, 20)HEDR,DATE1 ,DATE2,DATE3,CASE 
WRITEOt 110) ITER, JTER, EPS, DEL 
110 F0RMAT(///T1,'0R - COORDINATES ! UNTRANSFORMED )', 
1T46,*ITER = • ,T53, I4,T61,* JTER = * , T68 , I 4 , T 76 , • E PS = * 
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2,T82,E13.6,T99, 'DEL = ■ , T 105, E 1 3. 6, ///T7, • PHI V, T66 , 

3 * PS I * ) 

WRITE (3,40) IPS I! IP), IP = 2tNPl) 

DO 120 1 = 2, M P 1 

12 0 WRITE (3, 60) PHI (I ) , (R(I , J) , J=2,NP1) 

OUTPUT R COORDINATES INSIDE SLOTS 

WRITE (3, 20)HEDR,0ATEl,DATE2fDATE3fCASE 
WRITE (3, 130) ITER, J TER, EPS, DEL 

1 3 C FORMAT* ///T2 , * R - COORDINATES INSIDE SLOTS', 

1146, 'ITER = • ,T53, I4,T61,'JTER = ' ♦ T68 , 1 4, T76 , • EPS ^ 
2 , T8 2 , E 1 3 • 6 , T99 , 'DEL = 1 , T 1 0 5 , E 1 3 . 6 , / / / T7 , • PS I • , T66 , 
3' PHI • ) 

WRITE (3, 40) (PH I (IP) ,IP = 2,MP1) 

IF (FLAG! • GT. 0 .OR. FL AG3 .GT. 0) GO TO 150 
WRITE(3,140) 

140 FORMAT! //T2, 'LOWER SLOT') 

WRITE (3, 60) PSI (NSTLPI) , ( RSLOTL ( I P ) , I P = 2 , MP 1 ) 

1 5 C WRITE(3,160) 

1 6 C FORMAT* //T2 , 'UPPER SLOT') 

WRITE (3, 60) PS l (NSTUP1 ) , (RSLOTU( IP), I P=2 , MP1 ) 

IF (NUT • EQ • 1) RETURN 

OUTPUT OF X - COORDINATES 
WRITE ( 3,20) HE0R,DATE1 ,DATE2 ,DATE3,CASE 
WRITE (3, 170) ITER, J TER, EPS, DEL 
17 0 FORK AT ( ///T2 , ' X - COORDINATES ( UN TR AN SFORMED )', 

1 T 46 , • ITER = ' , T53, I4,T61 ,' JTER = • , T68 , 14 , T76, • EPS = 
2,T8 2,E13.6, T99, 'DEL = • , T105 , E 1 3. 6 • ///T7 , 1 PHI • ,T66, 
3' PSI ' ) 

WRITE! 3, 40) (PSI! IP) ,IP = 2,NP1) 

00 180 1=2, MP1 

16C WRITE (3, 60) PHI (I),(X( I , J ) , J =2 , N PI ) 

OUTPUT X VALUES INSIDE SLOTS 

WRITE! 3 ,20) HEDR, DATE1 , DATE2 , DAT E3 , CASE 
WRITE (3, 190) ITER, JTER, EPS, DEL 
19 C F0RMAT(///T2,' X - COORDINATES INSIDE SLOTS', T46, 

1 ' ITER = 1 ,T53, I4,T61, ' JTER = • , T68 , 14 , T76 , • EPS = ', 
2T82, E13.6,T99, 'DEL = T105, E 13.6, ///T7 ,• PSI ', T66, 

3 ' PHI ' ) 

WRITE! 3,40) ( PH KIP), I P= 2, MP1) 

IF ( FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) GO TO 200 
* WR I T E ! 3 , 140) 

WRITE (3,60) PS I (NSTLPI ), (XSLOTL(IP), IP = 2,MP1) 

200 WRITE(3, 160) 

WRITE* 3, 60) PSI (NSTUP1 ) , (XSLOTU! IP), IP=2,MP1) 

OUTPUT VELOCITY VECTOR ANGLES ( ALPHA ) 

WRITE (3, 20) HE0R,DATE1 , D A TE2 , DAT E3 , C A SE 
WRITE(3,210) ITER, JTER, EPS, DEL 
210 FORMAT ! / / /T 1 , ' OFLOW ANGLES ( ALPHA )', 146, 'ITER 

1T53, I4,T61, • JTER = ' , T6 3 , 1 4 , T?6 , ' E P S = , ,T82,E13.6, 
2T99 , 'DEL = • , T 105 , El 3 . 6 , ///T7 , ' PHI • ,T66 , • PSI • ) 
WRITE(3,40) (PSI! IP),IP = 2,NP1) 

DO 220 I=2,MP1 

220 WRITE (3, 60) PHI ! I ),( ALPHA! I ,J) , J=2,NP1 ) 
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OUTPUT FLOW ANGLES INSIDE SLGTS 

WRITE! 3,20) HEDR , DAT El , DATE2 , DAT E3 , CASE 
WRITE! 3,230) ITER, JTER , EPS, DEL 
23 C FORMAT ( ///T2 , ' FLOW ANGLES INSIDE SLOTS' » T46 , • ITER = 
1T53, I4.T61, ' JTER = ' , T68 » I A ,T76 , ' EPS = • , T32 , E 13-.6 , 
2T99 , 'DEL = ' , T 105, E 13. 6, ///T7, • PSI • ,T66, • PHI • I 
WRITE (3, 40) (PH II IP) , IP=2,MP1I 

IF (FLAGI .GT. 0 .OR. FLAG3 .GT. 0) GO TO 240 
WRITE(3,I40» 

WRITE(3,60)PSI (NSTLP1) , (ASLOTL! IP) , 1P=2,MP1) 

24C WRITE(3,I60) 

WRITE! 3 ,60) PSI (NSTUP1 ) , ( ASLOTU! IP), IP=2,MP1 ) 

OUTPUT X AND R - COORDINATES ( TRANSFORMED) 

IF RMULT > 0 

IF (RMULT .LE. 0.0) GO TO 2000 
YMl N=0.999E6 
DO 166 1=2, MP1 

IF(R(I,NP1) ,LT.YMIN)YMIN=R( I,NP1) 

166 CONTINUE 

RMUL T=RMULT/Y.M IN 
DO 250 J =2 , NP1 
DO 250 1=2, MP1 
R ( I , J) = RMULT * R ( I , J ) 

250 X(I,J) = RMULT * X(I,J) 

DO 260 1=2, MP1 

XSLOTL(I) = RMULT * XSLOTL(I) 

XSLOTU(I) = RMULT * XSLCTU(I) 

RSLOTL ( I ) = RMULT * RSLOTLU) 

26 C RSLCTU(I) = RMULT * RSLOTU(I) 

REWIND 5 
REWIND 12 

WRITEI5) ( !R( I , J ) , 1=1, MP2) ,J=1,NP2) 

WRITE (12) ( (X! I , J) , I = 1,MP2) ,J = 1 ,NP2) 
WRITE(3,20)HEDR,DATE1,DATE2,DATE3,CASE . 
WRITE(3,270)ITER, JTER, EPS, DEL 
270 FORMAT, !///Tl, • OR - COORDINATES ( TRANSFORMED )',T4o, 
1'ITER = * ,T53, I4.T61, ' JTER = ' , T68 , I 4 , T 76 , ' EPS = *, 
2T8?,E13.6,T99, 'DEL = • , T 105 , El 3 . 6 , / //T7 , • PHI • , T66 , 

3 ' PS I ' ) 

WRITE(3,40) ! PS l( IP ) , IP = 2,NP1 ) 

DO 280 1=2, MP1 

28 0 WRITE (3, 60) PHI (I ) , <R< I , J ) , J=2,NP1) 

WRITE! 3,20) HEDR ,DATE1 , DATE2 ,DATE3,CASE 
•WRITE (3, 290 ) ITER, JTER, EPS, DEL 
290 FORMAT! ///T2, ' R - COORDINATES IN SLOT ( TRANSFORMED )• 
1 , T46 , ' ITER = ' ,T53, I4,T6I, • JTER = ' , T68 , 14 , T76 , 

2'EPS = ' ,T32,E13.6,T99, 'DEL = • , T 105 , E 1 3 . fc , / //T7 , 

3 ' PSI ' , T66 , • PHI ' ) 

WRITE (3,40) (PHI! IP >, IP=2,MP1) 

IF (FLAGI .GT. 0 .OR. FLAG3 .GT. 0) GO TO 300 
W R I T E ( 3 , 140) 

WRITE (3,60) PSI (NSTLP1 ) , (RSLOTL < IP) , I P = 2 , M P 1 ) 

300 WR1TE13, 160) 

WP. ITE(3»60)PSI (NSTUP1 ) , <RSLOTU( IP ) , I P=2 » MPi ) 

WRI TE( 3, 20) HEDR, DA TE1 , DATE2,DATE3,CASE 
WRITE (3,310 ) ITER, JTER , EPS, DEL 
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31 C FORMAT* ///Tl , • OX - COORDINATES ( TRANSFORMED )',T46, 

1'ITER = » ,T53, I4,T6l, • J7ER = ', T68 » I 4 , T76 ,' EPS = ', 
2T82,E13.6,T99, 'DEL = • , T 105 ,E 13 . 6, / //T7, • PHI • , T66 , 

3 • PS I ' ) 

WRITE(3,40) (PSI( IP > , IP=2 ,NP1 ) 

DO 320 1=2, MP1 

32 C WRITE (3, 60 1 PHI (I) , (X( I , J ) , J = 2 , NP 1 ) 

WRITE (3,20) HEDR »DATE1 , DATE2 ,DAT E3 ,C A SE 
WRITE(3,330> ITER, JTER, EPS, DEL 
330 FORMAT ( / / /T 2 * ' X - COORDINATES IN SLOT ( TRANSFORMED )• 

1 » T46 , • ITER = ' ,T53, 14.T61, ' JTER = ' , T68 , 1 A , T76 , 

2 ' EPS = • ,T82,E13.6,T99, 'DEL = • , T 105 , E13 . 6 , ///T7 , 

3'PSI • ,T66, • PHI • ) 

WRITE (3, 40) <PHI(IP),IP=2,MP1) 

IF ( FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) GO TO 340 
WRITE(3»140) 

WRITE(3,fe0)PSI (NSTLP1) , (XSLOTL( IP) , IP=2,MP1) 

340 WRITE (3, 160) 

WRI TE(3,60)PSI (NSTUPl ) , ( XSLOTU( IP), IP=2,MP1) 

SUMMARY OF WALL STREAMLINE COORDINATES 

2000 I = MSTAGU ♦ 1 
J = NSTAGU +1 
Q(MSTUP1 , NSTUPl) = 0.0 
0( MSTLP1 ,NSTLP1 ) = 0.0 
QSL0TUIMSTUP1) = 0.0 
QSLOTL(MSTLPl) = 0.0 

SS = SORT ( (X( I , J ) -X C 1*1 ,J ) I’M X ( I , J ) -X I U1,J ) ) 

1 + (R< I ,J)-R< 1*1, J) )*( R( I , J )-R( 1+1 , J) ) ) 

SS = SS - SKI+1) 

SSS = SQRTI (XSL0TU(I)-XSL0TU(I+1) )*(XSLOTU(I)-XSLOTU( 1 + 1) ) 

1 ♦ (RSLOTU(I) -RSLOTUI H-l ) )*<RSLOTU( I ) -RSLOTU ( 1 + 1 ) ) 

SSS = SSS - S2 ( 1 + 1 ) 

MSTUP2 = MSTAGU + 2 
DO 345 I =MSTUP2 ,MP1 
SKI) = SKI) ♦ SS 
345 S2( I ) = S2( I ) + SSS 

WRI TE( 3, 20) HEDR, DATEl, DATE2.DAT E3 .CASE 
DO 350 1=2, MP1 
DO 350 J=2 , NP1 
P = Q ( I , J ) 

350 Q(I,J) = DEXP( P ) 

WRITE! 3, 360) 

360 F0RMAT(///T1,' 0 SUMMARY OF INLET WALL COORDINATES OF • , 

1T39, ‘DIFFUSER (UPPER WALL)* ,///T8, * X * , T27, * R * ,T47, ‘Q* , 

2T68, ‘S', T85, 'ALPHA'/) 

ICOUNT = 1 

DO 370 1=2, MSLUP1 

ICOUNT = ICOUNT 4- 1 

IF (ICOUNT .LT. 23) GO TO 370 

ICOUNT = 1 

WRI TE(3,20) HEDR.DATEl , DATE2 , DAT E3 , C ASE 
WRITE(3, 360 ) 

370 WRITE(3,375)X( I ,NP1 I »R( I ,NPl ) , 0 ( I , NP1 ) , S ( l I , 

1ALPHA ( I , NP1 ) 

375 FORMAT (/1X,E13.6,4(7X,E13.6) ) 

WRI TE( 3, 20) HEDR, DATE 1 ,DATE2,DATE3, CASE 
WRITE(3,380) 

38C FORMAT! ///T2, ' SUMMARY OF DOWNSTREAM WALL COORDINATES', 
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IT AO * * OF DIFFUSER (UPPER WALL) ’ , / // T8 , • X' , T27 , • R ' , TA7 , 

2*0* , T67 » ' S' * T9 5 , ' ALPHA* /) 

ICOUNT = 1 

DO 390 I = MSTUP I » MP1 

ICOUNT = ICOUNT ♦ 1 

IF (ICOUNT .LT. 23) GO TO 390 

ICOUNT = 1 

WRITE ( 3,20) HEDR,DATE1,DATE2,DATE3,CASE 
WRI T E ( 3 , 3 80 ) 

39C WR I TE ( 3 » 375 ) X( I,NSTUP1) , R < I ,NSTUPl ) ,Q( I.NSTUPl), 

1SKI ), ALPHA! I.NSTUPI ) 

WR I TE( 3, 20 )HEDR,DATE1 , DATE2.DAT E 3, CASE 
WR I T E ( 3 , AOO ) 

AO 0 F0RMAT(///T2,' SUMMARY OF WALL COORDINATES INSIDE OF ', 

1TA0, 'DIFFUSER SLOT (UPPER W AL\ ) • , / / /T 8 , * X* , T27 , • R * , 

2TA7,'Q* ,T67 ,' S ' ,T8 5, • ALPHA*/) 

ICOUNT = 1 

DO A 1 0 I = MSTUPI,MSLUP1 
ICOUNT = ICOUNT + 1 
IF (ICOUNT .LT. 23) GO TO AIO 
ICOUNT = 1 

WRITE(3,20)HEDR,DATEI,DATE2,DATE3,CASE 
WRITE! 3 »A 00 ) 

AIO WRITE(3*375(XSLOTU( I ) *PSLOTU(I) » OSLGTU ( I ) * S2 ( I > * 

1 ASL0TU( I ) 

IF (FLAG1 .GT. 0 .OR. FLAG3 . GT . 0) RETURN 
I = MSTAGL + 1 

J = NSTAGL + 1 

SS = SORT ( ( X( I ,J ) -X ( I + l, J ))*( X( I , J )-X( I-fl , J ) ) 

1 + (R( l ,J )-R( 1+1, J ) )*( R (I , J )-R( Kl, J ) ) ) 

SS - SS - SBAR 1 ( I 1 ) 

SSS = SORT ( (XSLOTL( I )-XSLOTL( I + 1 ) )#(XSLOTL( I )-XSLOTL( I ♦I) ) 

1 ♦ (RSLOTLf I)-RSLOTL( 1+1) ) *( RSLOTL ( II-RSLOT L ( I *1 ) ) ) 

SSS = SSS - SB AR2 ( I + I ) 

MSTLP2 = MSTAGL + 2 
DO A 15 I=MSTLP2,MPL 
SBAR 1 ( I ) = SHARK I ) + SS 
A15 SBAR 2 ( I ) = S3AR2 ( I ) + SSS 

WRITE! 3,20) HEDR.D ATE1.DAT E 2, DATE3, CASE 
WR I TE ( 3 » A20 ) 

A20 FORMAT! ///T2, ' SUMMARY OF INLET WALL COORDINATES OF *, 

1T39, 'DIFFUSER (LOWER W ALL )*,/// T8 ,* X *, T27 ,* R *, TA7 ,* Q • , 

2T67, 'S' , T85 , 'ALPHA* /) 

ICOUNT = 1 

DO A 30 1=2, MSLLPi 

ICOUNT = ICOUNT ♦ 1 

IF (ICOUNT .LT. 23) GO TO A30 

ICOUNT = 1 

WRI TE( 3, 20) HEOR,DAT El , DATE2 , DATE3 , CASE 
WRITE! 3, A20) 

A3 0 WRITE (3, 3 75 ) X ( I ,2) , R< I , 2) ,Q( I , 2) ,SBAR( I) , ALPHA! I ,2) 

WRITE (3, 20) HEDR.DATEl ,DATE2 ,DATE3 ,CASE 
WRI TE ( 3 , AAO ) 

A A 0 FORMAT! ///T2 ,' SUMMARY OF DOWNSTREAM WALL COORDINATES', 

1TA0,' OF DIFFUSER (LOWER W ALL )',/// T 3 ,' X ', T27 ,' R ', TA7 ,' Q • 

2 , T6 7 , ' S ' , T85, • ALPHA'/ ) 

ICOUNT = 1 

DC A 50 I = MSTLP1,MP1 

ICOUNT = ICOUNT + 1 

IF (ICOUNT .LT. 23) GO TO A50 
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ICOUNT = 1 

WRITE (3, 20 >HEDR,DAT El, DATE2.DAT E3,CASE 
WRI TE ( 3 » 440 ) 

45 C WRITE(3,375)X(I,NSTLP1),R(I ,NSTLP1 ) »Q( I ,NSTLP1 ) , 

1SBAR1 { I ) , ALPHA ( I, NSTLPl ) 

WRI TE( 3,20) HEDR, DATE 1 ,DATE2 ,D4TE3, CASE 
WR I T E I 3 , 460 ) 

460 FORMAT ( ///T2 * ' SUMMARY OF WALL COORDINATES INSIDE OF ', 
1T40, 'DIFFUSER SLOT (LOWER W ALL )',/// T 8 ,' X ', T27 ,' R • , 
2T47, »T67*'S',T85» 'ALPHA' /) 

ICOUNT = 1 

DO 47 0 I-MSTLP1.MSLLP1 

ICOUNT = ICOUNT + 1 

IF (ICOUNT .LT. 23) GO TO 470 

ICOUNT =1 . 

WRI TE< 3,20) HEDR.DATEl , DATE2 , DAT E3 , CASE 
WR I T.E ( 3 , 460 ) 

47 0 WRITE (3, 375 >XSLOTL( I ) ,R SLOTH I ) ,QSLOTL( I ) » $8 AR2 ( I ) , 
1ASL0TH I ) 

RETURN 

END 

C «=**#*#**** SUBROUTINE GRAPIC ********** 

C 

C THIS SUBROUTINE PROVIDES GRAPHICAL OUTPUT FOR MAIN 

C 

SUBROUTINE GRAPIC 
DIMENSION BUFR(500) 

DIMENSION DX(150) , ES IMP < 150 > , FS IMP < 150 ) ,GSIMP( 150) , 

1 GRAND ( 150) , GIFT ( 150 ), PHI ( 150) ,PSI ( 150) , 

2 DPSI( 150 ) ,DPHI ( 150) ,S( 150) , SBAR ( 150) , 

3 SI ( 150) , SBAR 1 ( 150) , S2( 150) ,SBAR2( 150) , 

4 HEDR< 15) ,EPHI ( 150 ) ,FPHI ( 150) ,GPHI ( 1501 , 

5 EPSI( 150) ,ASL0TL(150) ,XSLOTU( 150) ,XSLOTL( 150) 

6 ,FPSI ( 150) ,GPSI (1501 ,RSL0TU(150) ,RSLOTL( 150) , 

7 0SL0TU( 150 ) , OSLOT L ( 1 50 ) ,ASL0TU(150) 

DIMENSION 0(150,50) ,R( 150, 50), X(150, 50) 

DIMENSION XX ( 150) ,XY( 150) , L ABL ( 20 ) , RABL ( 20 ) , WABL ( 20 ) 
DOUBLE PRECISION P 

INTEGER FLAG1 , FLAG2 , FL A G3 , FL AG4 , C AS E , DAT E 1 , DAT E 2 , DATE 3 
COMMON / C2 / FLAG1 . FL AG2. FLAG3. FL AG4 . HEDR , UATE1 . 0ATE2 

1 
2 

3 

4 

5 

6 
7 

REWIND 4 
REWIND 5 
REWIND 12 
MP1 = M 1 
NP1 = N + 1 
MP2 = M + 2 
NP2 ■= N + 2 
MM 1 = M - 1 
NM1 = N - 1 
MSTUP2 = MSTAGU + 2 
. MSTUP1 = MSTAGU + 1 
MSTUM1 = MSTAGU - 1 
NSTUP1 = NSTAGU + 1 
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NSTUM1 = NSTACU - 1 


GRAP 

40 


MSTLP2 = MS TAG L ♦ 2 


GRAP 

41 


MSTLP1 = MSTAGL ♦ 1 


GRAP 

42 


M STL Ml = MSTAGL - 1 


GRAP 

43 


NSTLP1 = NSTAGL + 1 


GRAP 

44 


NSTLM1 - NSTAGL - 1 


GRAP 

45 


MSLUP1 = MS LOT U ♦ 1 


GRAP 

46 


MSLUM1 = MS LOT U - 1 


GRAP 

47 


NSLUPl = NSLOTU + 1 


GRAP 

46 


. NSLUM1 = NS LOT U - 1 


GRAP 

49 


MSLLP1 = MSLQTL ♦ 1 


GRAP 

50 


MSLLMl = MS LOT L - 1 


GRAP 

51 


.NSLLP1 = NS LOT L ♦ 1 


GRAP 

52 


NSLLMl = NSLOTL - 1 


GRAP 

53 


READ 14) ( (0( I «J ) »1=1»MP2) ,J=1»NP2) 


GRAP 

54 


READ <5M <R( I ,J > , I=1,MP2 ) ,J=1 ,NP2) 


GRAP 

55 


READ (12) l * X ( I, J ) , 1= 1,MP2 ) , J = 1,NP2) 


GRAP 

56 


CALL PLOTS! 8 UFR, 100) 


GRAP 

57 


DATA LABL ( 1 ) / ' , /,LA8L(2)/ i , /,LABL(3)/ > 

•/, 

GRAP 

58 


1 LA3L(4)/« '/ 


GRAP 

59 


-CAT A LABL ( 5 I / 1 X A • / , L A 8 L ( 6 ) / 1 X 1 S •/ 


GRAP 

60 


DATA LABL ( 7 ) / ' ' / , LABL ( 8 ) / ' '/,LABL(9J/' 

V, 

; GRAP 

61 


1 LABL (101/* •/ 


GRAP 

62 


DATA LABL ( 1 1 ) / * 1 / , LABL ( 1 2 ) / * */, 


GRAP 

63 


1 LABL ( 1 3 ) / • V ,LA8L( 14)/ 1 •/ 


GRAP 

64 


DATA LABL (15)/* Y-K 1 / , L ABL ( 1 6 ) / • AX I S 1 / 


GRAP 

65 


DATA L ABL ( 1 7 > / ' ' / , L ABL 1 18 ) / • •/, 


GRAP 

66 


l LABL ( 1 9 ) / 1 1 / , L ABL ( 20 ) / f •/ 


GRAP 

67 

c 



GRAP 

68 

c 

WRITE HEOR 


GRAP 

69 

c 



GRAP 

70 


, CALL SYMBOL ( 0» 5 tO • 0 »0 • 14 9 HEOR f 90 • 0 v 60 ) 


GRAP 

71 

c 



GRAP 

72 

c 

GRAPH 2-D AND BELL DIFFUSERS 


‘ GRAP 

73 

c 



GRAP 

74 


IF ( FLAG1 « LE • 0 .AND. FLAG3 .LE. 0 ) GO TO 1000 


GRAP 

75 


DC 10 J = 1 1 N 


GRAP 

76 


10 XX( J ) = R ( M P 1 , J + 1 ) 


GRAP 

77 


CALL ADJST( XX, N, 1 ,9 .0 ,0.0,SCAL, XMIN, XY , 1 ) 


GRAP 

78 


DO 20 1=2, MP1 


GRAP 

79 


DO 20 J = 2 , N P 1 


GRAP 

80 


Xil , J) = XI I , J ) / SCAL 


GRAP 

81 


20 R (I , J ) = R( I , J ) / SCAL 


GRAP 

82 


DC 30 I = MSTUP1 tMSLlJPl 


GRAP 

83 


XSLOTU(I) = XSLOTU(I) / SCAL 


GRAP 

34 


3C RSLGTU(i) = RSLOTU(I) / SCAL 


GRAP 

85 

c 



GRAP 

86 

c 

DRAW AXES FOP DIFFUSER PLOT 


GRAP 

87 

c. 



GRAP 

88 


- CALL COAXES ( 3.0,0. 5, X( MPl f 2 ) ,9.0, SCAL, SCAL, X( 2, 2 ) , 


GRAP 

39 


1 R ( 2 9 2 ) , L ABL , 0 ) 


GRAP 

90 


DO AO 1 = 1, M 


GRAP 

91 


DC AO J = 2 ♦ N P 1 


GRAP 

92 


X( I , J) = X( 1 + 1 ,J) 


GRAP 

93 


AG RU , J ) = P( 1+1 ,JI 


GRAP 

94 


DC 50 1=1, M 


GRAP 

95 


DC 5 0 J = 1 , N 


■ GRAP 

96 


X(I , J) = X(I,Jfl) 


GRAP 

97 


5 C R ( I 9 J ) = R < I , J + 1 ) 


GRAP 

98 

c 



GRAP 

99 

— 
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GRAPH STREAM LINES 

DO 70 J -1 f N STAGU 
DO 60 1=1 ,M 
XX( I ) = X< I , J) + 3.0 
6 C XY( I ) = R ( I t J) + 0.5 
70 CALL LINE<XX,XY,M, 1,0,0) 

DO 90 J = NSTUP1 * N 
00 80 1= 1 ♦ M SLQTU 
XX( I ) = X( I , J) + 3.0 
80 XY( I ) = R( I , J) ♦ 0.5 
90 CALL LINE ( XX , X Y , MSL OTU * 1,0,0) 

IOT = MSLUP1 - MSTAGU 
DQ 100 1=1, IOT 

XX(I> = XSLOTLU I+MSTAGU) + 3.0 
IOC XYfl) = RSL0TU( I+MSTAGU) + 0,5 
CALL LINE(XX,XY, IOT, 1,0,0) 

graph potential lines 

DO 120 1=1, MSTAGU 
DC 110 J= 1 , N 
XX( J ) = X< I , J) + 3.0 
lie X Y ( J ) = R ( I , J ) > 0.5 
120 CALL LINE(XX,XY,N*1,0,0) 

DO 190 1= MSTU P 1 , M 
DO 130 J= 1 , NST AGU 
XX( J ) = X ( l ,J) + 3.0 

130 XY ( J ) = RU , J) + 0.5 
14C CALL LI NE ( XX * X Y * NST AGU * 1 ,0,0) 

JOT = N - NSTUM1 
DO 160 l=M$TUPl ,MSLGTU 
DO 150 J=2 , JOT 
XX(J) = X( I , J+NSTUMi) ♦ 3.0 
1 5 C X Y( J ) = R( I , J + NSTUM1 ) + 0.5 
XX ( 1 I = XSL OTU ( I +1 ) + 3.0 
XY C 1 I = RSLGTUU+ll + 0.5 
160 CALL LINE(XX*XY, JDT, 1,0,0) 

XL I N = X ( M , 1 ) 4-9.0 
GO TO 2000 

GRAPH ANNULAR DIFFUSER 

100C CONTINUE 

A SH 1 = R ( MP 1 , N STL P 1 ) 

ASH = A SHI 

ASH2 = RSLOTLI MSLLP1) 

IF <ASH2 .LT. ASH1) ASH = ASH2 
IF ( R ( 2 , 2 ) .LT. ASH ) ASH = R(2,2> 
DO 170 J = 2 , NP1 
DO 170 1=2, MP1 
170 R ( I * J ) = R ( I, J) - ASH 
DO 173 I=MSTUP1,MSLUP1 
173 RSLOTU(I) = RSLOTUU) - ASH 
DO 176 I = MS TLP 1 , MSLLP1 
176 RSLCTLCI) = RSLOTL(I) - ASH 
JCT = NST AGU - NST AGL + 1 
TRA SHI = R ( MP1 , NSTUP1 ) 

. TRASH = TRA SHI 
TRASH2 = RSLOTU ( MSLU PI ) 
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IF (TRASH2 .GT. TRASH1) TRASH = TRASH2 

IF ( R ( 2 , NP 1 ) .GT. TRASH I TRASH = R(2,NP1I 

00 180 J= 1 , JOT 

180 XX(J) = R(MP1, J+NSTAGL ) ♦ ASH 

xxn ) = ash 

XX(JOT) = TRASH 

CALL ADJST( XX, JOT, 1 , 9 . 0 , 0. 5 , SC AL , XM I N , XY , 1 ) 

DO 190 1=2, MP1 
00 190 J=2 , NP 1 
X ( I , J) = X( I,J) / SCAL 
190 R ( I , J ) = R(ItJ) / SCAL + ASH - XMIN 

CALL COAXES (3.0,0. 5 ,X { MP1 , NSTLP1 ) ,9 .0 , SCAL , SCAL , X ( 2 , 2 ) , XMIN, 
1 LAB L , 0 ) 

00 200 I=MSTLP1,MSLLP1 
XSLOTL(I) = XSLOTL(I) / SCAL 
200 RSLOTUI) = RSLOTL(I) / SCAL + ASH - XMIN 
DO 210 I=MSTUP1 ,MSLUP1 
XSLOTUII) = XS LOT U ( I ) / SCAL 
210 RSLOTU(I) = RSLOTUI II / SCAL ♦ ASH - XMIN 
00 220 1=1, M 
DO 220 J=2 , NPl 
X( I , J) = X( 14-1 , J ) 

220 R(I,J) = R ( I +1 ,J) 

DO 230 1=1, M 
00 2 30 J = 1 , N 
X(I,J) = X( I,J+1) 

23C R( I , J) = R ( I , J + 1 ) 

GRAPH STREAMLINES 

00 250 J=1 , NSTLM1 
DO 240 1 = 1, MSL OTL 
X X C I I = X( I , J) ♦ 3.0 
24C XV(I) = R ( I ,J) +0.5 
250 CALL LINE (XX,XY,MSLOTL , 1,0,0) 

00 270 J = NSTAGL , NST AGU 
DO 260 1=1, MP1 
XXII) = X( I , J) + 3.0 
26C XY(I) = R < I , J ) + 0.5 
270 CALL LINE(XX,XY,M, 1,0,0) 

DO 290 J=NSTUP 1 , N 
DO 2 BO 1 = 1, MSL OTU 
XX( I ) = X ( I ,J) + 3.0 
28 C XYU1 = R ( I ,J) ♦ 0.5 
29 C CALL LINE! XX,XY,MSLOTU, 1,0,0) 

I GT = MSLOTL - MSTAGL + 1 
DO 300 1=1, IOT 

XXII) = XSLOTLI I+MSTAGL) + 3.0 
300 XY(I) = RSLOTL < I+MSTAGL ) ♦ 0.5 
CALL LINE (XX, XY, IOT, 1 ,0,0) 

IOT = MSLOTU - MSTAGU 4- 1 
DO 310 1=1, IOT 

XX(I) = XSLOTU( I+MSTAGU) + 3.0 
310 X Y ( I ) = RSLOTUC I+MSTAGU) + 0.5 
CALL LINE(XX,XY,IOT ,1 ,0,0) 

GRAPH POTENTIAL LINES 

MINN = MSTAGL 
MAXX = MSTAGU 
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IF (MSTAGU .LT. MSTAGL) MAXX = MSTAGL 
IF (MSTAGU .LT. MSTAGL) MINN = MSTAGU 
DO 330 1=1, MINN 
DO 320 J= 1 , N 
XX( J ) = X( I , J) + 3.0 
32 C X Y ( J ) = R ( I ,J) + 0.5 
33 C CALL L1N£(XX,XY,N, 1,0,0) 

MINN PI = MINN + 1 

IF (MSTAGL .EQ. MSTAGU) GO TO 415 
DO 410 I=MINNP1* MAXX 
IF (MSTAGL .GT. MSTAGU) GO TO 380 
DO 360 J= 1, NSTLM1 
XX( J ) = X ( I , J ) + 3.0 

36 C X Y( J ) = K ( I , J ) + 0.5 

XX(NSTAGL) = X SLOT L ( I ) + 3.0 
XY(NSTAGL) = RSLOTL(I) + 0.5 
CALL UNE(XX,XY, NSTAGL, 1,0,0) 

JOT = N - NSTAGL + 1 
DO 370 J= 1 , JOT 
XX(J) = X(I , J+NSTLM1) + 3.0 
37 C X Y ( J ) = R( I , J+NSTLM1 ) + 0.5 
CALL LINE(XX,XY, JOT, 1,0,0) 

GO TO 410 

360 DO 390 J= 1 , NST AGU 

XX ( J ) = X ( I ,J) + 3.0 

390 X Y ( J ) = R ( I ,J) + 0.5 

CALL LINE(XX,XY,NSTAGU, 1,0,0) 

JOT = N - NSTAGU + 1 
DO 400 J=2 * JOT 
XX(J) = X ( I , J+NSTUM 1 ) + 3.0 

400 XY ( J ) * R{ I , J+NSTUM1 I + 0.5 
XX(1) = XSLOTU(I) + 3.0 
X Y ( 1 ) = RSLOTU ( I ) +0.5 
CALL LINE(XX,XY, JOT, 1,0,0) 

4 1 C CONTINUE 
415 CONTINUE 

JOT = NSTAGU - NSTAGL + 1 
DO 430 I = MAXX» M 
DO 420 J=1 , JOT 
XX(J) - X ( I , NS TLM 1 + J ) + 3.0 

420 X Y( J ) = R ( I , NS TLM 1 + J ) + 0.5 
430 CALL LINE( XX,XY,J0T,1 ,0,0) 

JCT = NSTAGL - 1 
DO 450 I=MSTLP1*MSLQTL 
DO 440 J= 1 , JOT 
XX(J) •= X(l,J) + 3.0 
440 X Y( J ) = R( I , J ) + 0.5 

XX(NSTAGL) = X SLOT L ( I +1 ) + 3.0 
XY(NSTAGL) = R SLOTL ( I + 1 ) + 0.5 
450 CALL LINE(XX,XY, NSTAGL, 1,0,0) 

JOT = N - NSTAGU + 1 
DO 470 I=MSTUP1,MSL JTU 
DO 460 J=2 , JOT 
XXIJ) = X( I , J+NSTUM1) + 3.0 
460 X Y( J ) = R ( I , J+NSTUM1 ) + 0.5 
XX( 1 ) = XSLOTIM 1 + 1) +3.0 
X Y ( 1 ) = RSLOTU ( 1 + 1 ) + 0.5 
470 CALL LI NE (XX, XY, JOT, 1,0*0) 

XL I N - X(M, NSTAGL ) + 9.0 

C 
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GRAPH ANNULAR DIFFUSER IN PERSPECTIVE 

DC 48C J = 1 , N 
DO 480 I = 1 f M 

R ( I , J) = SCAL * R ( I » J ) + ASH 
460 X< I , J) = SCAL * XI I ,J ) 

JOT = NSTAGU - NSTAGL + 1 
DO 490 J= 1 , JOT 
490 XX(J) = R (M tNS T AGL + J ) 

X X ( 1 ) = 0.0 

XX(JOT) = TRASH + ASH 

CALL AOJST( XX, JOT, 1 ,9.0,0.5,ZCAL,XMIN, XY,1 ) 

DO 500 I -= 1 , M 
DO 500 J = 1 , N 
X( If J) = X( IfJ) / ZCAL 
500 R ( I , J ) = R ( I , J ) / ZCAL 
DC 510 I=MSTLP1,MSLLP1 

XSLOTL(I) = ( XSLOTL(I) * SCAL ) / ZCAL 

510 RSLOTL(l) = ( RSLOTLU) * SCAL + ASH ) / ZCAL 

DO 520 I-MSTUPl.MSLUPl 

XSLOTUm = ( XSLOTUI I) * SCAL ) / ZCAL 

520 RSLGTU(I) = ( RSLOTU(I) * SCAL + ASH ) / ZCAL 

CALL CGAXES(XLIN,0.5 ,X(M, NSTAGL ) ,9.0, ZCAL, ZCAL, 0.0, 0.0, LABL,0) 

GRAPH STREAMLINES 

OG 530 I = 1 , MSL OTL 
XX( I > = X(I f 1 ) f XL IN 
530 XY( I J = R(I ,1) ♦ 0.5 

CALL LINE(XX,XY,MSLOTL, 1,0,0) 

IQT * M • MSTAGL +1 
00 540 1=1, IQT 

X X ( I I = X( I+MSTLM1, NSTAGU ♦ XLIN 
540 XY( I) - R( I+MSTLM1 , NSTAGL) ♦ 0.5 
CALL LINE(XX,XY,I0T,1 ,0,0) 

IOT = MSLOTL - MSTAGL +1 
DO 550 1=1, IOT 

XX(I) = XSLOTL ( I+MSTAGL) ♦ XLIN 
550 XY( I) = RSLOTL ( I+MSTAGL ) ♦ 0.5 
CALL LINE(XX,XY, IQT, 1,0,0) 

DO 560 1 = 1, MSL OTU 
XX ( I ) = X ( I , N ) + XLIN 

56 C XY( I) = R(ItN) ♦ 0.5 

CALL LINE! XX , X Y , MS LOT U ,1,0,01 
IOT = M - MSTAGU + 1 
DO 570 1=1, IOT 

XX(I) = X( H-MSTUM1 , NSTAGU) + XLIN 
570 X Y ( I ) = RM+MSTUM1, NSTAGU) ♦ 0.5 
CALL LINE(XX,XY, IOT, 1,0,0) 

IOT = MSLOTU - MSTAGU > 1 
DO 580 1=1, IQT 

XX < I ) = XSLGTU( I+MSIAGU) + XLIN 
580 X Y ( I ) = RSL0T1.M I+MSTAGU) ♦ 0.5 
CALL LINE (XX, XY, IQT,1 ,0,0) 

GRAPH POTENTIAL LINES 

DO 590 J = 1 , N 

XX ( J ) = X( 1 , J) + XL IN 
590 X Y ( J ) = R(1,J) ♦ 0.5 
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CALL LINECXX,XY f NtltO,0) 

JOT = NSTAGU - NSTAGL + 1 
DO 600 J=1 9 JOT 

XX(J) = X(M, J+NSTLM1) ♦ XLIN 
600 XY(J) = RfM, J>NSTLM1 ) + 0.5 
CALL LINE(XX,XY, J0T,1 ,0,0) 

DO 610 1=1, M 
DO 610 J=1 , N 

X(I,J) = ( X(I,J) * ZCAL ) / SCAL 
610 R ( I , J ) = ( R ( I , J ) * ZCAL - ASH ) / SCAL 
XLIN = XLIN + X(M, NSTAGL) ♦ 9.0 

GRAPH VELOCITY PROFILES FOR DIFFUSER WITH 
NEUMANN BOUNDARY CONDITIONS 

2000 IF ( FLAG2 .LE. 0) GO TO 3000 

DATA WABL ( l ) / f •/,WABL(2)/ f •/•WABL(3)/ # */, 

1 WABL ( 4 ) /* W 

DATA WABL(5I/* Q A • / , W ABU 6 ) / • XI S •/ 

DATA WABL(7 ) /• • / , WA8L ( 8 ) / 1 */,WABL(9)/ # •/, 

1 WABL ( 1 0 ) / • •/ 

DATA WA&LU1)/* • / , WABL < 12 ) / • •/, 

1 WABL (13)/* 1 / , WABL ( 14 ) / * •/ 

DATA WABL (151/ • N A • / , W ABL ( 16 ) / • XI S # / 

DATA WABL (17)/* a / , WABL ( 18 ) / 1 •/, 

1 WABL ( 1 9 ) / • , /,wABL(20»/» •/ 

DO 615 J=2, NP1 
DO 615 1=2, MPl 

p = q c i , j > 

615 Q ( I , J) = DEXPI P) 

XX(NPl) = 0.0 
DO 620 J=1 , N 
62 C XXCJ ) = Q(2,J+l) 

CALL AD J ST ( XX, NPltl ,9.0 , 0 . 5 , SSC AL , XM I N , XY , 1 ) 

CALL CGAXE$( XL IN, 0.5, 12*0,9.0, SSCAL , SC AL , XM I N , 0 . 0, W ABL , 0 ) 
DQ 630 J=1,N 

XX(J) = ( Q ( 2 , J+ 1 ) - XMIN ) / SSCAL ♦ XLIN 
630 XY( J ) = R(1 ,J) ♦ 0.5 

CALL LINE (XX,XY,N,1 ,0,0) 

IF ( FLAG1 .GT. 0 .OR. FLAG3 .GT. 0) NSTAGL = 1 
JOT * NSTAGU - NSTAGL ♦ 1 
DO 640 J= 1 , JOT 

XX ( J ) = ( Q(MP1 ,J + NSTAGL) - XMIN ) / SSCAL + XLIN 
64 C XY(J) = RIM, J*NSTAGl-l ) ♦ 0.5 
CALL LINE(XX,XY, JOT, 1,0,0) 

MOVE PEN TO FINISHED POSITION 

3000 XSET = XLIN * 15.0 

CALL PL0T1 (XSET ,0.0, -3) 

WRITE(3,650) 

65 C F0RMAT(1H1,T2, •GRAPHICAL OUTPUT COMPLETED PROGRAM •, 

1T37, • SIGNING OFF* ) 

RETURN 

END 
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APPENDIX B - DIFFUSER DESIGN 


This appendix describes the procedure used in designing the 18-inch 
(46 cm) annular diffuser. One should realize that the final design contour 
is rarely if ever the result of a first trial program. Several reruns are 
nearly always necessary to achieve the desired geometry because such parameters 
as AR and L/W are not input directly even though specific values of these may 
be the ultimate goal. 

For the 18-inch (46 cm) annular diffuser, the aim was to design a diffuser 
with the following features: (l) an area ratio of about 3-1 (2) an inlet 

width of about 2 inches (5-08 cm) (3) an exit outer diameter of about 1.8-inches 
(46 cm) (4) an effective length for boundary layer growth upstream of the 
slot of about 3 inches (7.6 cm). I terns (l) and (4) were the same as the pre- 
viously designed 9-inch (23 cm) diffuser, while items (2) and (3) were 
twice as large. 

The design program was used to design both the accelerating nozzle 
(inlet section) and the d i f fuser i tsel f as one flow channel. Wall velocity 
vector magnitude ratio (Q.) increasing from 1.0 to 3-0 in the nozzle portion 
and from 3.0 to 3.5 from the narrowest portion of the channel to the slot were 
prescribed. By prescribing a favorable velocity gradient upstream of the slot, 
it was thought that slot widths ranging from the design width to several times 
the design width could be tested. (As the slot is widened, there is a 
tendency for the wall velocity just upstream of the slot to decrease, thus 
producing a possibly unfavorable velocity gradient.) Downstream of the slot, 

Q. along the wall was prescribed to.be constant at 0.-83. The Q. ratio of 3.0 
at diffuser inlet to 0.83 at exit is larger than the desired AR because suction 
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removes some of the flow and because the diffuser inlet and exit profiles 


were not precisely uniform. The value 0.83 was arrived at after several trial 

0 - FS) d. 

runs. (A good first trial for Q is 7-5 which is 0.95 in this case.) 

AK 

To get the appropriate values of <f> for the desired velocities, recognize that 

f 2 

<t >2 “ <f>| “ J Q dS 

s i 

After adopting a model such as Figure 21 corresponding values of Q and <j> can 
be found. 



Figure 21. Prescribed wall velocity as design input. 


The proper range of ^ can be determined by applying 


2 



0 f cos a j , 
= 2 J Q 
<)> * 


at either the inlet or the exit passage. At inlet, a 7-inch (17.8 cm) outer 
radius and a 5“ i nch (12.7 cm) inner radius were desirable since these dimensions 
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would enable the first three items of the previously listed characteristics to 
be attainable. Thus 

y a 3(7 2 - 5 2 ) = 36 

was a good first trial value. (For the final design, V was adjusted to 33-5-) 
Increments At were selected to be much smaller at each wall, and increments 
A<f> were selected smaller in the vicinity of the slot. This enabled better 
resolution to be obtained in the troublesome slot region and a minimizing of 
the error introduced by the wedge flow patch at the stagnation point. The pro- 
gram requires that A<j> = Aty in the grid network immediately around the stagnation 
point, and these were arbitrarily given a value of 0.05 for this program. The 
upper slot contained 7 streamlines (including wall and stagnation point 
streamlines) and 4 potential flow lines. The lower slot contained 6 streamlines 
and 4 potential flow lines. The number of streamlines contained in each 
slot and the exact value of each Aip in the vicinity of the slot must be 
adjusted so that the flow area represented by these AiJj's is sufficient to pass 
the design suction flow for that slot. If, for example, 2% suction is desired 
through the inner slot and 4% through the outer slot (these values were close 
to but not precisely the design suction rates), then 

Z A ip = .02 E Aip = .02V 

across inner slot total 

and Z A ijj = . 04 Z A <J; = .04v 

across outer slot total 

For the previously named first trial V of 36, this would yield 0.72 as 
the value of i}> for the lower slot stagnation streamline and 36.0 - 1.44 or 
34.56 for the upper slot stagnation streamline. The final design had O .65 as 
the lower slot stagnation streamline, 32.4 for the upper slot stagnation stream- 
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line, and 33-5 for ¥. This is equivalent to 0.65/33-5 or 1.94% lower slot 

suction and (33-5 “ 32.4)/33*5 or 3*28% upper slot suction and 5-22% total 
suction. The design suction percentage was determined by assuming the boundary 

layer growth to be the same as that along a flat plate, and the velocity pro- 
file to obey the 1/7 power law, and applying Taylor’s criterion to determine 
how much fluid must be removed. 

The total number of potential flow lines, M, was arbitrarily selected as 
127 (a maximum of 148 may be used), while 45 total streamlines (N) were used 
(maximum 48). MIDJ, the number of the ‘‘center streamline," was 21. MIDJ and 
the Aip values in the vicinity of this streamline may have to be adjusted from 
the first trial values to prevent the design from becoming extremely asymmetrical, 
possibly crossing the X-axis. This situation will result. in a divide check, 
an overflow, or an underflow, with subsequent program termination. 

Dirichlet boundary conditions were used at the flow channel inlet (not 
to be confused with the diffuser inlet) with a constant velocity magnitude 
ratio, Q, (not parallel) across the channel of 1.0 being specified. At the 
flow channel exit (which was also the d i f fuser exi t) , Neumann boundary condi- 
tions were specified to yield parallel (not necessarily uniform) flow across 
the exi t plane. 

The next two pages show the input data which was used to design the 
18-inch (46 cm) annular diffuser. Following the input data is the first 
page of output which summa r i zes -The control input data. Then the graphical 

output {figure 22) is shown reduced in size from the original output. Also 
shown is a listing of inner and outer wall radii. 
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THE FOLLOWING DAT A WAS INPUT TO DESIGN THE 18-INCHAN NULAR DIFFUSER 


2 1 8 1 G ANNULAR DIFFUSER DESIGN 070171 0011 
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CONTROL OATA 
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FLAG'! =0 TOLSYS = 0.000050 MXISYS = 10 IP3SYS = 11 

FLAG4 =1 TOLOl = 0.000050 MXIDE 5 50 FPROE = 51 

MSTAGU = 84 NSTAGU = 39 MSLOTO = 87 

MSTAGL * 84 NSTAGL = 6 MSLOTL = ~ST " 


OPTIONS 


AXIALLY SYMMETRIC 

DIRtCHLET B. C. AT INLET AMO NEUMANN 8. C. AT OUTLET 
WITH VELOCITY INPUT AT INLET 

ANNULAR DIFFUSER - NOZZLE 

GRAPHICAL OUTPUT 


o 

o 



3.20 U.OO 14 . ao 

X fix IS l*N INCHES) 


Figure 22. Graphical output showing streamlines and equi potent i a 1 lines and wall coordinates for 
18-inch (k6 cm) annular diffuser. 
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0.21 5395E 

01 

0. 152148E 00 


0 • 1058716 

01 

0.713184E 

01 

0.322997E 

01 

0.21 8497E 

01 

0.169524E 00 


0 .1 08951 E 

01 

0.713740E 

01 

0.323997E 

01 

0.2215886 

01 

o.iaeoiiE oo 


0.1 12012E 

01 

0. 7 1435 IE 

01 

0. 324997E 

01 

0.2246706 

01 

0. 207 168 E 00 


0 • l 15053E 

01 

0* 7 1502 1 E 

01 

0.325998E 

0! 

0. 227742E 

01 

0.227265E 00 


0.1 180736 

01 

0. 715751E 

01 

0.326997E 

01 

0 . 230804E 

01 

0.248460E 00 


0. 121068E' 

01 

0* 716545E 

01 

0.327997E 

01 

0.233858E 

01 

0.2708 726 00 


_? . 1 24036E 

01 

0.717405E 

01 

0.328997E 

01 

0.2369026 

01 

0 • 294692E 00 


0. 126977F 

01 

0.7183356 

01 

0.329997E 

01 

0.239937E 

01 

0.320056F 00 


_0_. 1 2988 7E_ 

01 

0.719339E 

01 

0.330997E 

01 

0.242962E 

01 

0. 3472096 00 


0. 132758E 

01 

0.72042 IE 

oi 

0.331997E 

01 

0.245979E 

01 

0.37 6464 E [00 


0.135591E 

01 

0.7215876 

01 

0 . 332997E 

01 

0. 248986E 

01 

0.408066E 00 


0. 1 38377E 

01 

0.722942E 

01 

0.3 33997E 

01 

0 . 2 5 1 9 8 5 E 

01 

0.4423986 _00 


0 • 1411 HE 

01 

0. 7241 9 IE 

01 

0.3 34997E 

01 

0. 254974E 

01 

0.479924E 00 


0. 143785E 

01 

0. 725643E 

01 

0 . 3 35997E 

01 

0.257955E 

01 

0.521236E 00 


0. 14638BE 

_0l_.. 

0. 727204E 

01 

0 • 3 36998E 

01 

0. 260926E 

01 

0.567117E 00 


0 • l 48907E 

J)l„ 

0. 7288846 

01 

0.3379978 

01 

0.2638896 

01 

0.618589E 00 


0.1513256 

01 

0. 730691E 

0,1 

0.338997E 

01 

0 • 266844E 

01 

0. 677078E 00 


0. 153619E 

01 

0. 7 326 37 E 

01 

0* 339997E 

01 

0.269789E 

01 

0.7445536 00 


0 .1557606 

01 

0.734732E 

01 

0. 34099 7E 

01 

0.272726E 

01 

0.8239586 00 


0.157704E 

01 

0^736986E 

01 

0.341998E 

01 

0 • 275654E 

01 

0.9197446 00 


0. 159392E 

01 __ 

0.7394 12 E 

01 

0.342996E 

01 

0 • 278574E 

01 

0. 1039QQE 01 


0.160736E 

01 

0.742018E 

01 

0. 343996E 

01 

0.2814856 

01 

0.1193616 01 


0. 161615E 

01 

0. 744842E 

01 

0.344746E 

01 

0 • 204389E 

01 

0.140277E 01 


0. 161876E 

01 

0.7 4809 6E 

01 

0. 345747E 

01 

0.2072856 

01 

0.169840E 01 


0.1615O7E 

01 

0. 7 52623E 

01 

0.346746E 

01 

0.290173E 

01 

0.213280E 01 


Q , t 1605 30E 

01 

. 0. 753676E 

pi 

0.347246E 

oi 

0.2916146 

01 

0.2 50805E 01 
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C.159233E 

01 

0.754251E 

01 

0. 347746E 

01 

0.293053E 

01 

3.2943706 

01 

0. 15781 LE 

01 

0.754202E 

01 

0.348246E 

01 

0 .294490E 

01 

3.341566E 

01 

' 0 .1 565 73 E . 

01 

0. 7535 08E 

01 

0. 34874 66 

01 

0.295924E 

01 

0. 389990f_ 

01 

0.155819E 

01 

0.752339E 

01 

0.349246E 

01 

0.297357F 

01 

3 . 437780E^ 

01 

SUMf VRY OF 

DOWNSTREAM wall coordinates 

OF DIFFUSER 

1UPPER 

WALL ). 


.... 


X 


R 


0 


S 


ALPHA 



0.167012E 

01 

0. 765309E 

01 

0. 1000C0E 

01 

0.0 


0.285718 E 

01 

0.1 70109E 

01 

0 . 7703 1 1 E 

01 

0. 830000E 

00 

0.546426E- 

-01 

0.1 1 1739E 

01 

0. 173 1938 

01 

0.775789E 

01 

' 0.8300006 

00 

0. 1 1 4ft 83E 

00 

0.10L841E 

01 

y;i'7 6 7 0 ^“ 

Oi 

O'. 7 809 3 5E 

01 

0. 8 30000E 

00 

0.175124E 

00 

0.9378 776" 

00 

0. 180540E 

01 

0. 785774E 

01 

0. 830000b 

00 

0.235365E 

00 

3.871287E 

oo 

0. 184630E 

01 

0.790338E 

01 

0. 830000E 

00 

0.295606E 

00 

3.8 148936 

do 

0 . 1889 1 OF 

01 

0. 794649E 

01 

0. 830000E 

00 

0.355847E 

00 

0.766357E 

00 

0. 1 97973 E 

01 

0.0026326 

01 

0. 830000E 

00 

0.476328E 

00 

3.685222F 

1)0 

0.207568E 

01 

0.8 099 20 E 

01 

0.8 30000E 

00 

0.596810E 

00 

0.6192816 

bo’ 

"5. 2175546 

01 

0.8 L6603E 

01 

0. 830000E 

00 

0 • 71 7291 E 

00 

3 • 564065E 

00 

0.22 7 840E 

01 

0.822757E 

01 

0.8 30000E 

00 

0.837773E 

00 

0.516799E 

00 

0 .243722E 

01 

0,8311176 

01 

0. 830000E 

00 

0. 101850E 

01 

0.456696E 

ob 

0 .260025E 

01 

0.838592E 

01 

0. 830000E 

00 

0. 119922E 

01 

0.4063 5 6E~ 

~bb 

0. 276643E 

01 

0. 845301 E 

01 

0. 830000E 

00 

0. 137Q94E 

01 

0.363449E 

oo 

~ 0~, 293 509 E 

01 

0. 851342 E 

01 

0. 830000E 

00 

0. 156066E 

01 

3. 326374E 

bo 

0.3 105706 

oi"' 

0. 856796E 

01 

0. 830000E 

00 

0.17413BE 

01 

3.293998E 

00 

0.3277886 

01 

0.861730E 

01 

0.8 30000E 

00 

0.192210E 

01 

3. 265491 E 

oo 

0.3451356 

Ol" 

0.866203E 

01 

0. 830000E 

00 

0.210283E 

01 

0. 240230 E~ 

"00 

0.362 5 87 E 

01 

0.8702626 

Oi 

0 . 8 30000E 

00 

0.228355E 

01 

0.217'724T 

00 

0.3801 26E 

oi 

0.873950E 

01 

0. 830000E 

00 

0.2464276 

01 

0. 1 9 7582 E 

00 

6 . 3 97 7 39 E 

oi 

0.877304E 

01 

0.830000E 

00 

0.264499E 

01 

3 • 1 79487E 

00 

0.4154136 

01 

0. 83035 7E 

01 

Q.830000E 

00 

0.282571E 

“bi 

3. 16 31816 

ob'”' 

0 .4331406 

oi 

0.8831 36E 

01 

0. 8 30000E 

00 

0.3006446 

01 

3. 148448E 

00 

0. 45091 OE 

01 

0.8856696 

01 

0.8 30000E 

00 

0.31871 6E 

01 

0. 135097E 

ob 

0 .468719E 

01 

0. 887976E 

01 

0.8 30000E 

00 

0. 336788E 

01 

3. 1229736' 

00 

0. 486561 E 

oi " 

0.8900776 

01 

0. 830000E 

00 

0 . 3548606 

01 

3. 111936E 

"oo 

0 • 5044 3 OF 

oi 

0.89 1 9 9 2 E 

01 

0.830000E 

00 

0 * 372932E 

01 

0. 1 01 0706 

bo 

0 75 22323'E 

"oT~ 

0.893736F 

01 

0.8 30000E 

00 

0.391004F 

01 

0. 926709 E"- 

-01 

0.540238E 

01 

0.895322E 

01 

0. 830000E 

00 

0.409076E 

01 

3 . 842395E- 

-01 

0.5581 7 OF. 

oi 

0.896764E 

01 

0. 8300C0E 

30 

0.427149E 

01 

^• 764982 E- 

-01 

0 .5761 19F 

01 

0. 898073E 

oi 

0. 830000E 

00 

0 . 44522 l E 

01 

3.69 3699?^ 

-01 

0 .5940 80E 

01 

0. 89926 IE 

01 

0. 8 30000E 

00 

0.463293E 

01 

0. 62 7882E- 

-01 

0.6^20556 

“oi" 

0.9003 34E 

01 

0.8 30000E 

00 

0.481365E 

01 

0. 5669006- 

-01 
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0.6300406 

01 

0.901302E 

6! 

0. 8 30000E 

00 

0.499438E 

01 

0. 510197E-01 

0.6480 33 E 

01 

0.902172E 

01 

0.8 300COE 

00 

0.5175L0E 

01 

0.4573536-01 

0.666034E 

01 

0.9029 5 l E 

01 

0. 630000E 

00 

0.535582E 

01 

0.4078266-01 

0.6040466 

01 

0.903642E 

01 

0*8 30000E 

00 

0.553654E 

01 

0. 3609346-01 

0.7020636 

01 

0.9042536 

01 

0*8 30000E 

•00 

0 . 57 1 7266 

01 

0. 316626E-01 

0.7200856 

01 

0.9047856 

01 

0* ft 300 00 E 

00 

0. 589798 E 

01 

0.2744366-01 

0.7381126 

01 

0.905242E 

01 

0.8300006 

oo 

0 . 607871E 

01 

0.2340456-01 

0.7561426 

01 

0.9056 30E 

01 

0.8300006 

00 

0.625943E 

01 

0. 195101 E-01 

0.7741746 

01 

0. 90 5946 E 

01 

0*8300006 

00 

0 .6440 1 *>E 

01 

3. 1573986-01 

0.7922136 

01 

0. 906 1 98 E 

01 

0.8300006 

00 

0.662087E 

01 

0. 1 2 0344E-0 1 

0.8 100676 

01 

0. 9 0638 0 E 

01 

0. 830000E 

00 

0.6801596 

01 

0.974582E-02 




























^SUMMARY OF 

”wauT 

COORDINATES INSIOE OF 01 

:ffuser SLOT (UPPER 

WALL) 



X_ 


' R 


0 


S 


ALPHA 










0. 1670126 

61 

0. 7658Q9E 

0! 

0.0 


0.0 


0.2857186 01 

0.1553536 

~oT 

0. 7631 34E 

01 

0.830000E 

00 

0. 1 1 96 16E 

00 

0. 3418486 01 

0,Tl 50 1786 

01 

0. 7602 1 76 

01 

0.830000E 

00 

0.1798566 

00 

0.3903856 01 

0.1469566 

01 

0.7549686 

01 

0.830000E 

00 

0.2400976 

00 

D • 444 i 856 01 




















SUMMARY OF INLET WALL COORO I NATES OF DIFFUSER (LOWER WALL) 


X R 0 S ALPHA 


. 

o 

* 

o 

0.391000E 

01 

0. 100000E 

01 

0.0 


0.2249476 

01 

-0.52 26146-01 

0.3961676 

01 

0. 128412E 

01 

0.8801816 

-01 

D.2Tr9~46F 

01 

-0 .8883486 -0 l 

0.401695E 

01 

0. 149755E 

01 

0.159424E 

00 

0.1989596 

01 

-0.1136826 00 

0.406892E 

01 

0.1653816 

01 

0.2229236 

00 

0.1867746 

01 

-0 • 1 304 19 E 00 

0.4118686 

01 

0. 179784E 

Oi 

0 • 2 80845 E 

00 

3.1757986 

01 

-0.140355E 00 

0.416647E 

01 

0.191427E 

01 

0.3347266 

00 

0 . 165572F 

01 

-0.1449366 00 

0.421239E 

01 

0.202216E 

01 

0.3855226 

00 

0.1561146 

"on ■ “ 

-0.144800E 00 

0.4256546 

01 

0.2108136 

01 

0.4339436 

00 

0. 1473626 

01 

-0.1408416 00 

0.4299006 

01 

0.218835E 

01 

0.4804856 

00 

3.1394546" 

01 

-0.1336106 00 

0.4339756 

oi 

0.2263066 

01 

0.525411E 

00 

3. 1320796 

01 

-0.1235186 00 

0. 4378 79 E 

01 

0. 233251 E 

01 

0.568925E 

00 

3.l2 5ll4f 

~ 0 \ “ 

-0.1109306 00 

0.4416S9E 

Of 

0.239693E 

01 

0.6111966 

00 

0.1184456 

in — 

-C. 9602256-01 

0. 445 1 76 E 

'61 

0.244696E 

01 

0.6524656 

00 

• 3 • H2266F 

in — 

-0.7924 52 ErOl 

0. 44858 86 

of 

0.250280E 

01 

0 .6928566 

00 

0 :iO 6459F 

ot 

-0.6068576-01 

0.4518406 

oi 

0.2546026 

01 

0.732453E 

00 

0 . 1008 796" 

oi 
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-0.4057916 

rOl 

0. 4549476 

01 

0. 25863 76 

01 

0.771419E 

00 

0.9572636 00 

-0.191262E 

-01 

0.457912E 

01 

C. 262400E 

01 

0. 8098 01 E 

00 

0.90 87 996“ W 

0. 35259 IE 

-02 

0-460741E 

01 

0. 2659C3E 

01 

0.847655E 

00 

0.863020E 00 

0.2723746 

-01 

0.4634 3 7E 

01 

0.269126E 

01 

0.885033E 

oo 

0.8197726 00 

0.5189106- 

-01 

0.4660Q7E 

01 

0.272190E 

01 

0.92197BE 

00 

0.7787626 00 

0.773776E 

-01 

0. 468454E 

01 

C. 27500 IE 

01 

0.9585256 

00 

D. 7 39717E 0~0 

0. 103615E 

00 

6'. 4 7 0784 E 

01 

0.27761 IE 

01 

0.9947 1 5E 

00 

0 • 702566E 00 

0.130536E 

00 

0.47 300 1 E 

01 

0.2300306 

01 

0.1 O3O50E 

01 

D.666868E 00 

0.1580776 

00 

0.4751086 

01 

0.2819416 

01 

0.106618E 

01 

0 • 633534E 00 

0.1861446 

00 

0.47711 7E 

01 

0. 283957E 

01 

0.110153E 

01 

0 . 602375E 00 

0 .’2146806 

•bo ‘ 

0. 4790266 

01 

0.2859256 

01 

0.11 3662E 

01 

0.571968E 00 

6 • 243651 6 

00 

0. 480840 E 

01 

0.2874036 

01 

0.1171496 

01 

0.543395E 60 r 

0 • 2 73003E 

00 

0. 48 2 567 E 

01 

0.2891366 

01 

0. 12061 86 

01 

3.5165726 00 

0 • 302662 E 

00 

0. 484207E 

01 

0.2907636 

01 

0.124067F 

01 

0.490305E 00 

0.3326^6 

c'o” 

0.485764E 

01 

0.291992E 

01 

0.127499E 

01 

0.4655046 00 

0.362884E 

00 

0.487243E 

01 

0. 293456E 

01 

0. 130915E 

01 

3.4419226 00 

0 ."3933 51 E 

ob ' 

“ .0.48864 5 E 

01 

0.294576E 

01 

0. 134317E 

01 

0.419361E 00 

0 .4240396 

oo"” 

0.4899 776 

01 

0.295922E 

01 

0.1 37704E 

01 

3.3978316 00 

0.454915E 

00 

0.4912396 

01 

0.2969676 

01 

0.1 41077E 

01 

0.3771586 00 

*0.4 8592 96 

00 

0.492434E 

01 

0. 2902436 

01 

0. 144438E 

01 

0.3572976 00 

0. 5 170 78 E 

00 

0.4935646 

01 

0.2992496 

OI 

0.I47785E 

OI 

0.3381576 00 

0 • 548330E 

00 

0. 494633E 

01 

0. 300499E 

01 

0.151 1 2CE 

01 

3.3196486 00 

0 • 5796866 

ob 

0.4956406 

01 

0.301498E 

01 

0.1544436 

01 

0.30 16896 00 

b.^TioVbF 

“oo 

0.4965896 

oi 

0. 302749E 

01 

0. 1577536 

01 

3.2841396 00 

0.642576E 

00 

0. 497480E 

01 

0.303748E 

01 

0-1610506 

01 

0.2669416 00 

6.67^ 10 3E 

00 

0*4983i 5E 

01 

0. 3049996 

01 

0.164336E 

01 

0.2500186 00 

0.7056726 

oo 

0. 499093E 

01 

0. 305999E 

01 

0. 167609E 

01 

0.2331356 00 

0 • 7372 71 E 

00 

0.499817E 

oi 

0.306998E 

01 

0.1708726 

01 

0.2167646 00 

0.7688736 

00 

0. 500487E 

01 

0* 307998E 

01 

0. 174124E 

01 

■ 3 • 230804E. 00 

0.800478E 

00 

0. 50 1 1 05E 

01 

0.309246E 

01 

0.177364E 

01 

' 0. 1846696 00 

6. 8 320 75 E' 

00 

0. 501 669 E 

01 

. 0.310248E 

01 

0* 1 80593E. 

01 

0.1682876 OO 

0 • 663653E 

00 

0. 502 l 806 

or 

0.3112496 

01 

0 • 18381 16 

01 

3 • 15 22906 00 

0.8952066 

do 

0. 5026 39E 

01 

0.312498E 

01 

0 . 107O17E 

01 

0. 135944E 00 

0.926719E 

“oo 

“"0.50304 36 

01 

0.3134996 

01 

Oi 1902116 

01 

3.1191516 00 

0.9581 92 E 

"00 

0.5033946 

01" * 

0.3144986 

01 

0.193396E 

01 

0. 102403E 00 

0.989616E, 

oo 

6.50369 OE 

01 

0.315499E 

01 

0. 196571E 

01 

0.85 6250E-01 

0. 102098E 

01 

0. 503933E 

01 

0.316748E 

01 

0 * 1997 346 

01 

0.6021256-01 

0. 105227E 

01 

0.5041196 

oi 

0.317747E 

01 

0. 2020866 

01 

0. 5OO063E-O 1 

0.108348 E 

oi 

0.5042476 

01 

0. 3 18748E 

01 

0* 206028E 

01 

3.316566E-01 

0 .ill 46 1 E 

01 

6"'56'43i7E 

01 

0 • 3 1974 76 

01 

0.2091606 

01 

0. 1274706^1 

0.1 14564E 

01 

0.5043276 

01 

0.3207466 

01 

0.2 12282E 

01 

-0.67 439T6^b? 

0. 117657E 

oi 

*1)750421766 

01 

0.3217486 

01 

0.2153956 

01 

-3.2676456-01 

0.1 20737E 

61 

7 0.50 4 1626 

01 

0.3229976 

01 

0.2184976 

01 

-0. 4792496-01 










— - 


_ 

-------- 

- - - " 

- - - 






0.123800E 

01 

0.50398 IE 

01 

0.3 2399 7E 

01 

O.221508E 

01 

-0.703674E- 

-01 

0. 126848E 

01 

0. 50373 1 E 

01 

0.324997E 

01 

0.2246 TOE 

01 

-0 .9364776- 

-01 

0.129879E 

oi 

0.503410E 

01 

0.325998E 

01 

0.227742E 

01 

-0. 1 l 8046E 

Ob — 

0.132890! 

01 

0.5030L5E 

01 

0.326997E 

01 

0 • 2 30804E 

01 

-0.143745T 

Oo 

0 3 58 78E 

01 

0. 50 2 54 2 E 

01 

0.327997E 

01 

0 . 233858E 

01 

-0. 17090 0 E 

bo 

" b v 138830E 

01 

0.501987E 

01 

0.328997E 

01 

0 . 236902E 

01 

-0. 199717E 

00 

0 « 14 1 7 71 E 

01 

0.5013476 

01 

0.329997E 

01 

0.239937E 

01 

-0.2303626' 

00 

1 7T4 4670E 

01 

-0.5006 rsr oi 

073309976 

01 

0 . 242962E 

01 

-0.2631 06TT 

00 

0 .147 5 30 E 

oi 

0.499794E 

01 

0.331997E 

01 

0.245979E 

01 

-0.2982346 

ob 

0.1 50346E 

01 

0.4988 70 E 

01 

0.332997E 

01 

0.248986E 

01 

-0. 3360496 

oo 

0. 153U1E 

01 

0. 497 839 E 

01 
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0.1617436 

01 
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0.830000? 

bo 

0.414466E 

oo 

-0 • 3957416 01 



APPENDIX C. DEVELOPMENT OF PERFORMANCE EQUATIONS 

The following assumptions were made in the derivation of the per- 
formance equations: 

(1) Tests were conducted under steady state steady flow conditions. 

(2) The temperature (and thus the internal energy) and the pressure 
were uniform across the inlet section and across the exit section. 

(3) The fluid density was constant. Since the maximum velocities 
encountered resulted in a Mach number of about 0.25, the resultant error 
should be negligible. 

Define the effectiveness c as follows: 


e = 


m (P - P .)/p 
e s ,e s , i 

m. (KE. - KE ) 
i i e 


(27) 


where m represents the mass flow rate and KE is the kinetic energy per unit 
mass flowing. 


r u . 2 

m. KE. = 4- (U. p) dA. (28) 

I I sJ A Z I I 

• ' ' ii 2 

( m. x ! m. . r u 

T) ™e KE e - \t) J -f- (U e^ )dA e <»> 


e A 


E = 


m (P - P .)/p 
e x s,e s, i 


i 

A. 

. i 


U.‘ 


2 vu i 


(U. p) dA 


- 0 ^ 


(30) 


(U p) dA 


e 


e A 


Now define a kinetic energy coefficient k: 


108 




so that the effectiveness expression becomes 



This was the expression used to compute the effectiveness e. The more 
conventional expression for the effectiveness is similar to equation 25 except 
that the ratio of mass flow rates is deleted. 

Thus 



Considering the total pressure loss, 




( 36 ) 



Values of k. (and k ) were computed from measured velocities as follows, 
i e 

Consider n velocities IL at n radii r^ incremented by a constant interval 
(where j = 1 , 2, . . . , n) . If we let the velocity IK be representative of the 
annular area between radii 



and 



then 


4 S&: 


A U 


I 


n U. 3 o 

dA = j 1 (z4 J ( r :+l 

j=l U 4 J 


r J-l + 2 r j (r J + l 


r M )] 


(37) 


where 


_ 1 n 2 
u s u j l 'U 


r j-i + 2r j (r j+ i 




(38) 


In equation (38) the radii r Q and r n+ | are equal to the inner and outer diffuser 

wall radii respectively. All other radii, rj to r n> are readily obtained since 

the increment between radii is a constant. Using equation (37), typical values 

obtained were 1.025 for K. and 1.011 for K . 

i e 


no 


E-7242 



